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1.0   INTRODUCTION 
1 . 1  Background 

The  term  "fluidized  bed",  as  it  is  now  used,  relates  to  a  technique  initi- 
ated more  than  twenty  years  ago  to  find  a  "better  catalytic- cracking  process  than 
the  fixed  bed  method  (lk) . 

A  fluidized  bed  is  a  collection  of  solid  particles  which  circulate  in  an 
enclosed  region  as  the  result  of  the  upward  flow  of  fluid,  either  liquid  or 
gas,  past  the  particles.  A  fluidized  state  occurs  over  a  finite  range  of  fluid 
velocities.  At  zero  fluid  velocity,  all  of  the  particles  are  in  a  static 
state  resting  on  the  bottom  of  the  containment  vessel.  As  the  fluid  velocity 
is  increased,  there  is  a  certain  rate  of  fluid  flow  at  which  the  bed  will  be 
expanded  to  such  a  point  that  the  particles  will  begin  to  move.   The  static  bed 
is  characterized  by  a  pressure  drop  across  the  bed  that  increases  with  in- 
creased fluid  velocity.  When  a  fluid  velocity  is  reached  that  causes  particle 
movement,  only  a  slight  change  in  the  pressure  drop  occurs  across  the  bed  as 
the  fluid  velocity  is  increased.  Further  increases  in  the  velocity  cause  the 
volume  of  the  bed  to  increase. 

At  high  fluid  velocities,  a  bed  fluidized  with  gas  usually  remains  largely 
in  a  compact  state.  A  bed  fluidized  with  a  liquid  expands  steadily  with  in- 
creasing fluid  velocity  until  a  state  of  complete  particle  isolation  is  reached. 
Maximum  expanded  bed  heights  for  gases  rarely  exceed  twice  the  static  bed 
height . 

At  low  fluid  velocities  slightly  above  those  for  minimum  fluidization,  the 
whole  bed  is  macroscopically  homogeneous.   This  state  is  sometimes  referred  to 
as  the  state  of  incipient  fluidization.   In  describing  the  physical  appearance 
of  the  bed,  this  form  of  flow  has  been  referred  to  as  particulate  fluidization. 


At  higher  fluid  velocities,  a  second  phase,  commonly  known  as  the  boiling 
bed  state,  often  occurs.   Bubble  formation  and  slugging  take  place  in  the 
fluidized  solids.   The  physical  appearance  of  the  bed  is  then  referred  to  as 
aggregative  fluidization. 

Whether  particulate  or  aggregative  fluidization  occurs  at  high  fluid 
velocities  relative  to  the  minimum  fluidization  velocity  depends  upon  the 
particle  size  and  the  ratio  of  solid  to  fluid  density  (7).   In  general,  the 
particulate  form  occurs  in  liquid  fluidized  systems,  and  aggregative  fluidi- 
zation occurs  in  gas  fluidized  systems. 

Particle  mobility  is  the  primary  difference  between  a  fluidized  bed  and 
a  fixed  bed.  Many  of  the  advantages  and  disadvantages  of  the  fluidized  bed 
technique  have  been  enumerated  in  the  literature  (7;  9,   10,  1*0. 

1.1.1  Advantages 

The  following  items  are  some  of  the  major  advantages  of  applications  of 
the  fluidized  bed  technique. 

1.  Turbulence  within  the  fluidized  system  agitates  and  diffuses  any  cold 
or  hot  spots,  providing  an  approach  to  isothermal  operating  conditions. 

2.  Turbulence  provides  a  more  uniform  solids  distribution. 

3.  A  fluidized  bed's  ability  to  handle  solids  as  if  they  were  a  liquid 
enables  continuous  operations,  such  as  addition  and  withdrawl  of  solids  from 
the  bed,  to  occur. 

k.      High  volumetric  heat  transfer  rates  can  be  conveniently  accomplished 
due  to  the  large  surface  area  available  for  heat  transfer. 

5.   High  heat  transfer  coefficients  from  the  bed  to  the  walls  of  the  con- 
taining vessel  are  possible  due  to  the  presence  of  the  fluidized  solids. 


6.  At  the  same  fluid  flow  rate,  smaller  pressure  drops  usually  occur  in 
a  fluidized  bed  with  no  restraints  than  in  a  packed  "bed. 

1.1.2  Disadvantages 

The  following  items  are  some  of  the  disadvantages  of  fluidized  "beds. 

1.  Only  solids  which  flow  freely  and  do  not  tend  to  agglomerate  may  be 
used. 

2.  For  certain  processes  a  temperature  gradient  may  be  desirable,  and 
this  cannot  be  accomplished  utilizing  the  fluidized  bed  technique  unless  stag- 
ing, partitions,  or  internal  heat  exchangers  are  incorporated. 

3-  Erosion  of  the  containing  vessel  and  internal  solids  may  occur. 
k.      Bubble  formations  rising  through  fluidized  solids  may  reduce  the 
available  contact  surfaces,  reducing  the  efficiency  of  the  process. 

1.2  Application  to  Nuclear  Reactors 

Fluidized  bed  techniques  have  been  applied  with  considerable  success  to 
many  processes  where  heat  removal  and  temperature  control  are  required.   These 
factors  are  important  in  a  nuclear  reactor,  and  several  proposals  of  the  appli- 
cation of  fluidized  bed  techniques  to  nuclear  reactors  have  been  submitted 
(8,  9>  10,  1^) •  Advantages  of  a  fluidized  bed  nuclear  reactor  design  include 
uniform  temperature  distribution,  uniform  burnup  independent  of  power  distri- 
bution, ease  of  fission  product  removal,  ease  of  fuel  fabrication  and  repro- 
cessing, and  low  fuel  holdup  (7,  10,  Ik) . 

Neutron  moderation  and  heat  removal  could  be  accomplished  by  the  fluid  in 
a  liquid  fluidized  bed  nuclear  reactor.  Heat  removal  would  be  limited  by  the 
fluid  velocity  required  to  maintain  fluidization  under  this  arrangement. 


Ordinary  control  rods  could  be  inserted  directly  into  the  fissionable  fluidized 
solids,  or  the  control  could  be  self- regulating.  For  self- regulation  the  bed 
could  be  designed  so  that  the  critical  mass  would  be  independent  of  reasonable 
fluctuations  in  the  bed  height,  but  a  scram  could  be  accomplished  by  shutting 
off  the  fluid  flow  (10).   Proposals  for  fluidized  bed  nuclear  reactors  should 
meet  the  following  requirements  (9): 

1.  The  fluid  must  have  enough  velocity  to  support  the  bed. 

2.  Heat  must  be  removed. 

3-  A  knowledge  of  the  expansions  and  contractions  of  the  bed  and  concen- 
trations of  the  solids  is  necessary. 

k.     Pressure  drops  across  the  bed  will  have  to  be  considered. 
5.  Erosion  due  to  particle  mobility  must  be  considered. 

1.3  Purpose 

Although  fluidized  bed  techniques  have  been  applied  with  great  success 
in  many  industrial  processes,  little  knowledge  has  been  gathered  of  the  actual 
movement  of  particles  in  fluidized  beds.   Before  it  is  possible  to  obtain  a 
basic  understanding  of  heat  and  mass  transfer  within  a  fluidized  bed  reactor, 
it  will  be  necessary  to  have  a  fundamental  understanding  of  the  motion  of  the 
particles  in  the  fluidized  bed.  Knowledge  of  the  degree  of  backmixing  and 
relative  velocities  of  the  particles  and  fluid,  along  with  other  character- 
istics of  the  fluidized  bed,  are  fundamental  requirements  to  an  understanding 
of  the  mechanisms  of  heat  and  mass  transfer. 

The  literature  reveals  hundreds  of  articles  pertaining  to  fluidized  beds, 
but  only  Gordon  (k)   and  Wilde  (13)  have  provided  information  on  the  statisti- 
cal behavior  of  individual  particles  within  the  bed.   It  would  seem  appropriate 


that  experimental  research  work  "be  carried  out  to  further  study  individual 
particle  motility.   Thus,  it  has  "been  the  purpose  of  this  thesis  to  study 
experimentally  the  movement  of  individual  particles  in  a  fluidized  bed. 

The  solids  and  fluid  used  in  the  study  were  stainless  steel  "ball 
bearings  and  water. 

l.k     Scope 

Water  velocities  ranging  from  minimum  fluidization  velocity  to  the 
velocity  required  to  expand  the  fluidized  bed  to  twice  its  static  height 
were  used.   One  ball  bearing  size  and  one  static  bed  height  of  ball  bearings 
were  utilized  in  the  experimental  study.  Particle  concentrations  and  mean 
vertical  speeds  were  obtained.   The  development  and  construction  of  the 
experimental  apparatus,  the  method  of  obtaining  the  experimental  data,  and 
a  statistical  analysis  approach  of  evaluating  the  experimental  data  are  in- 
cluded in  this  thesis.   Tracing  the  movement  of  a  radioactive  ball  bearing 
in  the  fluidized  bed  was  used  as  the  method  of  studying  individual  particle 
mobility. 


2.0  EXPERIMENTATION 
2.1  Experimental  Apparatus 

A  photograph  of  the  experimental  apparatus  is  shown  in  Fig.  1,  and  a 
schematic  diagram  is  presented  in  Fig.  2. 

The  apparatus  may  be  divided  into  the  following  categories:   (l)  the 
column  and  its  accessories,  (2)  the  gamma  ray  detection  equipment,  and  (3) 
the  measuring  and  recording  equipment. 

2.1.1  The  Column  and  Its  Accessories 

A  kS- inch-high,  3-0- inch- inside  diameter,  and  3- 5- inch- out side  diameter 
cast  acrylic  resin  tube  was  used  as  the  fluidizing  column.   It  was  fitted  with 
a  flange  on  each  end.  Figure  3  provides  a  close-up  view  of  the  column  and  bed 
of  solids.   The  outlet  of  the  column  was  fitted  with  a  wire  screen  to  prevent 
any  of  the  fluidized  solids  from  accidentally  escaping  the  column.   Tne  two 
flanges  were  fitted  with  pressure  taps  for  measuring  the  pressure  drop  over 
the  bed  of  solids.   The  lower  flange  was  bolted  to  the  distributor  assembly. 

Tne  distributor  assembly,  shown  in  Fig.  k,    consisted  of  a  calming  section 
of  3/8- inch  ceramic  balls,  housed  in  a  container  eight  inches  high  and  three 
inches  in  diameter.   A  brass  sieve  separated  the  calming  section  from  the 
fluidized  particles.   The  holes  in  the  brass  plate  were  3/32-inch  diameter, 
and  the  ratio  of  hole  area  to  total  area  of  the  sieve  was  0.3^3- 

As  shown  in  Fig.  2,  a  30-inch  mercury  manometer  was  used  to  measure  the 
pressure  difference  over  the  orifice  meter. 

The  water  flow  through  the  fluidized  bed  was  provided  by  a  closed  system. 
A  55-gs-llon  drum  coated  with  chr ornate  primer  to  prevent  corrosion  served  as 
the  water  reservoir. 


Figure     I.    The     Column     and     its     Accessories. 
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Figure    3.    Solids    Within    the    Fluidizing    Column, 
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A  Type  S.C.,  U.S.  Motor,  rated  at  5  horsepower,  220  volts,  three-phase, 
60  cycles,  and  1800  revolutions  per  minute,  powered  the  pump.   The  pump  was 
a  centrifugal  type,  with  a  2.5-inch  inlet  and  a  2.0-inch  outlet.   It  was  "built 
by  Fairbanks-Morse. 

The  structure  used  to  support  the  column  and  its  accessories  was  con- 
structed from  Uni strut. 

The  water  piping  system  was  constructed  of  standard  schedule  40  1^-inch 
diameter  galvanized  pipe. 

2.1.2  The  Gamma- Ray  Detection  Equipment 

Four  scintillation  probes  were  used  to  detect  the  radioactive  stainless 
steel  ball  in  the  fluidized  bed.   Two  probes  were  used  at  each  of  two  horizontal 
planes.   The  horizontal  planes  were  constructed  in  the  following  manner:   Lead 
bricks  of  2xUx8-inch  size  were  stacked  around  the  column,  with  l/8-inch  lead 
shims  placed  between  the  layers  of  bricks  at  each  of  the  two  horizontal  planes 
studied.   The  increased  output  from  the  scintillation  probes  behind  each  slit 
provided  a  means  of  detecting  the  movement  of  the  particle  each  time  the  radio- 
isotope passed  one  of  the  horizontal  window  slits  in  the  shielding. 

Figure  5  is  a  photograph  of  the  scintillation  probes  positioned  behind 
the  lead  shielding.  Part  of  the  lead  bricks  are  removed  from  the  top  of  the 
shield. 

Figure  6  illustrates  the  vertical  position  of  each  scintillation  detector 
level  and  the  distance  between  the  two  detector  levels,  A  Z. 

The  four  scintillation  probes  used  were  B.J.  Electronics  of  Borg- Warner 
Corporation,  Type  "JDk   Sodium  Iodide  (Thallium)  crystal,  Model  DP3,  serial  number 
B-168,  and  NE  numbers  287,  288,  3lU,  and  353. 
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Figure    6.     Schematic    of     Scintillation     Probes     and     Shielding. 
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2.1.3  The  Measuring  and  Recording  Equipment 

Figure  7  is  a  schematic  diagram  of  the  detection  and  recording  equipment. 

High  voltage  pover  supplies.   The  photomultipliers  of  the  two  scintillation 
probes  at  each  of  the  two  detector  levels  were  supplied  with  high  voltage  from 
B.J.  Electronics,  Borg-Warner  Corporation,  Model  DVI,  Detectolab,  NE  262  and 
NE  275  variable  high  voltage  power  supplies.  NE  275,  operated  at  about  10 50 
volts,  supplied  high  voltage  to  the  upper  detector  level,  and  NE  262,  operated 
at  about  1100  volts,  supplied  high  voltage  to  the  lower  detector  level.   The 
two  power  supplies  are  located  at  the  top  of  the  instrument  rack  in  Fig.  8. 

Ratemeters.  The  output  from  each  detector  level  was  fed  to  a  ratemeter, 
as  shown  in  Fig.  7*  The  two  ratemeters  are  shown  side-by- side  directly  below 
the  high  voltage  power  supplies  in  Fig.  8.  The  ratemeters  used  were  built 
similar  to  the  Baird- Atomic,  Inc.,  Model  ^32  A,  Research  Ratemeter.  NE  901, 
the  left-hand  ratemeter  in  Fig.  8,  was  used  for  the  lower  detector  level,  and 
NE  902  was  used  for  the  upper  detector  level.  The  time  constants  were  deter- 
mined experimentally  by  a  step  input  to  each  ratemeter. 

When  the  d-c  output  of  the  integrator  circuit  in  a  ratemeter  exceeded 
a  given  value,  firing  of  a  neon  bulb  caused  the  current  to  flow  in  the  final 
stage  of  the  ratemeter.   Two  outputs  were  provided  from  each  ratemeter.   One 
output  from  each  ratemeter  was  sent  to  a  two- channel  high-speed  recorder,  and 
the  other  output  from  each  ratemeter  served  as  the  input  to  a  Time  Interval 
Meter. 

Two- channel  recorder.   The  counting  rate  at  each  detector  level  was  moni- 
tored by  one  channel  of  a  dual  channel  recorder.   The  recorder,  pictured  on  the 
stool  in  Fig.  8,  was  an  Offner  Type  RS  Dynograph,  NE  733.   The  output  from  the 
lower  detector  level  ratemeter  was  the  input  to  the  left-hand  channel  of  the 
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Figure     8.      The     Experimental      Apparatus. 
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recorder.   The  output  from  the  upper  detector  level  was  the  input  to  the  right- 
hand  channel  of  the  dual- channel  recorder. 

Time  Interval  Meter.   The  Time  Interval  Meter,  shown  below  the  ratemeters 
in  Fig.  8,  was  a  Be ckman/ Berkeley  Model  7250R,  NE  763.  The  counting  rates 
from  the  ratemeters  were  fed  to  the  input  channel  amplifiers  of  the  Time  Inter- 
val Meter  for  pulse  shaping  and  amplification.  Whenever  a  sharp  change  in  the 
counting  rate  occurred,  a  positive  pulse  was  generated  at  the  output  of  the 
Time  Interval  Meter.   The  positive  pulse  output  of  each  input  channel  of  the 
Time  Interval  Meter  was  fed  to  a  single,  specially  designed,  transistorized, 
monostable  multivibrator  circuit,  KSU-NE  No.  009«   The  multivibrator  circuit 
is  shown  in  Fig.  ^.      The  output  from  this  circuit  drove  the  external  delta-t 
circuit  of  the  Technical  Measurements  Corporation  Multi-Channel  Analyzer 
Multi-Sealer  Logic  Unit.  All  other  outputs  from  the  input  channel  circuits 
of  the  Time  Interval  Meter  were  disconnected.   The  gate  times  were  stopped  by 
the  rapid  change  in  the  counting  rate.   To  restart,  the  oscillator  circuit  in 
the  Time  Interval  Meter  was  utilized.   To  provide  the  restart  pulse,  the  fol- 
lowing modifications  were  made  to  the  Time  Interval  Meter.   The  output  of  pin 
1  of  V^K^  in  the  oscillator  and  shaper  circuit  in  the  Time  Interval  Meter  was 
disconnected  and  reconnected  to  the  top  of  R206  through  a  ^70  mmfd  capacitor. 
The  circuit  previously  feeding  R206  was  disconnected  by  removing  C201.   These 
modifications  can  be  seen  by  observing  the  schematic  diagrams  in  the  instruction 
manual  for  the  Time  Interval  Meter.   The  Oscillator  and  Shaper  Schematic  Dia- 
gram and  the  Master  Schematic  Diagram  are  Beckman/Berkeley  Division  Drawing 
Numbers  7000  C  U19  and  7250  D  21,  respectively. 

The  oscillator  in  the  Time  Interval  Meter  also  served  as  a  fixed  frequency 
pulse  which  acted  as  a  time  base  generator  and  was  fed  into  the  amplifier  input 
of  the  Multi-Channel  Analyzer  Multi-Sealer  Logic  Unit.  A  four-decade  circuit 
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Figure      9.       Wiring      Diagram      for      KSU-NE      009      Multivibrator      Circuit. 
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with  an  output  jack  for  each  decade  was  added  to  the  output  of  the  oscillator 
and  shaper  circuit  so  that  a  variable  selection  of  frequencies  could  be  used. 
Only  one  pulse  frequency  was  used  to  obtain  the  experimental  data. 

Mult i- Channel  Analyzer.   The  Technical  Measurements  Corporation  256- 
channel  Model  CN-110  Analyzer,  NE  630,  utilizing  the  Multi-Sealer  Logic  Unit, 
served  as  a  counter  of  the  number  of  pulses  from  the  Time  Interval  Meter  oscilla- 
tor over  a  period  of  time.   The  analyzer  is  shown  in  Fig.  10.   The  Pulse  Height 
Logic  Unit,  NE  633>  as  shown  in  Fig.  10,  was  not  used.   Instead,  as  noted  above, 
the  Multi-Sealer  Logic  Unit,  Model  214,  NE  633,  was  used. 

The  output  pulse  from  the  mono- stable  multivibrator  circuit  was  connected 
to  the  external  delta-t  circuit  of  the  Multi-Sealer  Logic  Unit  and  served 
as  a  trigger  to  switch  channels  in  the  multi- channel  analyzer.  Each  channel,  one- 
by-one,  would  accumulate  and  store  the  number  of  oscillator  pulses  counted  be- 
fore another  rapid  rise  in  counting  rate  from  one  of  the  ratemeters  would  cause 
storage  to  be  transferred  to  the  next  channel.   Therefore,  the  number  of  pulses 
counted  from  the  Time  Interval  Meter  oscillator  was  a  measure  of  the  time  be- 
tween rapid  change  in  the  output  of  the  ratemeters,  resulting  from  the  radio- 
active source  passing  one  of  the  detector  level  slits. 

Digital  Print  Out.  A  Hewlet  Packard,  Model  561  B,  NE  632,  recorded  the 
information  stored  in  each  of  the  256  channels  of  the  Multi- Channel  Analyzer. 
The  printer  is  shown  directly  above  the  Multi-Channel  Analyzer  in  Fig.  10 

2.2  Materials 

Type  M40-C  stainless  steel  ball  bearings  of  reject  quality  obtained  from 
the  Hartford  Steel  Ball  Company  were  used  exclusively  for  the  solids  in  the 
fluidized  bed.   The  balls  had  a  sample  mean  diameter  of  O.I876  inch,  with  a 
variance  of  0.0003  (inch)2,  based  on  a  random  sample  of  30  stainless  steel 
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Figure     10.    Multi  -  Channel     Analyzer    and     Digital 
Printout     Units. 
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balls  removed  from  the  fluidized  bed.   The  mean  sample  weight  of  the  balls  was 
0.430^  gram. 

The  radiation  source  was  a  type  316  stainless  steel  ball  bearing,  irradiated 
at  a  flux  of  10^-2  neutrons  cnT^sec   in  the  lazy  susan  irradiation  facility  of 
the  Kansas  State  University  Triga  Mark  II  nuclear  reactor.   Irradiation  for  l|r 
hours  produced  a  Manganese- 56  source  strength  of  12  to  lh  millicuries. 

A  Type  316  stainless  steel  ball  bearing  was  used  because  of  the  2$  man- 
ganese content  by  weight  in  Type  316  stainless  steel.   The  diameter  of  the  iso- 
tope was  0.1875  inch,  and  its  weight  was  0.450  gram.   Industrial  Tectonics, 
Inc.  supplied  the  Type  316  stainless  steel  ball  bearings.  A  different  ball 
bearing  was  used  as  the  radioactive  source  each  day  data  were  obtained.   No 
detectable  variations  were  found  in  the  size  and  weight  of  the  individual  ball 
bearings. 

Figure  3  is  a  photograph  showing  the  stainless  steel  ball  bearings  inside 
the  fluidizing  column.   The  total  weight  of  ball  bearings  in  the  column  was 
2^.125  pounds,  and  the  static  bed  height  was  19^5  inches. 

2.3  Calibrations 

2.3-1  Half -life  Determination 

The  half-life  of  the  radioisotope  produced  upon  neutron  irradiation  in  the 
Kansas  State  University  Triga  Mark  II  nuclear  reactor  was  determined  by  measuring 
the  activity  one  yard  from  the  source  with  a  dose-rate  meter,  at  various  times 
after  removal  of  the  stainless  steel  ball  bearing  from  the  nuclear  reactor. 
Figure  11  is  a  plot  of  the  log  of  the  dose  rate  at  one  yard  from  the  source 
versus  the  time  after  removal  of  the  ball  bearing  from  the  reactor.  The  half- 
life  determined  from  this  figure  was  found  to  be  2.5  hours.  Manganese-56  has 
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a  gamma- ray  half-life  of  2.58  hours;  therefore,  the  activity  produced  was 
primarily  Manganese- 56. 

2.3.2  Calibration  of  the  Orifice  Meter 

A  length  of  1§- inch- inside  diameter  hose  was  connected  to  the  outlet  of 
the  column,  and  the  other  end  of  the  hose  was  placed  in  a  ^-gallon   barrel. 
Water  was  pumped  from  the  supply  tank,  through  the  fluidizing  column,  and  into 
the  55-Sall°n  barrel.   Two  levels  were  marked  off  in  the  barrel.   The  water 
flow  was  adjusted.   The  height  of  the  mercury  in  the  manometer  across  the  ori- 
fice meter  and  the  time  taken  for  the  water  level  in  the  barrel  to  pass  from 
the  first  mark  to  the  second  mark  were  recorded.  Figure  12  shows  a  plot  of 
gallons  per  minute  versus  the  manometer  reading.  Figure  13  converts  the 
manometer  reading  to  superficial  fluid  velocity  in  the  column. 

2.3.3  Radial  Effects 

The  scintillation  detectors  were  paired  off  according  to  their  operating 
voltage.  It  was  found  that  of  the  scintillation  probes  available,  NE  287  and 
NE  31k,   and  HE  288  and  NE  353  formed  the  best  pairs. 

The  radioactive  source,  in  its  vertical  movements,  would  rarely  be  moving 
in  the  center  of  the  column  during  a  data  run.   It  was,  therefore,  necessary 
to  determine  the  effect  of  radial  position  on  the  counting  rates  from  each 
pair  of  detectors  when  the  particle  passed  through  the  plane  of  each  detector 
level . 

The  source  was  placed  in  each  of  the  radial  positions  in  the  fluidizing 
column  as  shown  in  Fig.  Ik.     For  each  radial  position,  counting  rates  were 
obtained  at  different  vertical  positions  with  respect  to  the  detector  level. 
Counting  rate  versus  vertical  position  curves  were  plotted  for  each  radial 
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position  and  each  pair  of  detectors.   These  curves  are  shown  in  Figs.  15,  l6, 
17 _,  and  18.   These  figures  show  that,  with  the  arrangement  of  detectors  used, 
the  radial  position  of  the  radioactive  source  had  little  effect  on  the  out- 
put received  from  the  scintillation  probes.   Detectors  288  and  353  gave  a 
lower  counting  rate  than  detectors  287  and  31^-  at  the  same  high  voltage  setting; 
therefore,  during  the  data  runs  the  high  voltage  to  detectors  288  and  353  vas 
adjusted  to  a  higher  value  to  equalize  the  sensitivity. 

2.3'^  Gap  Width  and  Thickness  of  Shield  Determination 

It  was  necessary  to  determine  the  optimum  gap  width  and  shielding  thick- 
ness to  provide  a  high  counting  rate  when  the  source  was  in  front  of  the  gap 
and  a  low  counting  rate  when  the  source  was  not  in  front  of  the  gap.   The  gap 
width  and  shielding  thickness  were  determined  experimentally  by  measuring  the 
counting  rate  at  different  vertical  distances  with  respect  to  the  detector 
level  for  different  gap  widths  and  shielding  thicknesses. 

The  source  was  placed  in  a  plastic  container  and  suspended  into  the  water- 
filled  column.  It  was  moved  vertically  with  a  nylon  cord.  A  graduated  scale 
was  placed  alongside  the  column  and  the  radioactive  source  was  held  at  constant 
vertical  positions  while  one-minute  counts  were  obtained.  Figures  19,  20,  21, 
22,  and  23  show  the  counting  rates  versus  vertical  position  for  different  gap 
widths  and  shield  thicknesses.   The  gap  width  and  shield  thickness  selected 
are  those  shown  in  Fig.  23 • 

2.3*5  Calibration  of  the  Electronic  System 

As  the  source  passed  from  one  detector  level  to  the  other,  the  outputs 
from  the  ratemeters  were  recorded  as  a  pulse  on  the  dual- channel  recorder. 
The  right-hand  channel  of  the  recorder  would  map  the  output  from  the  upper 
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detector  level,  and  the  left-hand  channel  vould  map  the  output  from  the  lower 
detector  level.  Knowing  the  speed  of  the  recorder,  the  distance  between  pulses 
was  used  to  determine  the  elapsed  time  between  the  pulses.   If  the  right-hand 
channel  recorded  a  pulse,  and  the  next  pulse  recorded  was  by  the  left-hand 
channel,  the  particle  had  passed  from  the  upper  to  the  lower  detector  level. 
By  measuring  the  distance  between  the  two  pulses,  the  time  to  travel  from  upper 
to  lower  level  was  obtained. 

The  known  fixed- frequency  input  to  the  Multi-Sealer  Logic  Unit  of  the 
Technical  Measurements  Corporation  Multi-Channel  Analyzer  acted  as  a  timer. 
When  a  pulse  was  recorded  on  either  of  the  two  recorder  channels,  the  same 
pulse  provided  a  trigger  to  switch  channels  in  the  multi-sealer.   The  number 
of  counts  stored  in  each  channel  of  the  multi- scaler  was  converted  to  elapsed 
time  by  dividing  by  the  oscillator  frequency. 

The  reliability  of  the  measurement  of  elapsed  time  was  determined  experi- 
mentally.  This  was  done  by  moving  the  source  past  the  two  detector  levels  at 
precisely  known  time  increments  and  comparing  this  to  the  time  increments  ob- 
tained from  the  two  methods  outlined  above.   The  experimental  time  was  deter- 
mined by  finding  the  pulses  on  the  chart  recording  and  correlating  these  pulses 
to  the  correct  channel  locations  given  by  the  printout  of  the  multi- scaler 
memory  unit.   Once  the  channel  number  of  the  printout  tape  was  correlated  to  a 
pulse  on  the  chart  recording,  the  number  of  counts  stored  in  that  channel  was 
divided  by  the  frequency  of  the  oscillator  to  obtain  the  time  increment.   This 
gave  the  experimentally  measured  times  for  the  particle  to  pass  from  the  lower- 
to-upper  and  upper-to-lower  detector  levels. 

A  Zero-Max  motor  and  gearhead  unit  was  used  to  move  the  source  up  and  down 
in  the  column.  The  unit  was  capable  of  variable  speeds  from  0  to  400  revolutions 
per  minute.  A  pulley  was  mounted  at  the  top  of  the  column,  and  a  nylon  cord  was 
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passed  over  the  pulley  into  the  column.   The  source  was  suspended  on  the  end 
of  the  cord.   The  other  end  of  the  cord  was  attached  to  the  shaft  of  the  gear- 
head.  By  winding  the  cord  around  the  shaft  until  the  source  was  at  the  top 
of  the  column  and  turning  on  the  motor,  the  source  would  be  lowered  to  the 
bottom  of  the  column  and  returned  to  the  top  again.   By  this  technique  the 
source  was  lowered  and  raised  at  a  constant  speed  past  the  two  detector  levels. 
A  30-inch  section  of  the  column  was  marked  off,  and  the  time  for  the  source 
to  be  lowered  through  this  section  was  taken  as  the  best  estimate  of  the  time 
increment  and  was  taken  to  be  the  true  time. 

Experimental  runs  were  made  at  several  gearhead  speeds,  and  the  time  incre- 
ments required  to  pass  between  the  two  levels  obtained  by  the  combination  of 
the  two  methods  discussed  previously  were  compared  to  the  visual  observations. 

Several  sources  of  error  were  found.  The  multi-sealer  occasionally  skip- 
ped channels,  thereby  reducing  the  number  of  data  values  obtained  from  each 
data  run.  Due  to  mechanical  malfunctions  of  the  rotary  printing  head,  the 
digital  printer  often  printed  incorrect  numbers.  For  time  increments  of  less 
than  about  one  second  required  for  the  radioactive  source  to  pass  between  the 
two  detector  levels,  the  time  increments  obtained  from  the  experimental  appar- 
atus deviated  largely  from  the  actual  time  increments.   To  minimize  the  latter 
error,  the  following  technique  was  used. 

The  times  for  the  source  to  move  from  level  Z.  to  level  Z  and  from  level 

1  2 

Z  to  level  Zn  for  a  given  source  speed  measured  experimentally  were  averaged 
and  plotted  versus  the  true  time  for  different  true  times.   The  results  are 
shown  in  Fig.  2k.     At  short  time  intervals  the  deviation  from  a  one-to-one 
relationship  is  seen  clearly  in  this  figure.  Figure  25  gives  the  same  infor- 
mation, plotted  as  the  ratio  of  the  average  experimental  time  to  the  true  time 
versus  the  average  experimental  time.  A  significant  deviation  for  time 
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increments  of  less  than  one  second  is  shown  to  be  the  case. 

Figure  26  gives  the  ratio  of  experimental  time  to  true  time  as  a  function 
of  the  experimental  time  before  averaging  the  time  increments  obtained  experi- 
mentally for  the  radioactive  source  to  pass  upward  and  downward  between  the 
two  detector  levels  at  a  given  source  speed.  The  triangles  represent  the 
experimental  time  increments  obtained  for  a  given  source  speed  when  the  source 
was  moving  downward  between  Z  and  Z  ,  and  the  circles  represent  the  experi- 
mental time  increments  obtained  for  a  given  source  speed  when  the  source  was 
moving  upward  between  Z  and  Z  . 

The  time  increment  data  obtained  when  the  source  was  moving  downward  were 
fitted  with  three  straight-line  segments.  A  linear  regression  analysis  was 
utilized  to  obtain  the  equation  for  each  straight-line  segment.  A  similar 
procedure  was  used  to  fit  straight-line  segments  to  the  time  increment  data 
obtained  when  the  source  was  moving  upward.   The  equations  resulting  from  the 
linear  regression  analysis  of  each  straight-line  segment  are  given  in  Table  1. 

An  IBM  1620  computer  program  in  FORTRAN  was  developed  to  treat  all  of  the 
lower- to-upper  level  and  upper- to- lower  level  time  increment  data  obtained  from 
the  experimental  apparatus.  For  a  given  experimental  time,  the  ratio  of  ex- 
perimental time  to  true  time  was  obtained  from  the  appropriate  equation  shown 
in  Table  1.   To  obtain  the  true  time,  the  given  experimental  time  increment 
was  divided  by  that  ratio. 

The  computer  program  developed  to  correct  the  experimental  data  is  located 
in  Appendix  A. 

2.3*6  Time  Constants  of  the  Ratemeters 

The  time  constants  for  the  ratemeters  were  determined  by  feeding  a  step 
input  to  each  ratemeter  and  using  a  strip  chart  recorder  to  map  the  output  of 
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Table  1.  Equations  resulting  from  the  linear 
regression  analysis  of  the  data  from 
the  calibration  of  the  electronic 
system. 


Detector  level  1  to  detector  level  2 


Experimental  time                               Equation 
( seconds) 

0.00  to  0.12  ln(R)  =  1.006273  +  0.7186  ln(t) 

0.12  to  1.50  ln(R)  =  -0.21715  +  0.13258  ln(t) 

1.50  to  infinity  R  =0.85 

Detector  level  2  to  detector  level  1 

0.00  to  0.14  ln(R)  =  1.33479  -*•  0.65478  ln(t) 

0.14  to  infinity  R  =  1.1734 

R  =  Experimental  time/True  time 
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each  ratemeter.  The  procedure  vas  as  follows:   a  square-wave  input  was  fed  to 
each  ratemeter,  the  sensitivities  "being  so  adjusted  at  first  that  the  input 
would  not  "be  accepted;  the  sensitivity  of  each  ratemeter  was  then  increased 
to  accept  the  square-wave  input;  the  resulting  output  of  each  ratemeter,  map- 
ped by  a  strip  chart  recorder,  was  allowed  to  reach  steady  state;  the  sensi- 
tivity was  then  decreased  so  as  not  to  accept  the  input,  and  the  output  was 
allowed  to  return  to  the  original  no-input  value. 

The  time  required  to  reach  (l-l/e)  of  the  steady  state  value  is  equal 
to  the  RC  time  constant  for  the  rise  in  the  output,  and  the  time  required  to 
reach  (l/e)  of  the  steady  state  value  is  equal  to  the  RC  time  constant  for 
the  decay  of  the  output.  Figures  27  and  28  illustrate  the  method  of  obtaining 
these  time  constants. 

The  time  constants  obtained  for  the  two  ratemeters  in  the  experimental 
apparatus  are  given  in  Table  2. 

Table  2.   Time  constants  of  the  ratemeters. 

Ratemeter RC  rise  time RC  decay  time 

NE  901  0.492  second  0-300  second 

RE  902  0.500  second  O.336  second 

2.3.7  Pressure  Drop  Data 

The  pressure  drop  across  the  fluidized  bed  was  determined  by  recording 
the  pressure  difference  in  inches  of  mercury  across  the  column  for  different 
flow  rates  of  water.   Readings  were  taken  arbitrarily  as  the  flow  rate  was 
either  increased  or  decreased.  Figure  29  gives  the  pressure  drop  in  inches  of 
mercury  as  a  function  of  the  water  flow  rate  in  gallons  per  minute. 
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Figure    27.    Method    of   Determining    Rise    Time    of 
Rotameters. 
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Figure    28.    Method    of    Determining     Decay    Time    of 
Ratemeters. 
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2.h     Procedure 

To  determine  if  the  equipment  was  operating  properly,  the  nuclear  instru- 
mentation equipment  was  checked  out  prior  to  a  day's  run.   The  check-out  pro- 
cedure consisted  of  allowing  the  electronic  apparatus  to  warm  up  and  counting 
the  pulses  from  a  fixed- frequency  input  to  the  multi- scaler  for  one  second 
in  each  channel  of  the  multi-sealer.   The  digital  printer  was  then  used  to 
print  out  the  information  in  the  memory  of  the  multi- scaler.   If  the  number 
stored  in  each  channel  was  equal  to  the  frequency  input  in  cycles  per  second, 
the  equipment  was  assumed  to  he  operating  properly.   The  high  voltage  to  the 
photomultipliers  of  the  scintillation  probes  was  then  adjusted,  and  a  radio- 
active source  was  passed  in  front  of  the  scintillation  detectors  to  check  pro- 
per operation  of  the  dual- channel  recorder  used  to  measure  the  output  of  the 
two  ratemeters.   Next,  another  check  of  the  multi- scaler  was  made  by  observing 
the  cathode  ray  tube  to  determine  if  the  source  passing  in  front  of  the  scin- 
tillation probes  was  triggering  the  unit  properly.   If  the  system  was  operating 
properly,  a  Type  3l6  stainless  steel  ball  bearing  was  introduced  into  the 
Kansas  State  University  Triga  Mark  II  nuclear  reactor  and  was  irradiated  at 
100  kw  for  one  and  one-half  hours.  A  surface  dose  rate  of  about  2  r/hr  of 
gamma  radiation  was  produced.   The  source  was  then  inserted  into  the  fluid- 
izing  column,  and  the  pump  was  turned  on  to  allow  the  water  to  circulate  through 
the  system. 

The  high  voltage  to  the  photomultipliers  at  each  of  the  detector  levels 
was  adjusted  until  well-defined  pulses  with  minimum  background  noise  were  ob- 
tained on  each  channel  of  the  chart  recorder. 

The  distance  between  the  detector  levels  was  set  for  the  particular  data 
run,  and  the  water  flow  rate  for  the  run  was  properly  adjusted  by  setting  to 
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the  desired  value  the  mercury  height  in  the  manometer  across  the  orifice  in 
the  inlet  line  to  the  fluidizing  column. 

The  chart  recorder  was  then  started.  The  multi- scaler  was  started,  and 
the  exact  position  of  the  pens  on  the  chart  recorder  at  the  time  of  starting 
was  noted. 

During  the  data  run,  the  water  flow  rate  was  maintained  constant  by  making 
the  necessary  inlet  valve  adjustments  to  maintain  the  proper  height  of  mercury 
in  the  manometer. 

Observation  of  the  cathode  ray  tube  of  the  multi- channel  analyzer  pro- 
vided a  running  check  on  the  approximate  amount  of  time  left  on  a  particular 
run.  An  illuminated  dot  on  the  face  of  the  tube  monitored  the  channel  currently 
being  used  to  store  the  frequency  count.  As  soon  as  the  source  passed  one  of 
the  detector  levels,  storage  was  transferred  to  the  next  channel.  This  could 
be  visually  observed  by  watching  the  dot  jump  to  the  right  on  the  screen.  When 
the  256  channels  were  filled,  the  dot  would  then  move  to  an  arbitrary  position 
to  the  left  on  the  screen.   This  indicated  the  end  of  the  data  run,  and  the 
position  of  the  pens  on  the  chart  recorder  was  again  noted. 

The  manometer  reading,  the  distance  between  detector  levels,  the  run  num- 
ber, and  the  data  were  then  recorded  on  both  the  chart  recording  and  the  memory 
information  of  the  multi- scaler  printed  out  by  the  digital  printer. 

The  detectors  were  positioned  at  different  levels.   Because  of  the  strong 
gamma-ray  source  in  the  fluidized  bed,  changes  were  made  as  quickly  as  possible. 
The  flow  rate  of  the  water  was  adjusted  to  a  new  value,  and  a  new  run  was  be- 
gun. 

The  length  of  time  available  for  obtaining  the  data  for  each  day ' s  run 
was  determined  by  the  half -life  of  the  radioactive  source.  As  time  passed, 
in  order  to  maintain  well-defined  pulses  on  the  dual-channel  chart  recorder, 
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the  high  voltage  to  the  photomultipliers  was  increased  to  increase  detector 
sensitivity,  which  compensated  for  the  decreased  activity  of  the  source.   When 
the  high  voltage  was  increased  too  high,  the  background  noise  would  increase, 
resulting  in  difficulty  in  distinguishing  between  pulses  that  actually  occurred 
when  the  source  passed  a  detector  level  and  background  noise.   This  limited 
the  time  for  taking  data  to  about  seven  or  eight  hours  after  removal  of  the  ball 
bearing  from  the  reactor.  The  source  would  have  decayed  through  about  three 
half-lives  by  that  time  and  was  roughly  1/8  of  its  activity  at  the  time  of  re- 
moval from  the  nuclear  reactor. 
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3.0  ANALYSIS  OF  SINGLE  PARTICLE  EXPERIMENTS 
3.1  Experiment  Arrangement 

Figure  30  is  an  illustration  of  the  fluidized  column  with  specific  ver- 
tical positions,  Z. ,  noted.   Utilizing  the  experimental  apparatus,  the  time 
a  radioactive  particle  passed  level  Z.  was  recorded.   Only  two  levels  of  Z. 
could  he  investigated  at  a  time.   The  distance  between  the  two  Z^  levels 
was  the  distance  between  the  two  detector  levels,  A  Z. 

The  controlled  variables  were  the  fluidized  bed  height,  determined  by 
the  water  flow  rate,  the  quantity  of  particles  in  the  bed,  and  the  distance 
between  the  two  detector  levels. 

3.2  Quantities  Extracted  from  Data 

The  quantities  obtainable  from  a  single  particle  experiment  are  given 
below: 

t(Z^,  Z.)     Residence  time,  defined  as  the  time  required  for  a 

particle,  after  passing  level  Z.  and  moving  into  the 
region  between  Z.  and  Z.  and  not  returning  to  Z.,  to 

to  pass  out  of  the  region  at  Z.. 

J 

r(Z. )        Particle  traffic,  defined  as  the  number  of  times  a 

* 

particle  passes  level  Zi  per  unit  time. 
Pr(Zi,ZjjZi)  Reflectivity,  defined  as  the  probability  of  a  particle 

passing  level  Zi  and  not  reaching  level  Z.  before  re- 
el 

turning  to  level  Z^. 

l(Z. ,Z.,t)    Particle  concentration,  defined  as  the  number  of 

times  a  particle  passes  level  Z.  in  the  time  increment 

J 
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At  at  the  time  t  per  unit  time,  initially  passing 


level  Zi# 
m(Z.)        The  number  of  times  a  particle  passes  level  Z.  be- 
fore reaching  level  Z  . 

3.3  Determination  of  Fundamental  Quantities 
from  Experimental  Data 

Figure  31  is  a  simplified  illustration  of  a  data  trace,  showing  the  times 
at  which  a  particle  passed  levels  Z^  and  Z..   The  fundamental  quantities  are 
determined  from  Fig.  31  as  follows: 

t(z±,  z.)  -   0.175;  0.175;  0.20 
t(Zj,  Zi)  =  0.19;  0.125;  0.125;  0.18 

r(Zi)  =  16/4  -    k.O 
r(Zj)  =  iQ/k   =  4.5 

Obtaining  the  reflectivity  and  particle  concentration  data  requires  some 
explanation.   To  obtain  the  reflectivity  data  Pr,  m(Z)  must  be  considered. 
For  m  particles  passing  level  Z^,   in  general,  one  half  of  them  are  headed 
toward  Z.,  and  one  eventually  reaches  Z..   Therefore,  the  probability  of  that 
particle  reaching  Z.  in  m/2  tries  is  l/(m/2)  or  2/m.   The  probability  of  not 
reaching  Z*   but  returning  to  Z-j_  is,  therefore,  l-(2/m)  or  (m-2)/m. 

From  Fig.  31>  the  probability  of  reflection  of  particles  starting  at 
level  Z  is,  in  the  first  instance,  (K-2)/k,   since  the  radioactive  source 
passed  the  detector  level  Z^   four  times  before  reaching  level  Z,.   Tne  next 
time  the  source  was  detected  at  level  Z. ,  having  been  detected  at  level  Z. 
just  prior  to  detection  at  level  Zi,  it  passed  level  Z.  six  times  before  it 
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Figure   30.    Vertical     Arrangement    of    Detectors 
Adjacent     to     Fluidizing     Column. 
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Figure    31.    Simplified     Data     Trace. 
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reached  Z*,   giving  (6-2)/6.   The  third  time  the  source  was  detected  at  level 
Z.,  having  "been  detected  at  level  Z.  just  prior  to  detection  at  level  Z.,  it 
passed  level  Z^  twice  before  "being  detected  at  level  Z..   The  same  analysis 
is  used  to  determine  the  probability  of  reflection  of  particles  starting 
at  level  Z.. 

PrCz^ZjjZ^  =  (4-2) A  =  0.5;  (6-2)/6  -  0.667;  (2-2)/2  -  0. 

Pr(Zj,Zi;Zj)  =  (6-2)/6  =  0.667;  (2-2)/2  =  0;  (8-2)/8  =  0.75; 
(2-2)/2  =  0. 

To  obtain  the  particle  concentration  data,  the  number  of  times  a  particle 

passes  a  level  per  unit  time  within  a  time  increment  At   is  considered.  Also 

considered  is  the  level,  Zi  or  Z,,  at  which  the  particle  originated.  As  an 

example,  l(Z.,Z.,t)  and  l(Z^,Z.,t)  will  be  obtained  step-by- step.  For  con- 

venience,  4t  will  be  set  equal  to  1.   Since  the  second  pulse  at  level  Z. 

shown  in  Fig.  31  occurs  at  the  beginning  of  the  second  time  increment,  it  will 

be  used  for  the  initial  pulse.   Considering  l(Z.,Z.,t),  the  number  of  pulses 

J 

occurring  at  level  Z.  within  each  time  increment  At   with  the  time  scale  begin- 
ning at  the  initial  pulse  is 

l(Z.,ZJ,ti)=2,  l(Zi,ZJ,t2)=7,  and  l(Zi,Z(j,t3)=3. 

Considering  the  number  of  pulses  occurring  at  level  Z.  in  the  same  manner 
beginning  with  the  initial  pulse, 

l(Zi,Zi,t  )=7,  l(Zi,Zi,t2)=2,  and  l(Zi,Z.,t3)  =  6. 
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k.O     DATA  ANALYSIS 
k.l     General  Consideration 

Figure  32  shows  an  example  of  a  typical  recorder  trace,  with  arrows 
a*  ded  to  illustrate  the  upward  or  downward  movement  of  the  particle  as  it 
passed  the  two  detector  levels.  Figure  32  also  illustrates  a  typical  digital 
printer  recording  of  the  first  fifty  channels  of  the  multi-channel  analyzer. 
The  left-hand  column  of  numbers  indicates  the  channel  number,  and  the  right- 
hand  column  lists  the  number  of  pulses  from  the  Time  Interval  Meter  oscillator 
stored  in  each  channel. 

For  each  given  flow  rate  and  distance  between  detector  levels,  the 
recorder  trace  and  the  corresponding  digital  printer  recording  were  correlated. 
This  was  done  by  using  the  recorder  trace  to  determine  when  the  particle 
passed  either  the  bottom  or  the  top  detector  level  and  by  matching  the  time 
between  pulses  on  the  recorder  trace  to  the  numbers  recorded  by  the  digital 
printer.   The  dual  channel  recorder  was  always  operated  at  5  millimeters  per 
second,  so  the  graduations  on  the  recorder  trace  were  one  second  apart.   The 
pulses  from  the  Time  Interval  Meter  to  the  Multi scaler  Logic  Unit  amplitude 
input  were  always  adjusted  to  666   per  second;  therefore,  one  graduation  on 
the  recorder  trace  corresponded  to  the  number  666   on  the  digital  printer 
recording. 

The  channel  number  on  the  digital  print- out  recording  was  always  equated 
to  the  pulse  on  the  recorder  trace  that  stopped  the  accumulation  of  Time 
Interval  Meter  pulses  in  that  channel  and  transferred  the  storage  to  the  next 
channel.  Each  pulse  on  the  recorder  trace  that  caused  transfer  of  storage 
from  one  channel  to  the  next,  therefore,  was  identified  by  numbering  it  with 
the  proper  channel  number. 
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Figure     32.     Output  of    Recorder    and     Digital      Printer. 
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The  left-hand,  or  "bottom  detector  level  trace,  was  denoted  by  a  minus 
sign,  so  that  the  number  on  the  digital  print- out  recording  corresponding  to 
a  pulse  that  occurred  when  the  source  passed  the  bottom  detector  level  was 
identified  by  placing  a  minus  sign  beside  it.   The  number  on  the  digital  print- 
out recording  corresponding  to  a  pulse  that  occurred  when  the  source  passed 
the  top  detector  level  was  identified  by  placing  a  plus  sign  beside  it.   The 
digital  print- out  recording  with  a  plus  or  minus  sign  beside  the  number  stored 
in  each  channel  in  the  multi-channel  analyzer,  therefore,  contained  all  of 
the  information  obtained  from  the  recorder  trace  and  the  digital  printer  record- 
ing.  By  observing  the  plus  and  minus  signs,  it  was  easy  to  determine  when  the 
source  passed  from  one  detector  level  to  the  other.  For  instance,  when  the 
isotope  passed  from  lower  to  upper  detector  level,  the  sign  would  change  from 
minus  to  plus.   The  time  required  to  pass  from  lower  to  upper  detector  level 
was  found  by  dividing  by  666  the  number  of  pulses  counted  and  stored  in  the 
channel  denoted  by  the  plus  sign.   The  resulting  time  would  be  in  seconds. 

The  number  in  each  channel  location  of  the  digital  printer  recording, 
along  with  its  corresponding  plus  or  minus  sign,  was  then  transferred  to  an 
IBM  data  card.   To  reduce  the  amount  of  data  cards  used,  eight  of  the  positive 
or  negative  numbers  were  typed  per  card. 

Not  all  of  the  sets  of  data  obtained  from  the  experimental  apparatus 
contained  256  data  values,  principally  because  the  multi- channel  analyzer 
would  sometimes  skip  channels,  causing  no  numbers  to  be  recorded  in  the  skipped 
channels.   Therefore,  each  set  of  data  was  accompanied  by  the  number  of  data 
values  in  the  set.   This  value,  denoted  by  JMAX,  preceded  the  set  of  data  on 
an  IBM  data  card  by  itself. 

Utilizing  five  water  flow  rates  through  the  column  and  five  distances  be- 
tween detector  levels,  seventy- six  sets  of  data  were  obtained  from  the  experiments. 


55 


Table  3  is  a  listing  of  the  number  of  data  runs  obtained  for  each  given 
water  flow  rate  and  distance  between  detector  levels.   B^  is  the  water  flow 
rate  corresponding  to  the  height  of  mercury  in  the  manometer  that  measured 
the  pressure  drop  through  the  orifice  meter,  and  AZ  is  the  distance  between 
the  two  detector  levels. 

Table  h   is  a  listing  of  the  data  set  numbers  for  given  distances  between 
detector  levels  and  water  flow  rates. 

k.2     Statistical  Analysis  of  Data 

Sample  variation  was  found  to  be  useful  in  describing  the  frequency  distri- 
bution of  residence  times  obtained  from  the  experimental  data.   Before  any 
test  of  variance  can  be  used  on  a  sample  taken  from  a  population,  however,  the 
distribution  upon  which  the  test  of  variance  is  based  must  first  be  determined. 

The  simplest  distribution  for  continuous  data  is  the  normal  distribution, 
where  the  mean  and  median  are  the  same  and  the  distribution  is  symmetric  with 
respect  to  the  mean.   The  arithmetic  mean  is  defined  as  (ZLffX./^^),  and  the 

median  is  defined  as  the  value  equaled  or  exceeded  by  exactly  half  of  the  values 

2 
x.  in  the  given  list.  An  unbiased  estimator  of  the  population  variance  CT  '   is 


(i) 

n  -  l 

where  x.  is  the  variable  of  which  n  observations  are  made,  and  x  is  the 
average  or  mean  of  the  n  observations. 

A  measure  of  the  skewness  is  given  by 
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Table  3.  Number  of  data  sets  obtained  for 
each  B^,  A  Z. 


AZ  (inches)  £  1 

Bx  (2.50  in.  Hg)  2  3 

B2  (4.00  in.  Hg)  2  4 

B3  (4.85  in.  Kg)  2  4 

B4  (7.00  in.  Hg)  2  4 


B5  (8.50  in.  Hg) 


Ta"ble  4.      Data   set  nunoers   for  given  distances 
"between   detector  levels   A  Z  and  water 
flow  rates  B^# 
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AZ 
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1 

2 

4 

6 

Bl 

1   2 

3  4   5 

6   7   8   9 

10  11  12  13 

14  15 

^2 

16  17 

18  19  20  21 

22  23  24  25 

26  27  28 

29  30 

B3 

31  32 

33  34  35  36 

37  38  39  40 

41  42  43  44 

45  46 

B4 

47  48 

49  50  51  52 

53  54  55  56 

57  58  59  60 

61  62 

*5 

63  64 

65  ,66  67  68 

69  70.  71  72 

73  74 

75  76 

The  average  fluid  velocities  for  the  given  water  flow 
rates  through  the  fluidized  "bed  are  given  "below. 


B- 


Br 


B, 


B, 


Bc 


10,4  in. /sec. 
12.1  in. /sec. 
12.8  in. /sec. 
14.1  in. /sec. 
14.8   in. /sec. 


r 
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A  positive  value  of  S-5  indicates  that  the  mean  is  greater  than  the  median, 
and  a  negative  value  of  S-*  indicates  that  the  median  is  greater  than  the 
mean. 

Figure  33  shows  an  example  of  a  skewed  distribution  where  the  mean  is 
greater  than  the  median. 

The  function, 

tp-l  e-tA, 

f(t)  =    ri  ;  p       ■  (3) 


is  an  example  of  an  analytical  function  that  is  skewed  with  a  mean  greater 
than  the  median  for  different  values  of  the  two  parameters  p  and  If   (5).   It 
is  easy  to  observe  that  when  1)   —  2,  the  above  gamma  distribution  reduces  to 
the  well-known  /C     frequency  function  with  2p  degrees  of  freedom.  For  the 
number  of  observations,  n,  greater  than  5,  this  continuous  distribution 
closely  approximates  the  discrete  frequency  distribution, 

iz  =J2{0±  ~Hif    '  (4) 

(=1 

where  k  is  the  number  of  pairs  of  frequencies  to  be  compared,  0^   and  e^  denote 

the  i  th  pair  of  observed  and  expected  frequencies,  and  V^  o .  = /e-  =  n. 

The  number  of  degrees  of  freedom  is  k-1.0  (5)« 

z 

The  %    statistic  is  a  test  to  determine  how  well  certain  data  fit  a 

7. 

known  or  given  hypothesis  (ll).  K.  Pearson  established  the  distribution  of  jC 
for  any  number  of  degrees  of  freedom,  making  it  possible  to  test  whether  a 
sample  came  from  an  assumed  distribution  or  not.   This  made  it  possible  to 
calculate  the  probability,  P,  that  a  particular  value  of  %  would  be  exceeded, 
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The  commonly  accepted  value  for  P  is  0.05,  or  the  five  percent  level  of  accepta- 
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bility.   If  the  observed  value  of  /i    is  less  than  the  value  from  the  distri- 
bution for  the  given  degrees  of  freedom  and  P,  the  observed  value  is  not 
suspected  of  not  representing  a  sample  from  the  distribution  hypothesized. 
If,  however,  the  observed  value  is  greater  than  or  equal  to  the  y£   value 
from  the  distribution,  the  probability  of  this  occurring  is  less  than  0.05, 
and  the  hypothesis  being  tested  is  rejected,  or,  if  the  hypothesis  is  known  to 
be  correct,  then  the  observed  data  is  judged  to  be  nonrepresentative  of  the 
known  hypothesis. 

1+.3  Analysis  of  the  Residence  Times  Distribution 

A  frequency  distribution  of  the  residence  times  in  a  typical  data  set 
gave  the  frequency  function  f(t)  versus  the  time  t  illustrated  in  Fig.  33* 
The  two-parameter  gamma  distribution  was  used  to  develop  the  T  statistic, 
analogous  to  the  ^   statistic,  to  determine  the  data  that  were  not  repre- 
sentative of  the  f(t)  distribution.   It  was  concluded  that  the  experimental 
times  obtained  at  the  right-hand  end  of  the  distribution  in  Fig.  33  were  not 
representative  of  the  residence  times.   This  was  primarily  due  to  errors  in 
correlating  the  recorder  trace  to  the  digital  printer  output.  When  the  radio- 
active source  passed  the  two  detector  levels  at  high  speed,  one  of  the  re- 
sulting pulses  may  not  have  been  strong  enough  to  trigger  the  multi- scaler. 
The  counts  accumulated  until  the  multi- scaler  was  again  triggered  would  not  be 
representative  of  the  residence  times  distribution  but  would  be  far  greater. 
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The  data  points  that  fell  within  the  last  five  percent  of  the  area  under 
the  f(t)  distribution  were  rejected  on  the  "basis  of  the  criteria  mentioned 
above.   The  area  under  the  total  gamma  distribution  integrated  from  zero  to 
infinity  is  1.0.  The  area  under  the  last  five  percent  of  the  distribution  is 
0.05  and  is  equal  to  the  probability  P.   The  T  statistic  was  determined  by 

fT  tP-i  e-v* 

-  =   1-P   .  (6) 


;o    T(p)-yp 


In  each  data  set,  all  of  the  residence  times  that  were  greater  than  the 
T  statistic  were  rejected. 

An  IBM-1620  FORTRAN  computer  program  was  developed  to  determine,  for 
various  values  of  the  distribution  parameters  p  and  V   ,  the  area  under  the 
gamma  distribution  integrated  from  0  to  T.   The  computer  program  is  discussed 
in  Appendix  G.  A  Gaussian  quadrature  numerical  integration  technique,  dis- 
cussed in  Appendix  H,  was  utilized  to  determine  the  area  under  the  distri- 
bution for  combinations  of  p,z>  ,  and  T.   The  Gaussian  quadrature  approximation 

of   (f(x)dx    is 
J.  1 

rl  N 

^  f(x)dx  =     ^T[ak  f(*k)     •  (?) 


A  transformation  of  variables  from 


\  f(t)dt    to   \  f(x)dx  , 

0  J-i 

where  f(t)  represents  the  gamma  distribution,  gives 

/Jt>«  .  j^f-iZ   a^-lP-1  e"^-1* 

°  /(p)  vv  ^zj 

=  A   .  (8) 
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The  computer  program  gave  the  value  of  A  for  various  values  of  p,  "V   ,  and  T. 
Interpolation  of  the  answers  obtained  from  constant  values  of  p  and  1/  gave 
the  time  T  required  for  the  area  A  to  equal  0.95« 

Figure  3^  gives  the  results  of  the  T  statistic  corresponding  to  an  area 
(A)  under  the  curve  j  f(t)dt  equal  to  0.95,  for  given  gamma  distribution 
parameters  p  and  1)  . 

The  parameters  p  and  f  can  be  determined  from  the  experimental  data 

_  2 

by  utilizing  the  average,  t,  and  variance,  q-"    ,  of  the  residence  times. 


_     ft  f(t)dt     vvtA  o  4P^  jtPV^dt 

+  -  **-  Jo  .  (9) 
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O  Jn 


-,+P  -t/V 

t  =  pt;  -  2^ =  pz;  .         (10) 
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-^o Jo 

_{f(t)dt  Jf(t)dt 


*•*=""•  =^ SB — •        (ID 


From  the  definition  of  t,  and  since  T  is  a  constant, 
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Integrating  by  parts  gives, 


a       -ft 
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(13) 


The  average  residence  time  obtained  from  the  data  is 
n 

A—  n 

The  corresponding  unbiased  estimate  of  the  variance  is 


(14) 


n 


t-a2 


s2=\   (*i  -  *) 


<r  / 


n-1 


(15) 
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where  t.  is  the  i—  residence  time  and  n  is  the  number  of  residence  times  in 
each  set  of  data. 

The  most  probable  residence  time  is  determined  when 

d  f(t)  -  ,   x 

—  =  °  •  (16) 

This  is  equivalent  to 


P(ph,p    dt 


_*  (tP-i  e-Vv  j  = 


or, 


«  f(t)  =   e'W    (p-l)tp-2  -  fli  e-V-^  =  0 


e"^  tP  [(p-l)t-2  -  t"1^"1]  =  o  . 


Solving  for  the  most  probable  residence  time,  denoted  by  t  ,  one  has 

XT 


tp  =  (p-l)-^   .  (17) 

An  IEM-1620  FORTRAN  computer  program  was  developed  to  determine  the  aver- 
age residence  time,  t,  the  most  probable  residence  time,  t  ,  and  the  variance 
of  the  residence  time  data  for  each  water  flow  rate  and  distance  between  detec- 
tor levels.   The  program  also  calculated  the  gamma  distribution  parameters  p 
and  V  ,  determined  from  the  average  time  and  the  variance.   Once  the  gamma 
distribution  parameters  p  and  ~U   were  known,  Fig.  3^-  was  utilized  to  obtain 
the  T  statistic  for  the  data  set.  All  of  the  time  increments  in  the  data 
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set  that  were  greater  than  T  were  discarded,  and  the  data  were  rerun  through 
the  computer.  A  new  pair  of  gamma  distribution  parameters  was  obtained,  and 
a  new  T  statistic  was  determined  from  Fig.  J,k.     The  procedure  was  repeated 
until  none  of  the  residence  time  data  were  greater  than  the  T  statistic. 
The  average  time,  the  most  probable  time,  the  variance,  p,  and  V  ,  computed 
by  the  final  pass  through  the  computer,  were  tabulated  for  each  set  of  resi- 
dence time  distribution  data. 

The  larger  the  number  of  data  taken,  the  closer  the  sample  mean  value 
will  be  to  the  population  mean.  By  observing  a  plot  of  the  average  sample 
mean  versus  the  number  of  data  taken,  one  can  estimate  the  minimum  amount  of 
data  needed  to  give  a  sample  mean  representative  of  the  population  mean. 

Figures  35  and  36  give  representations  of  the  average  residence  time  as 
a  function  of  the  number  of  data  values  used  to  compute  the  average.  Appendix 
E  gives  a  discussion  of  Computer  Program  No.  5  used  to  calculate  the  average 
residence  time  as  a  function  of  the  number  of  data  values  used. 

Figures  37  and  38  give  sample  representations  of  the  variance  of  the 
residence  time  as.  a  function  of  the  number  of  data  values  used  to  compute  the 
variance.  Appendix  E  gives  a  discussion  of  Computer  Program  No.  6  used  to 
calculate  the  variance  as  a  function  of  the  number  of  data  values  used. 

Figures  39  through  h6   give  sample  representations  of  the  residence  times 
data  superimposed  upon  the  computed  gamma  distribution  using  the  parameters 
p  and  V   determined  from  the  experimental  data. 

Figures  47  through  56  give  the  gamma  distribution  parameters  p  and  T/ 
as  functions  of  the  distance  between  detector  levels. 

The  residence  time  data  used  in  the  final  pass  through  the  computer  were 
also  used  for  a  FORTRAN  program  (Computer  Program  No.  8)  to  compute  the  average 
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Number      of      Data      Values,      Data      Set      No.  2. 


68 


0.12 


V) 

c 
o 
o 
0 


E 


o 
o> 

c 

v_ 
o 

> 

< 


.                i 

T" 

1 

'                1        ..  .    1 

0.1  1 

- 

r            g 

,  o 

— © 

■■ 

- 

O.IQ 

• 

0.09 

. 

c 

> 

• 

Flow      Rate  =  B3 
A  2   =  2     Inches 

- 

0.08 

- 

1                        1 

1 

1 

I                  1                  | 

- 

10 


20 


30 


40 


50 


60 


70 


Figure;     36. 


Number     of       Data      Values 

Average       Residence      Time      as      a       Function      of 
Number       of      Data      Valuesr     Data      Set      No.     37. 


the 


69 


0.005 


0.004  . 


0.003  . 


0.002  L 


0.001  - 


0.000 


— r r               -i- 

i 

i                    i 
4.0  in,    Hg 

« — 

• 

AZ    =  V2     inch 

I                    i                    • 

Z  ,  to  2  2 

- 

■©■ 

u        *- 

<                    j                    i 

,   ,  i     - _.,.  i  

i 

10 


60 


70 


2  0  30  40  50 

Number    of    Data    Values 

Figure    37.    Variance    of    the    Residence    Times     Data    as     a    Function 
of    the    Number    of     Data     Values,    Data    Set    No.    16. 


70 


0.13 
0.1  2 

0.11 
0.10 
0.09 
0.08 
0.07 

0.06 
005 
0.04 
0.03 


; 


4.85  in.  Hg 


AZ    =   4  inches 


2  |  to  Z  2 


J. 


10  15  20  25 

Number      of     Data    Values 


30 


35 


Figure    38.    Variance    of    the    Residence.   Times    Data    as    a    Function 
of   the    Number    of    Data    Values,    Data   Set    No.   42. 


-   15 


2.5  In.  Hg 
AZ=    2  inches 
Z2 to  Z, 


2.4  3.2 

t  ( seconds) 


4.0 


4.8 


5.6 


Figure    39.    Frequency    Distribution     of     the    Residence    Time    Data 
and    Theoretical    Gamma     Distribution,    Data    Set   No.  8. 


72 


20 

I  S 

16 

14 

12 

10 

8 

6 

4 

2 

0 


«• 

I 

t 

4.0  in.  Hg 

A  Z  =    2  inches 

- 

22   fo  Zj 

•" 

• 

_ 

- 

o   / 

o 

°\ 

o 

- 

- 

X      o 

- 

- 

^o 

— 

L_ 

•°^^     o      c^- 

o 

o 

0.5 


t        (seconds) 


1.5 


Figure     40. 


Frequency     Distribution     of    the     Residence   Time    Data 
and  Theoretical    Gamma    Distribution,    Data    Set    Mo.   22, 


l 

T~ 

9 

- 

4.0  in.   Hg 
£Z=    4   inches 

- 

8 

Za  toZ, 

— 

7 

- 

- 

6 

_ 

/      ° 

-    5 

- 

- 

4 

- 

°\ 

- 

3 

- 

2 

- 

- 

I 
0 

Oy 

/Q 

i 

^ 1 

- 

1.0 


2.0 


t  (seconds  ) 


3.0 


Figure    41. 


Frequency     Distribution    of    the    Residence    Time     Data 
and    Theoretical    Gamma     Distribution,     Data    Set    No.   27. 


lh 


0.5 


4.85    ia    Hg 
AZ  =  4  inches 
21  toZi 


i.O 

t   (seconds) 


2.0 


2.5 


Figure     42.    Frequency    Distribution    of    the    Residence    Time     Data 

and    Theoretical    Gamma    Distribution,    Data    Set    No.  44. 


75 


1 

l" 

" V  '                           » 

36 

•m 

7.0  in    Hg 

- 

32 

A  Z  a  4inches 

- 

28 

Z2  to  Z, 

- 

24 

- 

20 

m 

• 

- 

16 

- 

/     ° 

°\ 

Vo 

- 

12 

- 

- 

8 

• 

- 

4 

o  ' 

ON. 

- 

0 

.!_ 

j,                 ^* 

j                           i 

0.5 


1.0  1.5 

t    (seconds) 


2.0 


2.5 


Figure    43.    Frequency     Distribution     of   the    Residence    Time    Data 
and    Theoretical    Gamma    Distribution,    Data    Set    No.  59. 


76 


20 


.    ,  r, 

— 1 

■     J-       -  — 1 

i  "             •  \ 

i 

18 

„ 

4.0  in.  Hg 

- 

16 

_ 

AZ  =    4  inches 

wm 

14 

- 

Z  |  to  Z2 

- 

12 

— 

10 

0 

- 

8 

°\ 

wm 

6 

V  ° 

- 

4 

-      / 

- 

2 

~  / 

'Xo     0 

- 

0 

1 

1 

.,  ,  i   T^r        : 

i             -      i 

' 

1.0 


3.0 


3.5 


2.0 
t  (seconds) 

Figure    44.    Frequency    Distribution    of    Residence    Time    Data 

and    Theoretical    Gamma    Distribution,     Data    Set    No.  28, 


T~ 

I 

7.0    in.    Hg 

4 

c 

LZ  -     2    inche  s 
Z,    toZ2 

— 

3 

- 

3  2 

<•- 

Q 

o\ 

- 

1 

n 

o 

i 

—  Oi 

i 

- 

0  0.5  1.0  1.5 

t    ( seconds  ) 

Figure     45.     Frequency     Distribution      of      the  Residence     Time      Data      and 

Theoretical     Gamma      Distribution,  Data      Set     No.    53. 


Tci 


20 


1 

1 

i  £ 

7.0  in.  Hg 

- 

16 
14 

- 

A  2  =    4  inchos 
Z\    to  2 2 

- 

12 

- 

10 

- 

8 

/       ° 

- 

- 

6 

, 

- 

4 

/ 

0& 

- 

2 

™      /"^ 

_ 

0 

1 

i 

1.0 


2.0 


3.0 


t  (seconds) 


Figure    46.    Frequency    Distribution    of    the    Residence    Time    Data 
and    Theoretical    Gamma    Distribution,    Data    Set    No.  58. 


19 


0.8  _ 


_  0.6 

•o 

■o 
c 
o 
o 
c> 
w 


*  0.4  . 


0.2  . 


0.0 


i              r* 

T' 

"1 

t 

i 

i 

- 

o 

Z,  to  Z2 
Z2  to  Z, 

- 

- 

□ 

^     Q 

- 

- 

- 

^     G 

- 

- 

B^^ 

Flow 

Rate  =  Bj 

- 

I                   i 

1 

• 

i 

t 

! 

3  4 

A  Z  (inches) 


Figure     47.      The      Gamma      Distribution      Parameter     v,    at     Flow     Rate     B|. 


to 


1 

II 

i 

"T~ 

.                     i                     i 

ao 

• 

B 

O 

O       2,  to  22 
□       22  to  2, 

— 

6.0 

o\ 

o^\ 

a 

mi 

4,0 

© 

© 

- 

B 

2.0 

Flow    Rote  =  Bj 

- 

0.0 

j 

i 

i 

1 

!  _                ...  1—                _J. 

0  12  3  4 

A  2  (inches) 

Figure     48.     The      Gamma      Distribution 
Flow       Rate       Bj. 


Parameters     P,     at 


~  0.15 

10 
"O 

c 
o 
o 

<D 

0> 


S  o.io. 


0.05. 


0.0 


1                     > 

i 

i 

i                   , ...... 

"■     1 — 

1 

O     2,  to  Z2 

□     Z2  to  Z| 

© 

- 

- 

0 

- 

^^"E 

0 

h    >;     0 

A> 

Flow   Rate  =  B2 

- 

^0 

i i     _ 

! 

t 

1                  » 

I 

0  I  2  3  4 

A  Z  (inches) 
Figure     49.     The      Gamma      Distribution 

Rate      B2. 


5  6 

Parameter     v% 


at     Flow 


62 


180 

160 

K0 

12.0- 

100 

80. 

6j0 
4.0 

2.0 

0.0 


V 

i                     i 

r~ 

1 

( 

i 

. 

n 

0     Zi  to  Z2 

_ 

- 

\° 

□     Z2  to  Z, 

- 

- 

□ 

- 

. 

El 

w 

- 

0 

^X® 

H 

G 

- 

- 

O 

Flow 

Roto   =    B£ 

- 

i 

>                   i 

i 

■  i 

• 

i 

3  4 

A  Z    (inches) 


Figure      50.     The     Gamma      Distribution      Parameter     P,     q^ 
Rate      B2« 


Flow 


63 


i 

i                      t 

* 

i 

i^ 

r» 

0.16 

o 

z2 

to 
to 

z 

2 
1 

- 

0.14 

- 

- 

£0.12 

c 
o 

^^^^^"            ^/ 

o  0.10 

0> 

^*^lL 

- 

*  0.0£ 

Q      >^ 

- 

0.06 

*© 

- 

0.04 

- 

a 

• 

0.02 

- 

& 

Q 

Flow 

Roto  =  B3 

- 

0.0 

J- 

-J ! 

, 

i 

t 

< 

Figure     51.      The      Gamma 
Rate      B3. 


14  5  6 

A  Z  Cinches) 

Distribution      Parameter     v. 


at     Flow 


bh 


• 

i 

t 

r" 

i                  r™ 

r™ 

«, 

o 

Z,    toZ2 

_ 

□ 

Z2  to  Z| 

- 

a 

• 

- 

•m 

©\ 

El 

B 

- 

- 

0 

<~-~^T 

Flow    Rate  =  B3 

- 

mm 

• 

I 

' 

j 

i 

»                         1 

' 

16.0 

14.0 
12.0 
10.0, 

8.0- 

6.0 
4.0 

2.0 

0.0 

0  12  3  4  5  6  7 

A  Z  (inches) 

Figure     52.     The     Gamma     Distribution     Parameter        P,      at      Flow 
Rate      B3. 


&5 


1 

i 

T 

i 

B 

1 

o.os 

- 

B 

- 

01 

o  0.07 
o 

- 

B 

0/ 

^ — © 

© 

- 

^ 

- 

•  0.05 

- 

© 

^     H 

0    Zj    to  Z2 
□     Z  2   to  Z  | 

- 

0.03 

1*  J 

• 

- 

© 

Flow 

Rate    =    B4 

- 

0.01 

j 

I 

,.  %  

i 

1...              -.i 

_      .!. 

)  I  2  3  4  5  6  7 

A  Z  (inches) 
Figure     53.    The     Gamma     Distribution     Parameter-   2/,     at      Flow 
Rate     B4. 


&6 


18.0 

16.0 

14.0. 

12.0 

10.0 

8.0 

6.0- 

4.0;. 

2.0 


T 


0.0 


o     Zi  ro  Z2 
□      z  2  to  z  J 


Flow     Rate  =  B4 


0  1  2  3  4  5  6  7 

A  Z   (inches) 

Figure     54.     The     Gamma     Distribution     Parameter     P,     at      Flow 
Rate     B4. 


67 


w_  w 

™i 

~"T~"                      f 

11 

i                   i 

» 

0.18 

O 

2,  <to  Z2 

'- 

0.14 

□ 

2  2  to  Z  j 

c    0.12 

o 

o 

5  o-io 

- 

& 

0.08 

- 

n 

m^^ 

- 

0.06 

- 

Q 

o 

- 

0.04 

- 

^Gj 

- 

0.02 

'if 

Flow 

Rate  =  B5 

- 

0.00 

L. 

»                    < 

» 

i                   > 

i 

Figure 


12  3  4 

A  Z  (inches) 

55.       The      Gamma      Distribution 
Rate      B5. 


Parameter      vy      at      Flow 


66 


i 

•                     i 

j 

l 

~ 

i 

24.0 

O  2|  toZ2 
E   Z2  to  Zj 

- 

21.0 

'  E 

- 

ia.o 

• 

- 

- 

15.0 

• 

E 

Flow 

Rate  =  B5 

- 

12.0 

. 

- 

9.0 

□ 

- 

6.0 

© 

G 

0 

-©— 

- 

3 

3.0 

- 

0.0 

i 

i                        i 

_  i 

i 

! 

1 

3 

L  Z 


incne-s 


Figure 


56.      The 
Rate 


Gamma      Distribution      Parameter      P,      at      Flow 

Br. 


89 


speed  of  the  tagged  particle  for  each  set  of  data  in  the  following  manner: 


=       N       -  (18) 


Z4Z  AZ      Y~   JL_ 

tj_  n  n      ^—      t.* 


The  output  of  Computer  Program  No.  8  is  given  in  Appendix  F. 

The  speed  of  the  tagged  particle  for  each  distance  between  detector 
levels  and  water  flow  rate  was  determined  by  the  following  three  methods: 


( v  '  ^7   =  I^iiv    •  (19) 


tT\  AZ  AZ  n  ,      . 

c'=i 


y\ 


v      = 


4Z 


*       V~"~    1 

~2_T7     '  (18) 


The  quantity  v(t )  is  the  speed  of  the  tagged  particle  corresponding  to  the 
most  probable  residence  time,  v(t")  is  the  speed  of  the  tagged  particle 
corresponding  to  the  average  residence  time,  and  v  is  the  average  speed  of 
of  the  tagged  particle. 

Figures  57  through  71  give  the  results  of  v(t  ),  v(t),  and  v"  as  functions 
of  the  distance  between  detector  levels. 

Figures  72,  73>  and  "Jk   give  straight-line  extrapolations  to  A  Z  —  0 
values  (obtained  from  Figs.  57  through  71 )  for  the  speed  of  the  tagged  par- 
ticle as  a  function  of  the  mean  fluid  velocity  through  the  fluidized  bed. 
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Figure     73.     Particle     Speed     Corresponding     to     the     Most 
Probable      Residence     Time     as     a     Function    of 
the     Mean      Fluid     Velocity. 
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h.k     Analysis  of  the  Particle  Traffic 

The  number  of  times  a  particle  passes  a  level  Z^  per  unit  time  is  de- 
fined as  the  particle  traffic,  r(Z  ),  and  was  determined  by  counting  the 
number  of  times  the  tagged  particle  passed  each  detector  level  in  each  data 
set  and  dividing  by  the  total  time  required  to  obtain  the  data  as  described 
previously. 

An  IBM- 1620  FORTRAN  computer  program  was  developed  to  determine  the 
particle  traffic,  and  it  is  discussed  in  Appendix  3. 

Figure  75  is  a  plot  of  the  particle  traffic  data  obtained  from  each 
detector  level  as  a  function  of  the  mean  fluid  velocity  through  the  fluidizing 
column. 

U.5  Analysis  of  the  Reflectivity  Data 

The  reflectivity,  Pr,  was  previously  defined  as  the  probability  of  a 

particle  passing  level  Z.  and  failing  to  reach  a  different  level  (Z.)  a  given 

distance  (AZ)  away  from  Z.  before  returning  to  Z  .  This  probability  is 

1  i 

represented  by  Pr(Z.,Z.;Z.)  and  is  equal  to  (m-2)/m,  where  m  is  the  number 

-1-  j  1 

of  times  the  particle  passed  Z.  before  reaching  Z..   The  reflectivity  is  equi- 
valent  to  the  reflection  coefficient  or  albedo  used  by  some  authors  (3)« 

The  transmissivity  is  the  probability  of  a  particle  passing  level  Z^  and 
reaching  a  different  level  (Z.)  a  given  distance  (Z\Z)  away  from  Zi  before 

returning  to  Z . .   This  probability  is  represented  by  Pr(Z.,Z.;Z  ).   The  sum 

1  1     J     j 

of  the  reflectivity  and  transmissivity  probabilities  must  equal  1.0;    therefore, 
■Pr(Z±,Zy,Z±)   si-  Pr(ZltZj|Zj)      .  (21) 
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The  reflectivity  and  transmissivity  probabilities  are  easily  derived  by 
noting  that  for  m  particles  passing  Z. ,  m/2  of  the  particles  are  moving  in 
the  direction  toward  Z.,  and  one  particle  eventually  reaches  Z..   The  pro- 
bability  of  transmission  to  Z.  is  l/(m/2)  or  2/m.  The  probability  of  re- 
flection  is  1  -(2/m)  or  (m-2)/m. 

The  following  analysis  was  used  to  attempt  to  express  analytically  the 

reflectivity  as  a  function  of  the  distance  between  detector  levels  (AZ)   and 

the  position  of  the  lower  detector  level  (Z  ). 

i 

The  probability  of  a  particle  passing  Z.    and  reaching  Z.  before  return- 

J 

ing  to  Z.  has  been  represented  by 

Pr(Zi,ZJ.;Z;j)  . 

Generalizing,  the  probability  of  a  particle  passing  Z.  and  reaching  Z  before 
returning  to  Zj_  is 

Pr(Zi,Z;Z)  . 

The  probability  of  a  particle  passing  Z  and  reaching  Z+AZ  before  returning 
to  Z  is 

Pr(Z,Z+AZ;Z4.<AZ)    , 

and  the  probability  of  a  particle  passing  Z.  and  reaching  Z+4Z  before  return- 
ing to  Z.  is 

Pr(Zi,Z+4Z;Z+AZ)  . 
These  two  probabilities  are  related  by 

Pr(Zi,Z+AZ;Z+ZZ)  =  Pr(Zi,Z;Z)Pr(Z,Z+AZ;Z+AZ) .    (22) 
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The  probability  of  a  particle  passing  Z.  and  reaching  Z  but  not  Z+AZ  be- 
fore returning  to  Z.  is 

Pr(Zi,Z;Z)[l-Pr(Z,Z+/lZ;Z+AZ)J  . 

Assuming  the  relationship 

l-Pr(Z,Z+AZ;Z+4Z)  =  fi  AZ      ,  (23) 

where  J!>   is  defined  as  the  reflection  coefficient,  and  substituting  it  into 
Equation  22  gives 

PrCZ^ZtdZ^+AZ)  =  Pr(Zi,Z;Z)(l-p/\Z)  .         (24) 

Rearranging  gives 

Pr(Zi,Z*4Z;Z.MZ)-Pr(Zi,Z;Z)  =   Qpr(z.>Z;Z) 
A  Z 

Taking  the  limit  as  AZ   — >•  0, 


dZ 

Integrating  gives 


dPr(Z*'Z;Z)  =   ^Pr(Zi,Z;Z)  . 


dPr(Zi'2f  -  .«  (    dZ   , 
Pr(Zi,Z;Z)      T  \ 


'^V 


In  Pr(ZifZ;Z)  r  -j?  (Z-Z^  +  K  (25) 


or, 


'r{Z±9Z;Z)    -    K  e")5^"2!^ 


(26) 
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When  Z  =  Z±,   Pr(ZjL,Z;Z)  =  1.0.   Therefore, 

Ksl. 

Z  is  limited  to  the  height  of  the  fluidized  "bed,  L.   The  probability  of  trans- 
mission from  Z.  to  L, 
1     ' 


Pr(ZifLfL)  =  e-^(L"Zi)   , 


is  never  zero  for  the  range  of  experiments;  therefore,  Pr(Z^,Z;Z)  ranges  from 
e"P^  "  i'  to  1.   The  probability  of  reflection,  therefore,  ranges  from  zero 
tol-e-^-Zi). 

Figures  "J6,   TJ,   and  78  give  illustrations  of  the  frequency  distribution 
of  Pr(Z^,Z.;Z. )  for  different  data  sets,  where  f(Pr)  is  the  number  of  occurrences 
of  each  value  of  Pr(z.  ,Z.;Z. ) .   Only  discrete  values  of  the  reflectivity  occur, 
since  the  reflectivity  is  defined  as  (m-2)/m,  where  m  represents  integers 

2,3)"',   etc.  Note  that  the  frequency  distributions  of  Pr(Z.,Z.;Z.)  compare 

1  J  -1- 

very  closely  for  Zj_  being  either  the  upper  or  lower  detector  level. 

The  average  m  for  each  set  of  data  was  determined,  and  the  average 
Pr(Z.^,Z ,;Zi)  for  each  set  of  data  was  defined  as 

Fr(Z19Zj;Z1)  =  (m-2)/G  .  (27) 

A  least  squares  analysis  was  used  to  fit  the  analytical  expression  Pr(Z  ,Z.;Z.) 
=  1  -  e~  fi(zyzi)   to  the  experimental  data,  Fr(Zi,Z.;Zi) .  The  following 
method  was  used. 

Since  the  probability  of  transmission  is  given  by 

Pr(ZltZj|Zj)  =  e-^(ZJ  *  Zi)   , 
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the  average  probability  of  transmission  determined  from  the  experimental  data 
is 

pHz^jZj)    =  1   -  v5{Z±,Zy,Z±)    .  (28) 

To  simplify  the  nomenclature,   let 

SrtZ^ZjiZj)    ■    P?     .  (29) 

For  each  water  flow  rate,  the  difference  between  experiment  and  correlation 
for  each  data  value  was  computed.  Forming  the  sum  of  the  squares  of  the 
deviations  gives  a  measure  of  how  well  data  values  fit  an  analytical  expres- 
sion as  a  whole.   Therefore,  the  sum  of  the  squares  of  the  deviations  of  the 
experiment  from  correlation  for  each  water  flow  rate  is 


all  data 


Z      [e-M^i).^]2  f 


k-  1 


(30) 


or 


6 


The  parameter  R    should  be  chosen  so  as  to  minimize  the  sum  of  the  squares 
of  the  deviations.  To  do  this, 


j|    £r  (_Mzk  -  in  Pr/]  =  o  . 


(32) 


Ill 


Solving  for  the  best  estimate  of  jQ    gives 


t     '       '         ■>'  i\  ■    ■       ■■■■!■ 

J3  2^-^  +    z_ 


2 

AZ^  In  Prk     =     0 


-  zz 


^  =    -=£= .  (33) 

Figure  79  presents  the  average  reflectivity  data  as  a  function  of  the 
distance  between  detector  levels  and  water  flow  rate,  along  with  the  theo- 
retical curves  determined  from  the  least  squares  analysis  of  the  data  to 
determine  the  best  estimate  of  R   .   The  lower  detector  level  was  Z. . 

Figure  80  presents  the  theoretical  reflectivity  curves  normalized  to 
the  reflectivity  at  the  top  of  the  fluidized  bed  as  a  function  of  the  dis- 
tance from  the  bottom  of  the  fluidized  bed  for  a  particle  initially  passing 
the  lower  detector  level. 

Table  5  presents  the  values  of  B   determined  from  the  least  squares 
analysis  of  the  reflectivity  data  for  the  water  flow  rates  studied. 

Table  5«   Reflection  Coefficient  determined  from 
the  least  squares  analysis  of  the  re- 
flectivity data. 

WATER  FLOW  RATE  fi 

B]_  (13-5  gal/min.)  0.300/in. 

B2  (17-0  gal/min.)  0.125/in. 

B^  (18.8  gal/min.)  0.087/in. 

B^  (22.3  gal/min.)  •  0.088/in. 

B^  (2^.0  gal/min.)  0.100/in. 
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k.6     Analysis  of  the  Particle  Concentrations 
4.6.1  Theoretical  Development 

A  knowledge  of  the  concentration  of  particles  with  respect  to  time  and 
position  for  the  given  source  condition  and  "boundary  conditions  is  of  par- 
ticular interest  in  the  study  of  fluidized  beds.   If  the  assumption  is  made 
that  the  concentration  of  particles  can  be  represented  by  a  solution  of  the 
diffusion  equation  applying  the  source  and  boundary  conditions  of  the  fluidized 
bed,  a  dispersion  coefficient  based  upon  the  experimental  data  can  be  obtained. 

The  diffusion  equation  for  no  absorption  may  be  written 

E  V2C  +  Q  =  dC/^t,  (34) 

(  -  Leakage  +  Source  =  Time  rate  of  change  of) 

Particle  Density 

where  E  is  the  diffusion  coefficient.  When  the  diffusion  is  one-dimensional 
this  reduces  to 


;2 


B 


^o   +  <*  =*C/*t,  (35) 


2  bz2 


where  Ez  is  defined  as  the  dispersion  coefficient.   The  Z-dependent  diffusion 
equation  must  be  solved  subject  to  the  following  boundary  conditions  and 
source: 


1£ 


«   0,  (36) 

Z  =  L 
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dz 


=  0, 


z  =  o 


(37) 


Ct  =  0  =  °' 


(38) 


at  t  =  0,  Q  particles  per  in  are  injected  at 

z  =  z0,  =  q  i(z-z0)  i(t) 


The  method  used  to  solve  the  diffusion  equation  for  the  given  source 
and  boundary  conditions  is  outlined  as  follows: 

Taking  the  Laplace  transform  with  respect  to  time, Equation  35  gives 


E. 


d2C 


2   dZ2 


-  sC   =   -q  $(Z-Z0)  . 


(39) 


Letting 


E. 


*2, 


and  substituting  Equation  40  into  Equation  39  gives 


(40) 


d2C 


dZ 


2   "   X*  = 


-  -1-  i(z-ZQ)  . 

Ez 


(41) 
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Taking  the  Laplace  transform  of  Equation  41  -with  respect  to  Z  yields 

p2C   -   pC(0)-    C'(0)-A2G     =    -   JL  e"PZo    .  (42) 

E„ 


Applying  Equation  36,  and  solving  C  gives  Equation  43. 


J  =  PC(O)   _  0,  e-^Zo 
P2-^   '  Ez(p2-X2) 


(43) 


The  first  term  on  the  right-hand  side  of  Equation  ^3  represents  the  homo- 
geneous or  steady- state  term.  The  homogeneous  solution  of  Equation  4-3  is 
given  "by  Equation  hk. 


C  =  A  eAZ  +  B  e~*Z 


(44) 


The  inverse  Laplace  transform  of  the  second  term  on  the  right-hand  side  of 
Equation  V3  is  determined  "by  utilizing  the  Second  Shifting  Theorem  (12). 
Any  function  f(x-a)  multiplied  by  H(x-a)  has  the  value  of  zero  for  x<a  and 
f(x-a)  for  x>a.   The  Second  Shifting  Theorem  is 


JL 


-1 


']■ 


e"as  f(s)   =  f(x-a)  H(x-a), 


(45) 


where  H(x-a)  is  the  unit  step  function  defined  above.  Applying  this  theorem, 
one  has 


Z 


-1 


Q  e~Pzo 


Mp2-*2) 


L  z 


Q   *-l 


E, 


£- 


>2-^2 


H(Z-Zft)  .   (46) 
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The  inverse  Laplace  transform  of  l/(p  -  A2)  with  the  argument  (Z-ZQ)  is  given 
"by  Equation  kj . 


t- 


>2-*2 


sinh  A(Z-Z0)   . 


(47) 


Substituting  Equation  ^7  into  Equation  k6  gives 


£ 


-1 


q   e"Pzo 
E,(p2-A2) 


=  -i-  sinh  a(Z-Z0)H(Z-Z0)    .  (48) 

*E„ 


Noting  the  identity 


sinh  A(Z-ZQ)  = 


*(z-z0)_  e-*U-z0) 


and  substituting  into  Equation  kS  one  can  obtain 


^ 


-1 


Q   e"Pzo 


Ez(p2-A2) 


Q 


2AE. 


:^Z-ZoLeA(Z-Z0) 


H(Z-Zft).    (49) 


For  Z>Z0,  only  the  second  term  of  Equation  49  is  applicable,  since  the  con- 
centration of  particles  beyond  the  fluidized  bed  height  (L)  is  zero.  Equation 
50  below  represents  the  solution  of  the  second  term  in  Equation  4  3  for  Z>  Z  . 


£ 


-1 


Q  e-Pzo~ 
^(p2-*2) 


_i-  e-^Z-Zo)H(Z.Z0)  . 
2AS2 


(50) 
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Since  the  particles  injected  at  ZQ  will  also  move  downward,  the  solution  of  the 
second  term  of  Equation  k-3   must  also  account  for  the  downward  movement.   The 
solution  is  the  same  as  for  Z>  Z  .   The  equation  applicable  is,  therefore, 
Equation  50,  except  the  Z  and  Z  are  interchanged  in  the  exponential  function, 
and  the  unit  step  function  H(Z-Z  )  is  changed  to  H(Z  -Z)  for  its  proper  oper- 
ation when  Z<  Z  .   Equation  51  "below  is  the  result. 


£ 


-1 


Q  e 


-pz0 


Ez(p2-*2) 


2^E 


4    -X(Z0-Z) 

—  e         H(Z-Z).   (51) 


The  general  solution  for  C  is  the  sum  of  Equations  kk,    50,  and  51* 


C=Ae  +Be   +  -— 


2XEZU 


->(Z-Zj 


••■   '0'E(Z.^)+e-^Zo-Z)H(Z0-Z)]  . 


(52) 
If  (Z-Z  )  is  changed  to  the  magnitude  of  .Z-ZQ  in  Equation  52,  the  operators 

H(Z-ZQ)  and  H(Z  -Z)  may  he  removed,  giving  Equation  53- 


-     AZ   ^  -*Z 

C  =  Ae   +  Be   + 


Q 


-A|Z-Z0| 


2XE, 


(53) 


The  constants,  A  and  B,  are  determined  by  applying  the  boundary  conditions 
given  in  Equations  36  and  37  as  follows: 


[-1  " 

[dZ  JL 


,  .  AL  „  -XL    Q   -a|L-Z0 

A*  e  -B*e e       =  0  . 

2EZ 


(54) 


UZ  J  0  T  2E„ 


(55) 


Solving  Equations  5^  and  55  for  A  and  B,  one  obtains 
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A  = 


2E, 


3L  e'*Zo 


2S 


-^e 


-A 


-*L 


^L 


-*L 


(56) 


Expanding  the  determinants  in  Equation  56  gives 


AQ 


A  = 


-A(L-Z0)    -ML+Z 


♦  e 


3 


A'(e 


AL    -Hn 


(57) 


-XL 
Factoring  e    in  the  numerator  of  Equation  57,  and  applying  the  following 

identities, 


cosh^ZQ    = 


sinh^L     = 


e^Zo   +    e~*zo 


e-XL      m    Q-XL 


(58) 


one  obtains 


q    e"^L   cosh^Z( 
2EZ  2  sinhAl 


(59) 
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E  = 


}e 


*L    _  Q   «-A(L-Z0) 


23, 


5L  e-*zo 


23. 


2  /\   sinh  !*L 


(60) 


Expanding  the  determinant  in  Equation  60,  and  applying  Equation  58  give 


B  = 


Q,  cosh  A(L-Z0) 
2EZ  3  sinh^L 


(61) 


To  obtain  the  final  solution  for  the  particle  concentrations,  C,  the  inverse 
Laplace  transformation  of  Equation  53  must  be  performed.   To  do  this,  Equation 
53  is  separated  into  the  sum  of  the  complementary  and  particular  parts  as, 


G  ~  Ccomp  +  Cp  * 


(62) 


where 


and 


Q 


'comp 


23 


Z  L_ 


r^(L-Z)cosh;iz   -^Zcosh>(L.z  ) 


^  sinh  XL 


,  (63) 


CP  = 


a-  A(z-zc) 


2XS, 


(64) 


For  the  complementary  part,  note  that  there  are  poles  at  A=0  and  n^f':i 

L 
The  inverse  Laplace  transform  for  the  case  of  simple  poles  is  given  by  the  sum 


of  the  residues  (12).   Thus, 
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n 
f(t)=  \   11m  (s-ak)  f(s)  est  .  (65) 

/     S  >  Q.T- 
K-l  K 

The  relation  of  (s)  to  /{        is  given  by  Equation  ^0,  "but,  for  simplicity,  /^ 
will  still  be  used  in  the  equations.   Operating  on  Equation  63  as  indicated  in 

Equation  65  produces  the  indeterminate  form  0/0  for  ^  -   0       and  n  ,     , 

L 
where 


f(3)=   e-^-^coshaZp  +  e"^ZcoshML-Z0)  ^  (g6) 

7\sinh  AL 

The  application  of  L' Hospital's  rule  will  produce  a  solution.  Let 

7(b)  =  JLLs!   ,  (67) 

D(8) 

where 


N 


(s)   =  e-^(L-Z)cosh^Z0+    e-^ZcoshA(L-Z0)    ,  {68) 


and 


D(s)   =   TK  sinh^L   .  (69) 

Equation  65  becomes 


f(t)r     V"      lim   (s-ak)   Ilil  est   . 
Z s  — »av  D(s) 


(70) 


Differentiating  the  numerator  and  denominator  of  (  s-av)  ■  ■ . 

D(s) 

with  respect  to  s  gives 
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^(S.ak)K(s)/^  B(8)  =    lT(s),(s-ak)II-(S)    < 
ds  ds  Df(s) 


(71) 


Therefore,   Equation  70  "becomes 


tf 


f(t)= 


K-/ 


lim 

s-^ak 


U(s)+(s-ak)K«(s) 


DT(s) 


est    , 


(72) 


where 


D»(s)  =     J^-1.  (Asinh^L) 
ds      dA 


(73) 


From  Equation  kO , 


dA 


ds  2^E, 


(74) 


Therefore,   Equation  73  becomes 


Dl(s)  =    — - —  (sinh^L  +AL  coshAL)    . 
2^EZ 


(75) 


The  residues  at  ^\  =  0  and     ■    are  given  "by  utilizing  Equations  63,  68, 

Xi 

and  75  in  Equation  72. 


At    A  =  Ot 
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'comp 


_    q 


7^-0 


E, 


sinhAL 


4  L   coshXL 


_     Q 


(76) 


where 


sinhAL 


X=0 


_d   sinh^L 
da 

d> 


as  L  cosh^L  —  L   * 


At   *  =  HiLi  , 


CO 


comp 


A  = 


_  ntfi     2E, 


n=  j      L     z 


'  .  £!Ii(L-Z)  ntfZ0 

e       L  cos  s^= — 

j-i 

_lTsinh  nfli       L   cos  nl3 
2EZ  L     n-rri  J 


e** 

«                             ifW    /n            ^Ox 

L          cos  n1T  (1-  -£-) 

1 

^sinh  nltt       Ti  ^     J 

2EZ 

niri                            J 

2   2^ 


where 


(77) 


cosh 


niri 


=  cos 


nff 


Noting  that 


sinh  ntfi    =     i   sin  n1T  =  0    , 


Equation  77  "becomes 


124 


CO 


comp 


a  = 


nffi     L 


M=l 


r    -   niri(L-Z)  n?fZ, 


L  cos 


cos  nit 


-  221  z  /,      Zo    ~         ^2. 

e       l  cos  nir(l-  £  ) 


cos  nTT 


-a. 


e     "T2"  ^z 


Substituting  the  identities 


and 


cos  ntr(l-x)  =   cos  nflr    cos  nffx  , 


.ie 


e   =  cos  £   +  i  sin0 
into  Equation  jQ   gives 


(78) 


(79) 


(SO) 


CO 


COEip 


20. 


n2ff2t 


5=  ±2.  \   cos  nfrz„  cos  HH2.  e"  T2 


*  = 


nua 


r\=i 


(31) 


The  inverse  Laplace  transform  of  the  particular  part,  Equation  6k,    obtained 
from  a  table  of  Laplace  transforms,  is  given  in  the  following  equation: 
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G   =   S e"  4S„  t   #  (82) 

*     2(E-fft)* 


The  complete  solution  for  C  becomes 


c  =   )   ccomp  +  cp  •   ■  (83) 


Defining  the  following  dimensionless  variables, 


X  = 


*0=^  ,  (84) 

and  substituting  them  into  Equations  j6,   8l,  82,  and  83,  and  rearranging,  one 
obtains 

t  ^2 

—   2  o,  (x"Xo) 


-—  =  1  +  2  )   e      cos(ntrx)cos(mrx0)+       -A 


n  =  i 

The  equilibrium  concentration  obtained  a  long  period  of  time  after  intro- 
duction of  the  quantity  of  tracer  Q,  is 


-£-  =   C^  ,  (86) 
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or 


f  -   2i  •  (87) 


4.6.2  Analysis  of  Experimental  Results 

The  concentration  of  tagged  particles,  C,  is  assumed  to  be  proportional 
to  the  number  of  times  that  a  tagged  particle  passes  a  horizontal  plane  per 
unit  time,  I,  after  it  has  been  introduced  into  the  fluidized  bed.  Since  C 
is  assumed  proportional  to  I, 


£_   =   -L.  .  (88) 


The  data  obtained  from  a  single  radioactive  particle  were  used  to 
determine  the  concentration  due  to  Q  particles  being  injected  at  one  time, 
A  tagged  particle  is  considered  to  be  injected  into  the  fluidized  bed 
at  a  detector  level  every  time  the  tagged  particle  is  detected  at  that 
detector  level.  Each  time  a  tagged  particle  passed  one  detector  level  a 
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time  scale  was  started,  and,  every  ti:,ie  it  p.  ;sed  the  other  detector  level, 
the  time  was  recorded.   This  was  done  for  a  finite  length  of  time.   These 
data  were  used  to  obtain  the  number  of  times  a  particle  passed  a  detector 
level  per  unit  time  as  a  function  of  the  time  since  it  was  first  detected 
at  a  detector  level.   Summing  Q  such  occurrences  in  the  set  of  data  will 
result  in  the  same  distribution  as  injecting  Q  particles  at  a  detector  level 
at  the  same  instant. 

Figure  8l  illustrates  a  modified  dual  channel  recorder  trace  for  a 
typical  set  of  data.   The  minus  sign  represents  the  pulses  obtained  from  the 
lower  detector  level,  and  the  plus  sign  represents  the  pulses  from  the  upper 
detector  level.   The  particle  concentration  data,  I,  is  obtained  as  described 
previously.   The  time  span,  W,  in  Fig.  8l,  is  chosen  of  great  enough  length 
to  obtain  the  equilibrium  concentration,  I  .   The  time  span,  W,  is  then 
divided  into  q  increments  of  time  ^t,  giving 

\  4tn-  =  W   . 


The  time  increments  do  not  necessarily  have  to  be  equal. 

The  time  span,  W,  is  begun  initially  with  the  first  pulse  in  the  data 
set,  which  may  be  either  positive  or  negative.   The  number  of  positive  or 
negative  pulses  within  each  time  increment  ^t  are  recorded  until  the  time  W 
is  reached.  At  this  point  the  time  span  W  is  then  initialized  at  the  next 
pulse  in  the  data  set,  which  may  be  either  positive  or  negative,  and  the  num- 
ber of  positive  or  negative  pulses  within  each  time  increment  ^t  are  recorded 
until  the  time  W  is  reached.   This  stepping- ahead  from  pulce-to-pulse  pro- 
cedure in  initializing  the  time  span  W  is  continued  until  the  last  data  pulse 
in  the  data  set  falls  within  the  last  time  increment  A t  at  the  time  W  from 
the  initialized  pulse. 
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Figure      81.      Modified      Two-  Channel       Recorder      Output 
Trace. 
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By  this  method,  four  sets  of  frequency  distributions  are  obtained.   The 
number  of  occurrences  of  the  radioactive  source  being  detected  at  the  upper 
(or  lower)  detector  level  per  unit  time  at  the  time  t,  initially  being  detec- 
ted at  the  lower  (or  upper)  detector  level  yield  the  frequency  distributions 
l(Z-,,Z  ,t)  and  l(Z  ,Z  ,t).   Also,  the  number  of  occurrences  of  the  radioactive 
source  being  detected  at  the  upper  (or  lower)  detector  level  per  unit  time  at 
the  time  t,  initially  being  detected  at  the  upper  (or  lower)  detector  level 
yield  the  frequency  distributions  l(Z2,Z  ,t)  and  l(Z-j_,Z,  ,t). 

Figures  82  through  85  represent  examples  of  the  particle  concentration 
frequency  distributions  obtained  from  the  experimental  data.   Computer  Pro- 
gram No.  k,    explained  in  Appendix  D,  was  used  to  obtain  the  experimental 
particle  concentrations. 

An  attempt  was  made  to  determine  the  diffusion  coefficient,  E  ,  for  each 
water  flow  rate  by  two  methods  of  applying  a  least  squares  analysis  to  the 
experimental  particle  concentration  data.   The  first  method  produced  dis- 
persion coefficients  for  each  water  flow  rate,  but  the  second  method  was 
unsuccessful.  An  attempt  to  explain  the  reason  for  lack  of  success  in  deter- 
mining dispersion  coefficients  by  the  second  method  will  be  given  in  Section 
5.^.  An  explanation  of  each  method  will  first  be  given. 

Method  1.   The  first  method  minimized  the  sum  of  the  squares  of  the  errors 
in  the  time  corresponding  to  the  maximum  experimental  I/I^  for  five  distances 
between  detector  levels  for  a  given  water  flow  rate. 

The  maximum  value  for  the  correlated  C/C^  occurs  at  the  time  correspond- 
ing to 

i{C/C-)-  =   0   .  (89) 


it 


iJU 


120 


4  6 

time  (seconds) 


°  1MM 


□ 


IMP 


2.5  in.Hg 

A  Z  =  2inches 


& 


8 


10 


Figure      82.      Particle      Concentration      Data,      at      Flow       Rate 


Bj,       Particles      Initially      at      Lower       Detector 
Level. 


140 


120 


100    _ 


4  6 

time     (seconds) 


3 


Figure      83.      Particle      Concentration       Data,       at       Flow      Rate 
Bj,      Particles'     Initially      at      Upper      Detector 
Level. 
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Figure     84.     Particle      Concentration       Data,       at  >    Flow      Rate'     B5, 
Particles      Initially      at  :    Lower      Detector      Level. 
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Figure     85. 


Particle       Concentration       Data,       at      Flow      Rate 
Particles       Initially       at       Upper       Detector      Level. 
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13* 


Therefore,  from  Equation  85,  one  has 


CO 


(C/Cj-1 


e  "  cos(ntfx)cos(nTfx   ) 


Tl=i 


1       e-(x-x0)2/4r 


(90) 


Differentiating  with  respect  to  time  gives 


(C/Cc*)-1 


L    2 


^t 


=  0=V~  -nV51  cos(n77x)cos(n7Tx0)e"n2^2^ 


L   e-'^o'2/4tt    L5(x-x0)2e-^-^4r,      (9l) 


1  l^/o  i       JA        ■*"/? 


Solving  for  the  time  corresponding  to  (C/C^,  )TriDV,    one  obtains 


max' 


t    = 


(aE^Uz)2 


r        ^  Uz)74Ezti 


81T    e 


_    L3(E^t)"3/2 


)      n  cos(ntrx)cos(n-&x0)  e  +1 


nr| 


where 


(92) 


4Z    =  (x-xn)L      . 


(93) 


For  a  given  flow  rate  of  water  and  five  distances  "between  detector  levels, 
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4Z,    one  1 


y~    >Tn     *bn  -    ^expLi'       _       0 


b%, 


bT, 


(94) 


where  (t        )     is  the  time  corresponds :.  i  to  (i/l,*,  ) 

v  exp'm  ^  '  °°   'max 

The  particle  concentration  frequency  distributions  obtained  from  con- 
sidering the  source  plane  at  the  lower  detector  level,  l(Z,  ,Zp,t),  for  all 
data  sets  were  normalized  and  combined  for  each  distance  between  detector 
levels  and  water  flow  rate,  and  are  shown  in  Figs.  86  through  90.   The 

experimental  values,  (t   )  were  obtained  from  these  figures.   Thus, 

exp  m  °  ' 


where 


arm 


[^-(texp)m]2=    2  &ai- (*expV]  ^     >  (95) 
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inaccuracy  was  the  fact  that  the  time  constants  of  the  ratemeters  were  greater 
than  one- third  of  a  second.   Observing  Figs.  57  through  71,  except  in  the 
cases  of  high  fluid  velocity  and  a  small  distance  between  detector  levels 
where  the  residence  times  were  very  short,  an  exponential  relationship  is 
seen  to  exist  between  particle  speed  and  the  distance  between  detector  levels. 
Because  of  errors  in  measurement  of  short  residence  times,  only  the  data  ob- 
tained when  the  distance  between  detector  levels  was  greater  than  one  inch 
were  used  for  an  extrapolation  to  zero  distance  between  detector  levels  to 
produce  the  results  shown  in  Figs.  72,  73>  and  7^- 

Figures  72,  73>  and  7^  clearly  illustrate  that  the  speed  of  the  particles 
within  the  fluidized  bed  is  greater  for  upward  movement  (Z-,  to  Z  )  than  for 
downward  movement  (Z  to  Z  ).   Since,  on  the  average,  the  same  number  of 
particles  move  upward  as  downward,  this  means  that  greater  than  fifty  percent 
of  the  cross- sectional  area  of  the  fluidized  bed  must  be  utilized  for  downward 
particle  movement. 

5.2  Particle  Traffic 

Figure  75  presents  the  particle  traffic  data  as  a  function  of  the  average 
fluid  velocity.   The  data  r(Z  )  were  obtained  at  one  vertical  height  for  all 
of  the  data,  but  the  data  r(Zp)  were  actually  obtained  at  five  different 
vertical  heights.  No  trend  dependent  upon  vertical  height  was  obtained,  so 
all  of  the  data  at  Z  were  combined  and  are  shown  in  Fig.  75  as  a  single  data 
value  at  each  water  flow  rate  studied.   It  was  concluded  that  r(Z)  had  no 
dependence  on  Z  over  the  range  of  vertical  heights  studied.   However,  the 
particle  traffic  was  measured  only  in  the  middle  sections  of  the  fluidized 
bed,  and  thus,  no  data  were  available  to  determine  the  effect  on  particle 
traffic  measurements  close  to  the  top  or  bottom  of  the  fluidized  bed. 


1^ 


5-3  Reflectivity  Data 

Figure  79  presents  the  experiments   reflectivity  data  "based  on  an  average 
m(Z.  )  for  each  data  set.  Also  shovr.  :,re  the  correlated  probability  of  re- 
flection curves  (l-e  '        )  determined  frorr  b!   data  "by  least  squares  analyses, 
limited  to  the  height  of  the  fluidized  bed  for  each  water  flow  rr.te.   Since 
no  data  were  obtained  for  A  Z  greater  than  six  inches,  no  comp.-ricon  is  avail- 
able between  experiment  and  theory  for  large  A  Z.   For  the  lowest  water  flow 
studied,  B-^,  data  were  obtained  over  a  larger  percentage  of  the  distance  from 
Z^  to  L  than  for  the  other  flow  rates  studied,  and  the  data  appear  to  support 
the  theory,  as  shown  in  Fig.  79* 

5.^  Particle  Concentrations 

Two  methods  of  determining  the  dispersion  coefficient,  E  ,  were  given 
in  Section  k.6,   "based  upon  a  diffusion  theory  model.   The  first  method  pro- 
duced results  reported  in  this  section.   The  second  method  gave  no  results. 
Before  attempting  to  determine  the  reason  for  this,  the  correlated  distribution 
(C/Coj)  was  calculated  for  the  five  water  flow  rates  studied  by  assuming  differ- 
ent values  of  the  dispersion  coefficient  and  comparing  the  resulting  distribu- 
tions for  each  dispersion  coefficient  with  the  experimental  data  (I/Iq,)  deter- 
mined from  Figs.  86  through  90.   The  effect  on  the  correlated  distributions 
as  a  result  of  varying  E^  for  each  water  flow  rate  was  as  follows.  As  Ez 
was  increased,  the  maximum  value  of  C/C^  was  decreased  for  each  distance 
between  detector  levels,  and  the  maximum  values  occurred  at  lesser  time  incre- 
ments.  The  dispersion  coefficients  that  produced  the  correlated  distributions 
that  came  closest  to  representing  the  experimental  distributions  for  each 
water  flow  rate  were  then  tabulated. 
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It  was  noted  that  for  time  increments  less  than  about  0.3  second,  the 
correlated  (C/Cgc)   distribution  did  not  compare  very  closely  with  the  experi- 
mental (i/loo)  distribution,  even  though  the  distribution  for  time  increments 
greater  than  0.3  second  were  closely  matched.  This  would  indicate  that  the 
data  for  short  time  increments  (less  than  0.3  second)  were  incompatible  with 
the  rest  of  the  data  in  each  experiment,  and,  in  attempting  to  determine  a 
dispersion  coefficient  by  Method  2  using  all  of  the  data  in  each  experiment, 
results,  if  any  could  be  obtained,  would  be  highly  questionable. 

Figure  9±   presents  the  experimental  (i/l^)  data  for  the  lowest  water 
flow  rate  studied,  with  the  correlated  (c/Cco)  distribution  for  a  dispersion 
coefficient  of  1  in  /sec  superimposed.   The  data  below  about  0.3  second  are 
not  compatible  with  the  correlated  distribution.  Figure  92  presents  the 
correlated  (C/C^,  )  distribution  for  a  water  flow  rate  of  18.8  gal/min., 
using  a  dispersion  coefficient  of  10  in2/sec.   Comparison  of  this  figure  with 
Figure  38  illustrates  the  marked  difference  below  a  0.3  second  time  increment, 
although  the  dispersion  coefficient  used  gave  the  best  representation  at 
greater  time  increments. 

Table  6  gives  the  dispersion  coefficients  determined  as  outlined  above. 
Table  6  also  gives  values  of  dispersion  coefficients  determined  from  a  micro- 
scopic viewpoint,  as  discussed  below. 

In  a  medium  where  only  scattering  reactions  occur,  expressions  for  par- 
ticle current  density  have  been  developed  (3)-  Assuming  only  Z-dependence 
and  that  the  particles  move  with  average  speed  v,  the  following  current 
densities  result. 


t  =  c5  .   _i_|J£l     ,  (101) 

+     4       6Z  I  iZ  JO    '  ' 
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Table  6.   Dispersion  Coefficient  Ez  as  Determined 
by  Single  and  Multi-Particle  Experiments. 


Dispersion  Coefficient 
(in^/sec) 


Water  Flow  Rate  (gal/min) 


13.5     17.0     18.8     22.3     24.0 
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Method  1 
Method  2 


13.6     11.9     24.0     23.6     23.6 


Estimated  from 

Correlated 

Distribution 


1.0 


8.0     10.0     11.0     12.0 


Determined  from 
Single  Particle 
Experiments 


6.7    25.3    38.3    ^1.6    41.7 
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and 


4     eZ      z 


o 


(102) 


where  2f  is  the  probability  of  interaction.   Applying  Equations  101  and  102 
to  fluidized  beds,   ZT  is  equivalent  to  the  reflection  coefficient  (B),    for 
Q     is  an  interaction  probability,  as  implied  in  Equation  26,  and  is  deter- 
mined from  single  particle  experiments. 

In  experiments  with  more  than  one  particle,  determination  of  J        and 
J„   cannot  be  made;  only  the  net  movement  of  particles  through  the  area 
studied  can  be  made.  The  net  movement  is  given  by 


Jnet  ~"   J+  "  J- 


(103) 


Substituting  Equations  101  and  102  into  Equation  103  give 


net 


3/?  j_^Z_ 


(104) 


In  diffusion  theory,   Equation  10 4  is  Fick's  law  of  diffusion, 


net  z  ** 


(105) 


where  E  is  the  diffusion  coefficient.  From  Equations  10 4  and  105, 


JS_  = 


(106) 
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Thus,  we  have  a  means  of  comparing  results  based  on  single  particle  experi- 
ments with  multi-particle  experime;  ui  results.   The  q~  entities  v  and  R   were 
determined  from  microscopic  invest    ions  (Sections  K.:   and  ^-.5),  and  Ez 
has  been  determined  from  macrosc     investigations. 
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6.0   CONCLUSIONS 

The  results  of  this  investigation  have  led  to  the  following  conclusions 
about  fluid-solids  fluidized  beds. 

Frequency  distributions  of  resic     times  agree  very  veil  with  gamma 
distribution  correlations  based  on  the  parameters  p  and  1)   determined  from 
the  experimental  data.   The  gamma  distribution  parameter  i)   increases  with 
increasing  distance  between  detector  levels,  while  the  parameter  p  tends  to 
have  a  minimum  value  at  about  three  inches  distance  between  detector  levels 
for  all  water  flow  rates  studied,  as  illustrated  in  Figs.  ^7  through  56. 

The  logarithm  of  the  vertical  speed  of  the  particles  tends  to  increase 
linearly  with  decreasing  distance  between  detector  levels  for  all  water  flow 
rates  studied,  as  illustrated  in  Figs.  57  through  Jl. 

The  vertical  speed  of  the  particles  increases  with  increasing  fluid 
velocity,  as  illustrated  in  Figs.  72,  73*  and  Jk. 

More  than  fifty  percent  of  the  cross- sectional  area  of  the  fluidized 
bed  is  used  for  downward  particle  movement. 

About  fifty  residence  times  per  data  set  are  sufficient  to  study  particle 
speeds. 

In  the  middle  sections  of  the  fluidized  bed,  the  particle  traffic  in- 
creases with  increased  fluid  velocity,  as  illustrated  in  Fig.  75 • 

A  frequency  distribution  of  the  probability  of  reflection  for  downward 
particle  movement  is  the  same  as  for  upward  particle  movement,  as  illustrated 
in  Figs.  76,  77,  and  78.   The  probability  of  reflection  increases  with  in- 
creased distance  between  detector  levels  for  all  water  flow  rates  studied, 
as  illustrated  in  Fig.  79. 

The  analytical  model  developed  to  determine  the  concentration  of  particles 
with  respect  to  time  and  distance  from  the  source  plane  gives  distributions 
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that  are  in  good  agreement  with  the  experimental  f:  ?ouency  distributions  of 
particle  concentrations  for  distances  between  det   .or  levels  greater  than 
one  inch,  as  illustrated  in  Fig.  91* 

The  determination  of  particle  (    jrsion  properties  by  microscopic  and 
macroscopic  investigations  has  given  particle  dispersion  coefficients  with 
differences  of  less  than  an  order  of  magnitude  for  all  water  flow  rates 
studied. 

The  study  cf  individual  particle  mobility  by  the  method  developed  in  this 
thesis,  utilizing  the  apparatus  used  in  this  thesis  to  obtain  the  experimental 
data,  led  to  the  following  conclusions. 

Residence  times  obtained  from  the  experimental  apparatus  which  were  less 
than  0.2  second  were  not  reliable. 

Measurements  at  the  top  and  bottom  of  the  fiuidized  bed  should  be  made 
to  determine  how  end  effects  affect  the  particle  traffic. 

For  low  water  flow  rates,  a  data  correlation  was  obtained  to  enable  a 
prediction  of  particle  reflection  characteristics. 

Experimental  measurements  with  greater  distances  between  detector  levels 
and  higher  water  flow  rates  are  needed  to  produce  more  particle  reflection 
data  to  compare  with  the  analytical  expression  for  particle  reflectivity  de- 
veloped in  this  thesis. 

Analyses  of  the  residence  time,  reflectivity,  and  particle  concentration 
data  led  to  the  conclusion  that  the  data  obtained  for  distances  between  de- 
tector levels  of  less  than  one  inch  were  not  representative  of  the  rest  of  the 
data,  and  this  was  attributed  to  the  limitations  of  the  ratemeters  used  in 
the  experimental  study.  Better  instrumentation  sensitivity  and  recovery  times 
are  needed  to  obtain  data  at  high  water  flow  rates  and  small  distances  be- 
tween detector  levels. 
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Additional  experiments  with  different  fixed  bed  heights  and  different 
particle  sizes  and  densities  using  t  radioactive  tracer  technique  should 
prove  worthwhile. 

A  large  amount  of  labor  was  exp  ^ed  in  reducing  an  correlating  the 
raw  data  from  the  dual-channel  re^    r  trace  and  digital  printout,  and  trans- 
ferring the  data  to  IBM  cards.   It  is  suggested  that  should  further  experi- 
ments occur,  a  method  should  be  developed  that  will  enable  more  data  collec- 
tion at  reduced  labor  by  automatically  recording  the  raw  data  in  a  form  suit- 
able for  digital  computation. 

The  preliminary  data  obtained  in  this  thesis  indicate  that  the  use  of 
the  methods  of  data  collection  and  analysis  developed  and  presented,  is  fea- 
sible in  assisting  in  the  determination  of  the  behavior  of  individual  particles 
in  fluidized  beds. 


c£ 
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APPENDIX  A 


Description  and  Explanation  of  the  IBM 
1620  Computer  FORTRAN  Program  No.  1  for 
Correcting  the  Time  Increment  Data 


This  computer  program  was  written  to  correct  the  experimental  time 
increments  obtained  when  the  tagged  particle  passed  from  the  upper- to-lower 
and  lower-to-upper  detector  levels.   The  corrections  were  applied  on  the 
basis  of  Table  1.   Table  7  defines  the  various  symbols  utilized  in  this 
computer  code. 


Table  7*   Input  data  and  variables  required  for 
the  IM  1620  Computer  Program  No.  1. 


SYMBOL  EXPLANATION 


JMAX  Number  of  time  increment  data  for  given  set 

D(j)  •     Individual  time  increment  data  in  given  set, 

J— 1  to  JMAX,  In  pulses  per  channel  x  10"^ 
p(j)  Individual  time  increment  data  converted  to 

time  in  seconds 
I  Time  increment  data  set  in  the  process  of 

correction 


First  the  program  reads  the  experimental  data  D(j)  and  the  number 
of  individual  data  values  JMAX  in  the  Set  I.   The  program  then  converts 
the  data  to  time  in  seconds  P(J). 

The  program  determines  which  data  represents  upper-to-lower  and  lower- 
to-upper  detector  level  time  increments.   By  applying  the  equation  necessary 
from  Table  1,  each  time  increment  for  both  the  upper-to-lower  and  lower-to- 
upper  detector  level  is  then  corrected  to  true  time. 

After  each  time  increment  in  the  set  has  been  so  treated,  the  program 
punches  IBM  cards  containing  the  data  set  I,  the  number  of  data  values  in 


the  set  JMAX,  and  the  individual  values  of  corrected  time  incre- 
ments . 

Input  to  Computer  Program  No.  1 

The  data  obtained  from  the  fluidized  bed  apparatus  was  transferred 
to  IBM  cards  and  served  as  the  input  to  Computer  Program  No.  1.  Eight 
individual  time  increments  were  typed  per  card.  A  plus  or  minus  sign 
preceded  each  time  increment  to  designate  detection  of  the  tagged  par- 
ticle at  either  the  upper  or  lower  detector  level  respectively. 

Seventy- six  sets  of  data,  with  water  flow  rates  and  distances  be- 
tween detector  levels  for  each  data  set,  are  given  in  Table  k.      Table  8 
gives  the  input  data  to  Computer  Program  No.  1.  Each  data  set  is  indi- 
cated by  Data  Set  I.   The  first  number  shown  in  each  data  set  is  the 
number  of  time  increments  in  the  set  (JMAX).   Each  data  card  is  read  into 
the  computer  as  FORMAT  (l4,  Ik,   8E7-3)-   The  first  four  columns  give  the 
card  number.   The  next  four  columns  give  the  subscript  J  of  the  first 
time  increment  D(j)  on  each  data  card.  J  —  1,  2,    ...,  JMAX.   The  indi- 
vidual time  increments  are  recorded  as  the  number  of  pulses  from  the  Time 
Interval  Meter  x  10~3  stored  in  each  channel  of  the  multi-channel  analyzer 
before  rapid  change  in  the  counting  rate  caused  transfer  of  storage  to  the 
next  channel. 

Output  of  Computer  Program  No .  1 

The  output  of  Computer  Program  No.  1  is  given  in  Table  9-   The  output 
was  punched  as  FORMAT  (ik,    Ik,   8F8.3).  Each  data  set  is  given  as  Data 
Set  I.   The  first  number  shown  in  each  data  set  is  the  number  of  time 
increments  in  the  set  (JMAX).   The  first  four  columns  per  card  indicate 
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Table   8.  Input   to   Computer  Program  17o.   1 

Data  Set   6 


232 

10001 

+  00.001 

+00491 

+01371 

+00356 

-00209 

+00077 

-00  197 

-00  0  77 

20009 

-00168 

-00170 

-00300 

-00113 

-00648 

-00576 

-07350 

-00139 

30017 

-06000 

-00312 

+00280 

-0  0  961 

-00700 

-CC100 

-00140 

_  r-  <*>  k,  c  II 

40025 

-C0656 

-00431 

-00200 

-00233 

+  03.330 

+04332 

+00812 

+00  5  86 

50033 

+00278 

-00212 

-00358 

-00579 

-05909 

-00934 

-03928 

-0  0  507 

60041 

-00253 

-0  0  2  2  0 

-00409 

-00162 

-00185 

-00407 

-00121 

-0  0  79  0 

70049 

-00134 

-00300 

-008/  1 

+0  0150 

+00363 

+00  5  59 

+00606 

+00  263 

80057 

+00333 

+00672 

+00124 

+0  0061 

+00299 

+00267 

+00152 

+0  0  2  24 

90065 

+00147 

+00585 

+00389 

-0  0  207 

-00400 

-00375 

-00300 

-0  37  5  8 

100073 

+00300 

-00400 

-00579 

-00192 

-0  0  40  0 

-CO  200 

- 0 0 2  00 

-Or;  42  2 

110081 

-00259 

-00747 

+00200 

+003C0 

+00233 

+00141 

-00426 

-00113 

120089 

-00190 

-00279 

-00067 

-16147 

+00300 

+99999 

+00156 

+00211 

130097 

+00486 

+00593 

-00158 

-00400 

-00300 

-00172 

-03370 

-24476 

140105 

+00160 

+00700 

+00152 

+00299 

+00  508 

-00322 

+00603 

+00636 

150113 

+  00625 

+00571 

+00724 

+00395 

+0ol84 

-00167 

+C0195 

+C4331 

16  0121 

+01152 

+01204 

-00174 

-00676 

-00519 

—  00  348 

-02068 

-0  0  5  52 

170129 

-00132 

-00215 

-00889 

+0  0r29 

+  0vJ  561 

+00414 

-00440 

-14596 

180137 

-00077 

-00274 

-00107 

-00135 

-00003 

-00136 

-21461 

-0014  7 

190145 

-00398 

-00300 

-00490 

+0  02  0  0 

-0  0  600 

-00240 

-00347 

-0  3126 

200153 

+C0900 

+00C79 

+00835 

-00311 

-0  0  938 

+00937 

+00168 

+0014  8 

210161 

-07133 

-00334 

-00332 

-0097  1 

-04928 

-00425 

-00195 

+00376 

220169 

+00232 

+18067 

+00299 

+00120 

+00568 

-01010 

-00097 

+00427 

220177 

+00269 

+00329 

+00136 

+00278 

+00885 

+00452 

+07963 

+0  1218 

240185 

+00524 

-00497 

-00421 

-04827 

-00141 

-00238 

-00244 

+01064 

2  5  019  3 

+04504 

+00514 

-00572 

-00684 

-01793 

-00313 

-00021 

O  1  ct\J  O 

260201 

-00752 

-00240 

-00109 

-00320 

-00386 

+00080 

+00125 

+10097 

270209 

+01200 

+00560 

+00444 

+00490 

+00597 

+01183 

+15436 

+00310 

280217 

+01184 

+02718 

+00989 

+00C99 

-00644 

-00171 

+0C73C 

-00  334 

290225 

-00069 

-0C252 

-00188 

-02753 

+00128 

+0C559 

+00187 

+C03C9 
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Table  8  (cont.) 

Data  Set  7 


232 

10001  -00001  -00104  -00472  -00595  -00248  -00338  -00357  -00053 

20009  -00522  -01229  -00572  -00499  -01043  -00260  -00599  -00122 

30017  -20055  -00440  -00469  -00277  -00131  -00431  -00389  -00691 

40025  -00074  -01183  -00616  -00601  -00761  -00148  -01074  +00696 

50033  +00220  +00372  +00232  +00687  -00355  -00354  -00085  +00423 

60041  +00344  +00586  -00300  -00100  -01394  -01103  -01778  -04202 

70049  -00871  -r-0370  -00359  -08466  +00439  +97715  +00487  -00669 

80C57  +00531  +00352  +0CC91  +00717  +00151  +00140  +00474  +01705 

90065  +00700  +00968  +00200  -02027  -00600  -00200  -0H090  -00090 

100073  -00086  -00642  -00183  +00346  -00595  -14012  +00072  +00524 

110081  +01200  +00529  +28747  +00659  +01399  -00751  +00373  +00247 

120089  +CT244  +00188  +00712  +00123  +00367  +01336  +00311  +00245 

130097  +00481  -00645  +00392  -003^0  -00188  -00146  -0020O  -38400 

140105  -00600  -00300  -01293  -01507  -00120  -00807  -02559  -00392 

150113  +01499  +00830  +00157  +00860  +0J247  +00253  +06908  +00546 

160121  +01124  +00292  +00315  +00191  +00526  +00786  +00826  -00337 

170129  -00270  -00210  -00392  -00246  -00213  -00368  -00090  +00448 

180137  -00141  +00281  +00655  -  '0831  +00839  +00285  +00097  +00457 

190145  +00545  +00205  +16883  +00220  +02912  +00331  +00502  +00068 

200153  -00962  -00214  -08572  +00387  +00644  +00243  +00454  +00914 

210161  +00869  +00451  -00184  -0-  274  -0^519  -01339  -04105  -00424 

220169  -00544  -00098  +00443  +00313  +00144  +00114  +01232  +00447 

230177  +00689  +00315  +00124  +00^72  +00251  +00455  -0H267  +00576 

240185  -07247  +00560  +01171  +00520  +00177  +C0376  +00434  +00145 

250193  +00094  +00526  +00067  +00535  +C0684  +00470  +00513  +00673 

260201  +00172  -00072  -00285  -00200  -01060  -00093  -00422  -00547 

270209  -00118  -00931  +00873  -00114  -00315  -00315  -00077  -01381 

280217  +00065  -00494  -00662  -00381  -00146  -0C351  -00813  -00056 

290225  -00692  -0C074  -06216  -00368  -07730  -00244  -00168  -17801 


Table  9.   Output  of  Computer  Program  ITo.  1 


Data  Set  6 
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232 

► 

6 

1 

.001 

.736 

2.056 

.534 

-.267 

.199 

-.251 

-.115 

6 

9 

-.252 

-.255 

-.450 

-.169 

-.972 

-.864 

-11.025 

-.20  8 

6 

17 

-9.000 

-.468 

.585 

-1.22d 

-1.050 

-.150 

-.210 

-.891 

6 

25 

-.984 

-.646 

-.300 

-.349 

5.876 

6.49  3 

1.218 

.  3  7  r-< 

6 

33 

.417 

-.271 

-.537 

-.863 

-8.863 

-1.401 

-5.892 

-.760 

6 

41 

-.379 

-.330 

-.613 

-.243 

-.277 

-.610 

-.181 

- 1 .18  5 

6 

49 

-.20  1 

-.45  0 

-1.261 

.  340 

.544 

.8  38 

.90^ 

.3Q4 

6 

57 

.499 

1.008 

.186 

.091 

.448 

.400 

.228 

.3  36 

6 

65 

.2  2U 

.877 

.583 

-.264 

-.600 

-.562 

-.450 

-13.182 

6 

73 

.621 

-.511 

-.868 

-.283 

-.60  0 

-.300 

-.300 

-.62? 

6 

81 

-.388 

-1.120 

.43  7 

.450 

.349 

.211 

-.544 

— .  169 

6 

89 

-.285 

-.418 

-.100 

-24.22<, 

.62  1 

149.998 

.2  34 

.  -/  M   \J 

6 

97 

.729 

.889 

-.201 

-.60o ■ 

-.450 

-.258 

-5.05  5 

-36.714 

6 

105 

.360 

1.050 

.228 

.448 

.762 

-.All 

1.147 

.9  54 

6 

113 

.937 

.856 

1.086 

.592 

.276 

-.213 

.427 

6 

121 

1.728 

1.806 

-.222 

-1.014 

-.77  8 

-.522 

-3. 102 

-.828 

6 

129 

-.198 

-.322 

-1.33  3 

1.016 

.841 

.621 

-.562 

-21.8C4 

£, 

137 

-.115 

-.411 

-.1  60 

-.202 

-.004 

-.2  04 

-32.19  1 

-.22^ 

6 

14  5 

-  .  5  °  7 

-.450 

-.735 

.437 

-.767 

-.360 

-.  520 

-4.6^9 

6 

.15  3 

1.611 

.118 

1.25  2 

-.39  7 

-1  .407 

1.669 

.252 

.222 

6 

16] 

-9.182 

-.501 

-.496 

-1.456 

-7.392 

-.637 

-.292 

.7  56 

6 

169 

.34  8 

27.  100 

.44  3 

.180 

.652 

-1.291 

-.  145 

•  c  -+■  *-* 

6 

177 

.40  3 

.493 

.204 

.417 

1.327 

.678 

11.944 

1.327 

6 

185 

.786 

-.635 

-.631 

-7.240 

--.211 

-.432 

-.366 

1.877 

6 

193 

6.75  6 

.771 

-.731 

-1.026 

-2.689 

-.469 

-.031 

-11.109 

6 

201 

-1 .  128 

-.360 

-.163 

-.480  . 

-.579 

.197 

.137 

15. 145 

6 

2o9 

1.80  0 

.840 

.666 

.735 

.895 

1.782 

24.729 

•  465 

6 

217 

1.776 

4.077 

1.483 

.148 

-.823 

-.256 

1.344 

-.490 

6 

225 

-.10  3 

-.378 

-.282 

-4.129 

.296 

.838 

.280 

.463 

Taole   9    (cont.) 


Data  Set  7 


iCh 


232 

> 

7' 

1 

-.001 

-.156 

-.708 

-.892 

-.372 

-.507 

-.535 

-.0  79 

7 

9 

-.783 

-1.843 

-.858 

-.746 

-1.564 

-.390 

-.  398 

-.183 

7 

17 

-30.082 

-.660 

-.703 

-.415 

-.196 

-  .  646 

-.56  3 

_  i  r  '  a 
X   .  ■  J  ->  o 

7 

25 

-.111 

-1.774 

-.924 

-.90  1 

-1.141 

-.222 

-1.611 

1.269 

7 

33 

.330 

.558 

.348 

1.030 

-.453 

-.531 

-.127 

.837 

7 

41 

.516 

.379 

-.383 

-.]  5U 

-2.091 

-1.654 

-2.667 

—  (■     '■>  n  "> 

7 

49 

-1.30  6 

-.555 

-.538 

-12.699 

.86  4 

146.572 

.73  0 

-.355 

7 

57 

1*020 

.523 

.136 

1.07  5 

.226 

.210 

.711 

2.557 

7 

65 

1.050 

1.452 

.300 

-2.591 

-.900 

-.300 

-.  135 

-.  133 

7 

73 

-.129 

-.963 

-.274 

-  .703 

-.760 

-21.018 

.195 

.766 

7 

81 

1.800 

.793 

43.  120 

.96b 

2.09  8 

-.96  0 

.75  0 

.376 

7 

69 

•  366 

.282 

1.068 

.  184 

.350 

2.004 

.466 

.3  67 

7 

97 

.721 

-.624 

.7S3 

-.363 

-.282 

-.219 

-.300 

-57.600 

7 

105 

-.900 

-.450 

-1.939 

-2 • 260 

-.16  0 

-1.210 

-3.8  36 

-  •"<  ^ 

7 

113 

2.645 

1.245 

.235 

1.290 

.370 

.3  79 

10.362 

.819 

7 

121 

1.636 

.433 

.472 

.266 

.789 

1.179 

1.2  39 

-.43  0 

7 

129 

-.405 

-.315 

-.588 

-.369 

-.319 

-.552 

-.  135 

O  ft  O 

.  O    w  *_. 

7 

137 

-.  180 

.537 

.932 

-1.062 

1.516 

.427 

.145 

.5-5 

7 

145 

.317 

.3  07 

25.324 

.330 

4.368 

.496 

.75  3 

.10  2 

7 

153 

-1.229 

-.321 

-12.858 

.775 

.966 

.364 

.681 

1.371 

7 

161 

1.303 

.676 

-.235 

-.411 

-.778 

-2.008 

-6.157 

- .  636 

7 

169 

-.816 

-.147 

.871 

.469 

.216 

.171 

1.845 

.  6  7  0 

7 

177 

1.033 

.472 

.186 

.553 

.376 

.682 

-.341 

1.094 

7 

135 

-9.264 

1.068 

1.756 

.780 

.265 

.564 

.631 

~>   1  "7 

.   4^.   A.    . 

7 

193 

.141 

.789 

.100 

.802 

1.026 

.705 

.769 

1.009 

7 

201 

.258 

-.122 

-.427 

-.300 

-1.590 

-.139 

-.633 

-.620 

7 

209 

-.177 

-1.396 

1.569 

-.145 

-.472 

-.472 

-.115 

-2.071 

7 

217 

.189 

-.631 

-.993 

-.571 

-.219 

-.526 

-1.219 

7 

225 

-1.338 

-.111 

-9.324 

-.552 

-11.595 

-.366 

-.252 

-2  6.701 
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the  data  set.   The  next  four  columns  indicate  the  subscript  J  of  the 
first  time  increment  D(j)  on  each  data  card.   The  individual  time 
increments  are  given  in  seconds,  with  eight  time  increments  per  card. 
A  plus  or  minus  sign  preceded  each  time  increment  to  designate  detection 
of  the  tagged  particle  at  either  the  upper  or  lover  detector  level 
respectively. 


Lorcio   lHnp;ram  for  Computer  Program  l\o,    1 


LOO 


P(J)  =  -P(J)/a 


ICY 


33 


Calculate 

R=    Exp{- 0.21715  + 

0.13268    log   P(J) 

n 


P(J)  = 
P(J)/R 


T 


44 


Calculate 

R=  Ex?(l.006273 


0.7186 


lognP(J) 


22 
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FORTRAN  Source  Program  for  Computer  Program  ITo.  1 


PROGRAM  FOR  CORRECTING  ALL  DATA 

DIMENSION  D(300) »P( 300) 
101  FORMAT  (  1 4  ♦  1 4  t  E7 . 3  >  E7 . 3  » E7 . 3  »  E~ . 3  >  E7 . 3 ♦ E7 . 3  >  E7 . 3  >  E7 . 3 ) 
103  FORMAT ( 15) 
10  5  FORMAT ( I  4 , 14 , FQ • 3 » FS . 3 vF8 . 3 > F8 . 3 » Fu .3 > F8.  3 > F8 . 3 »FS  .  3 ) 

READ  103»!< 

1=1 

200  READ  103»JMAX 
PUNCH  103»JMAX 

201  J=l 

202  READ  1C1,M,N,D(J) »D( J+l) »D(J+2) »D( J+3) »D( J+4) >D( J+5) »D( J+6) »D(J-r7) 
J  =  J+8 

IF(JMAX-J) 1»1»202 

1  J  =  l 
P(J)=1.5*D(J) 
IF(P(J) )3,2,2 

2  J=J+1 
IF(JMAX-J)203»4»4 

4  P(J)=1.5*D<J> 
IF(P(J) )5>2>2 

5  P(J)=-P(J) 
IF(P(J)-0.14)7>6»6 

6  P( J)=-P(J)/1.1734 
GO  TO  3 

7  X=LOGF(P( J) ) 

Y = 1.3 3479+ (0. 65478 )*X 

R=EXPF(Y) 

P< J)=-P< J)/R 

GO  TO  3 

3  J=J+1 

IF( JMAX-J)203»8»8 

8  P( J)=1.5*D( J) 
IF(P( J) )3»9,9 

9  IF(P( J)-1.5) lltl0»10 

10  P< J)=P(J)/0.85 
GO  TO  2 

11  IF(P(J)-0.12)12»13»13 


ifc 


13  X=LGGF(P( J)  ) 

Y =-0.2 17 15+10. 1326  3  )*X 

R=EXPF(Y) 

P( J)=P( J) /R 

GO  TO  2 
12  X  =  LOGF(P( J)  ) 

Y=l . 006273+ ( 0.7186 )*X 

R=£XPF(Y) 

P(J)=P(J)/R 

GO  TO  2 

203  J  =  l 

204  PU,\CM105,K,J»P(  J)  »P(  J  +  l )  ,P(J+2)  »P< J+3)  >P( J+4)  >P(  J  +  5  )  >P( J+6)  »P(  J  +  7) 
J  =  J  +  8 

IF (JMAX-J) 300 ♦ 204*204 
300  K  =  K  +  1 

GO  TO  200 
END 
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APPENDIX  B 


Description  and  Explanation  of  the  IBM 
1620  Computer  FORTRAN  Program  No.  2  for 
Computing  the  Average  Number  of  Times  a 
Particle  Passed  a  Detector  Level  per 
Unit  Time 


Tnis  computer  code  was  written  to  determine  the  average  number  of 
times  per  unit  time  a  tagged  particle  in  the  fluidized  bed  passed  a  given 
horizontal  plane.   This  calculation  was  done  by  counting  the  number  of 
times  the  tagged  particle  passed  each  of  the  two  detector  levels  in  a 
given  set  of  data  and  dividing  each  number  by  the  total  time  required  to 
obtain  the  data.  Figure  75  illustrates  the  results  obtained  from  this  pro- 
gram.  Table  10  defines  the  various  symbols  utilized  in  this  computer  code. 


Table  10.   Input  data  and  variables  required  for 
the  IBM  lo20  Computer  Program  No.  2. 

SYMBOL  EXPLANATION 

JMAX  Number  of  time  increment  data  for  given  set 

D(j)  Individual  time  increment  data  in  given  set, 

J  =1  to  JMAX,  in  seconds 

M  Data  set 

SUMT  Total  time  required  to  obtain  data  set 

ZT  Number  of  times  particle  passed  upper  level 

ZB  Number  of  times  particle  passed  lower  level 

RT  Average  number  of  times  particle  passed  upper  level 

RB  Average  number  of  times  particle  passed  lower  level 


Input  to  Computer  Program  No.  2 

The  input  to  Computer  Program  No.  2  is  the  output  of  Computer  Pro- 
gram No.  1.   The  input  is  read  into  the  computer  as  FORMAT  (ik,   Ik,   8F8.3) 
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Output  of  Computer  Program  ITo.  2 

Table  11  gives  the  output  of  Computer  Program  No.  2.  The  left-hand 
column  designates  the  data  set  number.  KB  is  the  average  number  of  times 
the  tagged  particle  passed  the  lover  detector  level  per  second,  and  RT  is  the 
average  number  of  times  the  tagged  particle  passed  the  upper  detector  level 
per  second.   SUMT  is  the  sum  in  seconds  of  the  absolute  values  of  the  indi- 
vidual time  increments  for  the  given  data  set  and  is,  therefore,  the  time 
required  to  obtain  the  data  set. 


Table  11.    Output  of  Conp.  .cr  Program  No.  2. 
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SET  NO 

RS 

1 

.16672560 

2 

.09450548 

3 

.11695540 

4 

.13426900 

5 

♦05934301 

6 

.21515310 

7 

.20535470 

8 

.08815150 

9 

.08581589 

10 

.09155243 

11 

.08233181 

12 

.14307640 

13 

.17278970 

14 

.11348069 

15 

.12455340 

16 

.14613030 

17 

.22857180 

18 

.18476280 

19 

.24350410 

20 

.11517251 

21 

.11542609 

22 

.16808850 

23 

.15844310 

24 

.18090110 

25 

.11236883 

26 

.20963420 

27 

.19083230 

28 

.11988491 

29 

.17872270 

30 

.14200264 

31 

.14953520 

32 

.23925280 

33 

.14225290 

34 

.16932390 

35 

.11938770 

36 

.15294970 

37 

.16924440 

38 

.13034090 

RT 

.15737650 
.09093924 
.13592120 
.15090410 
.10981867 
.17178390 
.18515580 
.08958486 
.09765256 
.08176820 
.10328900 
.09745790 
.08223032 
.17369490 
.14173320 
.13569240 
.21664630 
.16501710 
.19590180 
.13469327 
.18275790 
.21087470 
.17740220 
.18374990 
.13349801 
.22038460 
.18121050 
.11403686 
.14808460 
.08520158 
.12059290 
.21947990 
.15554760 
.15485180 
.16259470 
.17008000 
.20117730 
.13491430 


SUMT 

641.77300000 

1121.62800000 

949.07900000 

841.59400000 

1466.05300000 

599.57300000 

594.09400000 

1395.32500000 

1351.73100000 

1430.87400000 

1336.05700000 

964.51900000 

972.28000000 

863.53300000 

931.32700000 

766.43900000 

503.12400000 

709.01700000 

546.19200000 

1024.55000000 

831.70100000 

654.41700000 

738.43500000 

702.04100000 

1041.21400000 

558.11500000 

623.58400000 

1025.98400000 

733.33600000 

1091.52900000 

829.23600000 

505.74100000 

752.18100000 

690.98300000 

379.48700000 

817.26200000 

626.31300000 

874.62900000 


Table  11  (cont.) 


17: 


SFT  NO 


RB 


RT 


SUNT 


39 

.18596680 

.17155080 

'♦0 

.16316460 

.18204320 

'♦1 

.22216290 

•  18700040 

42 

.16976470 

.14678900 

43 

.17520630 

.18650990 

4  4 

.15603420 

.18248070 

45 

.15801098 

.26793160 

46 

.14800190 

.18873630 

47 

.16941140 

.16941140 

48 

.12674220 

.13371890 

49 

.22402360 

.20124160 

50 

.18539720 

.13161786 

51 

.21051470 

.20383170 

52 

.13754290 

.16810310 

53 

.20392410 

.16000200 

54 

.22753960 

.20653590 

55 

.20490950 

.23590420 

56 

.13238130 

.22510120 

57 

.17511360 

.15210600 

58 

.17245750 

.13303860 

59 

.17930070 

.13517940 

60 

.17708730 

.20710210 

61 

.21499230 

.20514530 

62 

.21355070  ' 

.21355070 

63 

.17144320 

.18948990 

64 

.20055220 

.24352770 

65 

.19473420 

.21453770 

66 

.17229080 

.13959330 

67 

.14466620 

.12314560 

68 

.21476790 

.18792190 

69 

.18971670 

.18647360 

70 

.14801120 

.18397660 

71 

.19706826 

.20104940 

72 

.16156630 

.17725230 

73 

.16260130 

.16830670 

74 

.17690540 

.17090860 

75 

.15796580 

.17403010 

76 

.18569990 

.16855840 

693.67200000 
741.58200000 
625.66700003 
733.43700000 
707.73700000 
756.24400000 
582.23800000 
736.47700000 
637.50100000 
860.01300000 
526.73000000 
706.59100000 
598.53300000 
735.20900000 
637.49200000 
571.32900000 
530.74400000 
608.61500000 
782.34900000 
811.79400000 
680.42100000 
666.33300000 
609.32400000 
599.33900000 
664.94300000 
558.45300000 
605.95400000 
795.16700000 
836.40800000 
595.99200000 
616.70900000 
722.91800000 
502.36400000 
637.50900000 
701.10100000 
667.02300000 
746.99700000 
700.05400000 


Lo^ic  Diagram  for  Computer  Program  No.  2 
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Str:rt; 

Road 

Data  (J),  J  MAX 


END 


5c  t 

J  =  1 

SUMT  =  0.0 

X  =  0.0 


Calculate 

RT  =  ZT/SUMT 

R3= Z3/SUMT 


1 


Punch 

M,  RB,  RT, 

SUMT 


J=  J  +1 


H 


ZT  =  0.0 
Z  B  =  0.0 


ZB  =  I.O  +  ZB 


X=  -X 


SUMT  = 
X  4  SUMT 


X  = 

Data  (J) 


ZT  =  1.0+ZT 
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FORTRAN  Source  Program  for  Computer  Program  ITo.  2 


C      PROGRAM  FOR  DETERMINING  THE  NUM3ER  OF  TIMES  A  PARTICLE  PASSES 
C      THE  TOP  AND  BOTTOM  LEVEL   ANS  IN  SECONDS 
DIMENSION  D(300) 

101  FORMAT ( U»I4»F8.3»F8.3»F8-3»F8.3»F8.3>F8*3»F8«3»F8.3) 

102  FORMAT ( 15 ) 

10  3  FORMAT ( I4,F15.8>F15.8>F15.8»F15.8) 
1=1 

200  READ  102»JMAX 
J  =  l 

201  READ  lOl»MtN»D(J)tD(J+l)  »D(  J+2)  »D( J+3)  »D{  J+4)  »D(  J  +  5  )  »D(  J  +  6)  »D(J-*-7) 
J  =  J  +  8 

IF(JM/X-J)202»202»201 

202  J=l 
SUMT=0.0 
ZB=0.0 

ZT=0.0  •  < 

x=o«o 

I  X=D(J) 
IF(X)203,203,204 

203  ZB=1.0+ZB 
X=-X 

GO  TO  205 

204  ZT=1.0+ZT 

205  SUMT=X+SUMT 
J  =  J  +  1 
IF(JMAX-J)2»1»1 

2  RT=ZT/SUMT 
R3=Z3/SUMT 

IF (SENSE  SWITCH  3) 10 v 11 
10  PRINT  103»M,RB»RT»SUMT 
GO  TO  12 

II  PUNCH  103,M,RB»RT»SUMT 
12  1=1+1 

GO  TO  200 
END 
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APPENDIX  C 

Description  and  Explanation  of  the  IH'l 
lo20  Computer  FORTRAN  Program  No.  3  for 
Performing  the  Fluid! za-oi on  Statistical 
Analysis 

This  program  was  written  to  determine  the  individual  and  average 
times,  variances,  and  skewness  of  the  experimental  data  obtained  from 
the  fluidized  bed  apparatus. 

Figure  93  is  an  illustration  of  the  different  time  increments  that 

can  be  obtained  from  a  set  of  experimental  data. 

The  average,  variance,  and  skewness  of  a  given  time  increment  are: 
H 

*  =    X   *iA  do) 

S2=  ]T(t±   -  t  )2/  flST-l)  (20) 

(-1 

By  utilizing  this  program,  individual  values,  averages,  variances,  and 
skewness  of  the  time  increments  t(BB),  t(BT),  t(BTL),  t(TT),  t(TB),  and 
t(TBL)  can  be  determined  from  the  experimental  data.   The  other  time  incre- 
ments shown  in  Fig.  93  can  be  determined  by  combining  the  output  of  this  code  in 
the  following  manner: 

t(T'TBB)  =  t(TT)  +    t(TB)  f  t(BB) 

t(BBTT)  —  t(BB)  -f  t(BT)  -f-  t(TT) 

t(TBT)   =  t(TB)  -f  t(BB)  +  t(3T) 

t(BTB)   ^  t(BT)  +  t(TT)  -f-  t(TB) 

t(BTT)   =  t(BT)  4-  tC^) 

t(TBB)   =  t(TB)  +  t(BB) 


J7-7 


pulses   f r .  ,n    upper 
detector    level 


Pulses    from    lower 
detector     level 


time    — 


Figure     93.      rime      Increments      AvcrliGble      From      a     Set 
of     Experimental      Data. 
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Sense  switches  on  the  IBM  control  console  route  the  computer  through 
the  program  to  yield  the  output  desired.   Sense  Switch  3  OH  gives  console 
typewriter  output  and  Sense  Switch  3  OFF  gives  punched  card  output.   Sense 
Switches  1  and  2  give  the  following  output: 

Sense  Switch  1     Sense  Switch  2  Output 

OFF  OFF  Individual  values,  t(BB), 

t(BT),  t(BTL),  t(TT), 
t(TB),  and  t(TBL) 

ON  OFF  Individual  values,  t(BT) 

and  t(TB) 

ON  ON  Individual  values,  t(BB), 

t(BT),  t(TT),  and  t(TB); 
average  values,  variances, 
and  skewne  s s  of  t ( B3 ) , 
t(3T),  t(BTL),  t(TT),  t(TB), 
and  t(TBL) 

OFF  ON  Individual  values,  average 

values,  variances,  and 
skewness  of  t(B3),  t(BT), 
t(B'IL),  t(TT),  t(TB),  and 
t(TBL) 


Table  12  defines  the  various  symbols  utilized  in  this  computer  code, 


Table  12.   Input  data  and  variables  required  for 
the  IBM  1620  Computer  Program  No.  3- 
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SYMBOL 


XPLAKATION 


JMAX 
D(J) 

I 
IT 

BB,  BT,  BTL 
TT,  TB,  TBL 
BBA,  BTA,  BTLA 
TTA,  TBA,  TBLA 
BBXR,  BTXR 
TTXR,  TBXR 
BTLXR,  TBLXR 
BBYR,  3TYR 
TTYR,  TBYH 
BTLYR,  TBLYR 


Number  of  time  increment  data  for  given  set 

Individual  time  increment  data  in  given  set, 

J  -  1  to  JMAX,  in  seconds 

Data  set 

Number  of  t(BT  or  t(TB)  in  data  set 

t(BB),  t(BT),  t(BTL) 

t(TT),  t(TB),  t(TBL) 

Average  values  of  t(BB),  t(BT),  and  t(BTL) 

Average  values  of  t(TT),  t(TB),  and  t(TBL) 

Variances  of  t(BB)  and  t(3T) 

Variances  of  t(TT)  and  t(TB) 

Variances  of  t(BTL)  and  t(TBL) 

Skevness  of  t(BB)  and  t(BT) 

Skewness  of  t(TT)  and  t(TB) 

Skewness  of  t(BTL)  and  t(TBL) 


Input  to  Computer  Program  No. 

The  input  to  Computer  Program  No.  3  is  the  output  of  Computer  Program 
No.  1.   The  input  is  read  into  the  computer  as  FORMAT  (i^,  Ik,   8F8.3). 

Output  of  Computer  Program  No.  3 


Tables  13  and  lk   give  the  output  of  Computer  Program  No.  3«   The  out- 
put given  in  Table  13  was  obtained  by  setting  Sense  Switch  1  ON  and  Sense 
Switch  2  OFF.   The  results  were  the  residence  times  t(BT)  and  t(TB)  of  the 
tagged  particle  as  it  passed  from  one  detector  level  to  the  other. 

Table  Ik   gives  the  output  of  the  program  when  Sense  Switch  1  was  ON 
and  Sense  Switch  2  was  ON.   The  results  give  the  individual  time  increments 
of  t(BB),  t(BT),  t(TB),  and  t(TT)  and  the  average  values,  variances,  and 
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Ta"ole    15.       Output  of  Ccputer  Program  Ho.    3. 


1  4  7  .02  7 

1  •  4 1  .  0  3  5 

1  49  .^63 

1  5  0  .120 

1  bl  .^90 

1  52  .065 

1  53  .273 

1  54  .107 

1  55  .122 

1  56  .109 

1  57  .096 
53 

2  1  .243 
2  2  .342 
2  3  .326 
2  4  .152 
2  5  .141 
2  6  .130 
2  7  .127 
2  .  8  .189 
2  9  .180 
2  10  .125 
2  11  .212 
2  12  .173 
2  13  .136 
2  14  .193 
2  15  .212 
2  16  .314 
2  17  .141 
2  '  18  .368 
2  19  .193 
2  20  .199 
2  21  .334 
2  22  .156 
2  2  3  .151 
2  24  .130 
2  25  .193 
2  27  .143 
2  28  .318 
2  29  .344 
2  30  .132 
2  31  .136 
2  32  .304 
2  33  .256 
2  34  .393 
2  35  .223 
2  36  .146 


52 

1    1 

.2  19 

1    2 

.223 

1     3 

.071 

1     A 

.096 

1    5 

.  123 

1    6 

.2  55 

!    H 

•  U72 

1    9 

•  U6  7 

1   1  v 

.  u  6  1 

1   1  1 

•  UG4 

1   1? 

.132 

1    14 

.  1  06 

1   1  5 

.14] 

1   1  6 

.  173 

1   1  7 

.  1  40 

1   ]  8 

•  U  8  3 

1   19 

.  1  08 

1   2  0 

.  102 

1   21 

.226 

1  22 

.0  9  7 

1   23 

.114 

I   2  6 

.  114 

1   2  7 

.  ]  00 

1   2" 

.109 

1   2  9 

.119 

1    3  ' ' 

.113 

1   ^  1 

.265 

1   32 

.256 

1   3  4 

.035 

1   35 

.124 

1   36 

.o27 

1   37 

.119 

1   38 

.217 

1   39 

.191 

1   40 

.078 

1   41 

.152 

1   42 

.173 

1   43 

.027 

1   44 

.237 

1   45 

.069 

1   46 

.227 

38 

•  200 

c 

3  9 

.2  37 

2 

41 

.404 

2 

4  2 

.  160 

2 

4  3 

.112 

2 

4  4 

.2  98 

2 

45 

.125 

2 

46 

.193 

2 

47 

.189 

2 

43 

.258 

2 

49 

.414 

2 

50 

.179 

2 

51 

.220 

2 

52 

.112 

2 

53 

.162 

2 

54 

.149 

2 

55 

.193 

2 

56 

.229 

36 

3 

1 

.110 

3 

3 

.088 

3 

5 

.071 

3 

6 

.093 

3 

7 

.154 

3 

8 

.162 

3 

9 

.106 

3 

12 

.^78 

13 

.103 

3 

14 

.130 

3 

16 

.121 

3 

17 

.027 

3 

13 

.027 

3 

19 

•  C51 

3 

20 

.164 

3 

21 

.  0  9 1 

3 

22 

.126 

3 

23 

.109 

3 

24 

.103 

3 

26 

.106 

3 

28 

.114 

Ta"ble  14.  Output  of  Computer  Program  17o.  3. 
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Data  Set     t(ii) 
31      3.23701920 
31     15.38825900 
31    135.90179000 


t(ij) 

.28257692 
.10038160 
.14487333 


t(ijL) 


3.50498070 

15.04018300 

.34.94166000 


32  8.05963460 
32  136.96316000 
32   3198.2905UOOO 


.23053846 
-07123977 
.06782903 


8.26719230 

136*49918000 

3199.70450000 


33  2.51266660 
33  26.91388100 
33    561.69701000 


2.731C1660 

30.41322800 

512.40540000 


3.68613330 

36.37274400 

763.84357000 


34  3.54749150. 
34  43.44949600 
34    899.12203000 


1.02400000 

18.66181200 

538.8125300.0 


4.25988130 

55.45750800 

1239.27460000 


35  1.72640290 
35  8.17724780 
35     63.60642500 


6.31608950 
146.51537000 

4793.33400000 


3.35334320 

33.68341000 

807.80359000 


36  4.18843930 
36  83.69991500 
36   2816.15760000 


.31922727 

.22380333 
.63194858 


4.33853030 

85.39653300 

2S63.51490C00 


37  6.20802040 
37  45.24609100 
37    555.72408000 


.44144897 
.07653613 
.01793786 


6.60061220 
45.49183100 
.66.78697000 


38  9.46933330 
38  174.97080000 
33   3877.21040000 


.30862500 
.02258699 
.00257611 


9.74420830 

175.37152000 

3898.11080000 


39  4.51019560 
39  27.86508600 
39    276.48140000 


.29517391 
•02092481 
.00365211 


4.77421730 

28.15774400 
285.51908000 


40  8.12597770 
40  105.52295000 
40   3312.57340000 


.41760000 

.65326877 

3.09779020 


8.53124440 

106.66721000 

3379.24220000 


41  5.10441660 
41  47.04307400 
41    908.10965000 


.66652777 
.07792275 
.03519333 


5-75194440 

47.69543100 

929.79262000 


42  8.71800000 
42  103.18281000 
42   1583.16050000 


.62762162 
.14275687 
.07150368 


9.31221620 

104.45862000 

1599.52760000 


J-G2 


skewness  of  t(BB),  t(BT),  t(ETL),  t(TT),  t(TB),  and  t(TEL). 

152  residence  times  sets  were  obtained  from  the  j6   Data  Sets  of  input. 

Tlie  second  column  in  each  set  of  columns  in  Table  13  gives  the  subscript 
assigned  to  each  individual  residence  time  in  each  given  residence  times  set. 
The  third  column  in  each  set  of  columns  gives,  in  seconds,  each  individual 
residence  time  in  each  given  residence  times  set. 

The  first  row  of  an  even- numbered  output  set  in  Table  1^  gives,  in 
seconds,  the  average  values  of  t(TT),  t(TB),  and  t(TBL),  respectively.   The 
second  row  gives,  in  (seconds)  ,  the  variances  of  t(TT),  t(TB),  and  t(TBL), 
respectively.   The  third  row  gives,  in  (seconds)  ,  the  skewness  of  t(TT), 
t(TB),  and  t(TBL),  respectively. 

The  first  row  of  an  odd- numbered  output  set  in  Table  Ik   gives,  in 
seconds,  the  average  values  of  t(B3),  t(BT),  end  t(BTL),  respectively. 
The  second  row  gives,  in  (seconds)  ,  the  variances  of  t(BB),  t(BT),  and 
t(BTL),  respectively.   The  third  row  gives,  in  (seconds)  ,  the  skewness  of 
t(BB),  t(BT),  and  t(BTL),  respectively. 


Logic  Diagram  for  Computer  Program  No.  3 
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Start 

Read    in    Dote  (J) 

and     J MAX 


iSot 

2DS  =  Q 
STS  sQ 

3B=    Q 

.31=    Q 

P  =1  .0 

BTLS»Q 

'esj-^ 

|Q  =  0.0 

Y 

«— <QJMAX>J  ^>-*-(  1 

/^V  No  y 

Set 

Point  = 
Datc(J) 

iQh 


o-j 


Sot 

J5    | 

1=2 
R  =  I.O 


t  BTX  = 
E7LX 


B  BY  =  Q 
3TY  =  Q 
B7LY  =  Q 


BBD=  Q 
2TD  =  Q 
BTLDsQ 

1 


SB=  Q 
BT  =  Q 
BTL  =  Q 


B3V  =Q 

BTV  =  Q 
3TLV=Q 


1  J  =  J  ♦  I    J 


R  =  R  + 1 .0 
BTX=BTD+BTX 

BBY=BBV+E3Y 
37Y=B7V  +  STY 
S7LY  =  37LV  +  B7LY 


B3X=33D-i-E3X 
BTD  =  (ST-3TA)'" 
37LX=E7LD+37LX 
3TV  =  (37-37A)° 


33D=  (3B-S3A)2 
BSY  =  (33-33A)5 
E7LD=(37L-37LA)21~S" 
S7LV=(37L-37LA)3 


37L  = 

B3  +BT 


Punch 

I,  BBXR,BTXR, 
37LXR, B3YR, 
37 YR,  B7LYR 


B7XR=  B7X/R 
BTLXR=S7LX/R 
B7LYR  =  37LY/R 
BBYR  =  BBY/R 
BTYR  =  37Y/R 


I 


R=R-2.0 
BBXRsBBX/rH 


Punch 

I,  B3X,B7X 

E7LX,  B3Y.B7Y, 

B7LY 


<-  I  =   13 
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TTS 

=  TTS 

+  TT 

TBS 

=  TES 

+  TB 

TBLS  =  T3LS,1 

+  TBL 
P=P  +  1.0 
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Sat 
J=  I 
1=2 
R  =  I.O 


TTX  =  0 
7  3X  =  Q 
73LX=0 


TTY=  Q 
T3Y  =   0 

T3LY=  Q 


TTO  =  Q 
T  3=0 
T3LD=Q 


r 


RsR  +  1.0 

7BV=(7B-7BA)° 
TTY=TTV+TTY 
7BY=73V  +  7BY 
73LY  =  7BLV+73LY 


73D  =  (73  -73A)' 
77X=77D  +77X 
7BX  =7SD-i-73X 
7ELX=7BLD-1-7BLX 


TBLD  =  (73L-7BLAV 
7BLV=(73L-7BLA)2 
77D  =  (77-77A}| 
77V=(77-77A)3 


Punch 

I,  TTXR.TBXR, 

7BLXR,  77YR, 
7BYR,  73LYR 


73XR=  TBX/R 

TBLXRsTBLX/R 

TTYR=TTY/R 

73YR=7EY/R 

TBLYR=TBLY/R 


R=  R-2.0 
77XR=77X/r[ 


Punch 

IJ7X.73X, 

7BLX,T7Y,73Y, 

7BLY 


16? 


FORTRAN  Source  Program  for  Computer  Program  2To.  3 


DIMENSION 

D(300) 

FORMAT ( 14, 

I4,F8.3i 

iF3.3 

,F8 

.3, 

F8 

FORMAT ( 14, 

F15.8, 

F15. 

8, 

F15 

.8 

FORMAT ( 15 ) 

FORMAT ( 14, 

14,  F7. 

.3) 

1=1 

READ  103,JMAX 

101  FORMAT ( I4»I4»F8.3fF3.3,F8.3»F8.3»F8.3,F8«3,F8.3»r8.3) 

102  F0RMAT(I4,  F15.8,  F15.8,  F15.8,  F15.8) 
103 
104 

200 

201  J=l 

202  READ  101,M,N,D( J) ,D( J+l ) »D( J+2) ,D( J+3) ,D( J+4) ,D( J+5) ,D( J+6) ,D( J+7) 
J  =  J  +  8 

IF( JMAX-JJ 1,1,202 

1  J  =  l 

2  N  =  l 

3  P  =  1.0 
Q  =  0.0 

4  5BS=Q 
5TS  =  Q 
5TLS=Q 

5  53  =  Q 
BT  =  Q 
BTL=0 

6  POINT=D(J) 
IF(POINT) 10,10,7 

7  J=J+1 

IF( JMAX-J )  19,6,6 

10  J=J+1 

I F( JMAX-J) 11,13,13 

11  IF(SENSE  5WITCH  1)19,203 

203  IF(SEN5E  SWITCH  3)300,801   . 

800  PRINT  10?, I ,BB,BT,8TL 
GO  TO  19 

801  PUNCH  102,1 ,BB,BT,BTL 
GO  TO  19 

13  TERM=D(J) 

I F( TERM)  14,15,15 

14  UB=BB-TERM 
GO  TO  10 

15  BT=TERM 

I F(BB) 16,17,16 

16  BTL=BB+BT 

17  IF(SENSE  SWITCH  2)18,204 

204  PUNCH  104,  I ,  N,  BT 
N  =  N  +  1 

GO  TO  5 
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IS  BBS=BBS+B8 
BTS=BTS+BT 

STLS=STLS+BTL 

P=P+1.0 

IF  (  SENSE  SWITCH 

PRINT  102>I,BB,3T 

60  TO  5 

PUNCH  102»  I*  B3»  3T 


205 

217 

19 

20 

21 

206 

218 

22 
23 


3  )  205»  217 


GO  TO  5 

IF(SENSE  SWITCH  2)20,37 

IF(P-1.0)22,21,22 

IF  (  SENSE  SWITCH  3  )  206,  218 

PRINT  102,  I,  BBS*  BTS,  BTLS 

GO  TO  24 

PUNCH  102,  Iv  BBS,  BTS,  BTLS 

GO  TO  24 

P=P-1.0 

BBA=(3BS/P) 

5TA=(BTS/P) 

BTLA=(BTLS/P) 

IF  (  SENSE  SWITCH  3  )  207,  219 
207  PRINT  102,I,B3A,BTA,BTLA 

GO  TO  24 
219  PUNCH  102,  I,  BBA,  BTA,  BTLA 

24  J  =  l 

R  =  1.0 

33X  =  Q 

BTX  =  Q 

3TLX=Q 

33Y=Q 

BTY=0 

BTLY=Q 

25  33D=Q 
BTD=Q 
BTLD=Q 
BBV=Q 
BTV=Q 
BT!_V  =  Q 
BB  =  Q 
BT  =  Q 
3TL=Q 

26  POINT=D(J) 
IF(POINT)28,28,27 

27  J  =  J  +  1 
IF(JMAX-J)34,26,26 

28  J=J  +  1 
IF(JKAX-J)34, 29.29 

29  TERM=D(J) 
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30 
31 
32 


33 


34 

35 

208 


220 


36 


209 


221 

37 
45 


IF(TERM)3G,31,31 

B3=BB-TERM 

GC    TO    23 

BT=TERM 

IF(3B;32>33>32 

BTL=B3v3T 

63D=(BB-B3A)  **2 

t3BV=<B6-bBA)  **3 

BTLD=(BTL-BTLA)**2 

BTLV=(BTL-BTLA)**3 

BBX=BBD+BBX 

BTD=(BT-BTA)**2 

BTLX=BTLD+BTLX 

R  =  R  +  1.( 

BTV=(0T-BTA)**3 

BTX=BTD+BTX 

BBY=BBV+BBY 

BTY=BTV+BTY 

BVUY=BTLV+BTLY 

GC  TO  25 

IF(R-2. 

IF  (  SE 

PRINTIU 

PRINT  1 

GC  TC  3 

PUNCH  1 

PUNCH  1 


GC  TC  3 
R=R-2.0 
BBXR=BB 
BTXR=BT 
BTLXR=B 
BBYR=BB 
BTYR=BT 
BTLYR=B 
IF  (  SE 
PRINT  1 
PR  IN'.  1 
GC  TC  4 
PUNCH  1 
PUNCH  1 
GC  TC  4 
J  =  l 
1  =  1+1 
N=l 
P  =  1.0 
TTS  =  Q 
T3S=Q 
TBLS=Q 


0) 36,35,36 

NSE  SWITCH  3  )  208,  220 

2,1 ,BBX,BTX»3TLX 

U2, I  ,BBY,BTY*3TLY 

7 

02,  I»  BBX,  BTX,  BTLX 

02,  I,  BBYt  BTY,  3TLY 

7 


X/R 
X/R 
TLX 
Y/R 
Y/R 
TLY 
NSE 
0?  . 

02, 

5 

02, 

02, 

5 


/R 


/R 

SWITCH  3  )  209,  221 
I >33XR,BTXR,BTLXR 
I,BBYR,BTYR,BTLYR 


I,  bby;^, 


BTXR, 
BTYR, 


3TLXR 
BTLYR 


1  cr> 


46 

TT  =  Q 

T3  =  Q 

TBL  =  Q 

47 

PCINT=D( J) 

IF(PCINT)4Ss>4S»49 

4S 

J  =  J+1 

IF(JMAX-J)58»47»47 

49 

J=J+1 

IF( UMAX- J) 50,52,52 

50 

IFCSENSE  SWITCH  1)51»210 

210 

PUNCH  102,I»TT,TB,TBL 

51 

60  TO  58 

52 

TERM=D( J) 

IF  (TERM)  54,53 »53 

53 

TT=TT+TERM 

GC  TC  49 

5  4    TERM=-TERM 
T3=TERM 
IF(TT)55,56»55 

55  TBL=TT+TB 

56  IFCSENSE  SWITCH  2)57,211 

211  PUNCH  104,I,N,TB 
N  =  N+1 

GO  TC  46 

57  TTS=TTS+TT 
TBS  =  TBi.  +  TB 

1   TSLS  =  TBLS  +  TBL 
P=P+1.0 
IF  (  SENSE  SWITCH  3  )  212,  222 

212  PRINT  102,  I,  TT,  TB 
GC  TO  46 

222  PUNCH  102,  I,  TT,  TB 
GO  TO  46 

58  IFCSENSE  SWITCH  2)59,75 

59  IF(P-1»0)61,60,61 

60  IF  (  SENSE  SWITCH  3  )  213,  223 

213  PRINT  102,  I,  TTS,  TBS,  TBLS 
GO  TO  62 

223  PUNCH  102,  I,  TTS,  TBS,  TBLS 
GO  TO  62 

61  P=P-1»0 
TTA=TTS/P 
TBA=TBS/P 
T3LA=TBLS/P 

IF  (  SENSE  SWITCH  3  )  214,  224 

214  PRINT  102,I,TTA,TBA,TBLA 
GO  TO  62 

224  PUNCH  102,  I,  TTA,  TBA,  TBLA 
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62 


63 


64 


65 


66 


67 


68     TT  = 


=  1 
=  1 
TX 
BX 
BL 
TY 
DY 
BL 
TD 
BD 
BL 
TV 
BV 
BL 
T  = 
B  = 
BL 
CI 
F( 
=  J 
F( 
=  J 
F( 
ER 
F  ( 


.0 

=  0 

=  0 

X=Q 

=  0 

=  0 

Y-Q 

=  Q 

-Q 

D  =  0 

=  0 

=  0 

V  =  Q 

0 

Q 

=  Q 

N  T  -  0  (  J  ) 

P0INT)65,65»66 

+  1 

JMAX-J )72»64,64 

+  1 

JMAX-J) 72*67,67 

M=D( J) 

1ERM)69»68,68 

TT+TERM 

TO  66 


69  TERM=-TERM 
TB=TERM 
IF(TT)70,71,70 

70  TBL=TT+TB 
TBLD=(TBL-TBLA)**2 
TBLV=(TBL-T3LA)**3 
TTD=(TT-TTA)**2 
TTV=(TT-TTA)**3 

71  R=R+1«0 
TBD=(TB-TBA)**2 
TTX=TTD+TTX 
T3X=TBD+T3X 
TBLX  =  T3LD-i-TBLX 
TBV=(TB-7BA)**3 
TTY=TTV+TTY 
TBY=TBV+TBY 
7BLY=TBLV+TBLY 
GO  TC  63 

72  IF(R-2. 0)74, 73.74 


1*2 


73 
215 


225 


74 


216 


226 


76 


IF  (  SEN 
PRINT  10 

PRINT 
60  TO 
PUNCH 
PUNCH 
60  TO 


S£  SWITCH  3  )  215, 
2, I ,7TX,TBX,TBLX 
2,I,TTY,T3Y,T3LY 


225 


10 

76 

io; 

102,  I,  TTY,  T5Y,  TBLY 
7  6 


I,  TTX,  TBX,  T3LX 


75 


R-R-2-0 

TTXR=TTX 

T3XR=T3X 

TBLXR=TB 

TTYR=TTY 

TBYR=TBY 

TBLYR=TB 

IF  (  SEK 

PRINT  10 

PRINT  iO 

CO  TO  76 

PUNCH  10 

PUNCH  10 

TTA=Q 

TBA  =  Q 

TBLA=0 

TTXR=Q 

TBXR=Q 

TBLXR=Q 

TTYR=Q 

T3YR=Q 

TBLYR=Q 

B3A  =  Q 

BTA=Q 

BTLA=Q 

3BXR=Q 

3TXR=Q 

3TLXR=Q 

3BYR=Q 

3TYR=Q 

3TLYR=Q 

1=1  +  1 

60  TO  20 

END 


/R 
/R 

LX/R 

/R 

/R 

LY/R 

SE  SWITCH  3  )  216,  226 

2,1 ,TTXR,TBXR,T3LXR 

2,I,TTYR,TBYR,TSLYR 

2,  I,  TTXR,  T3XR,  TBLXR 
2,  I,  TTYR,  TBYR,  TBLYR 
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APPENDIX  D 


Description  and  Explanation  of  the  IBM 
1620  Computer  FORTRAN  Program  No.  k   for 
Determining  the  Concentration  of  Tagged 
Particles  from  Individual  Particle  Data 


This  computer  program  determines  the  values  of  IMM,  HIP,  IPP,  and 
IPM  from  the  experimental  data.  Figures  82  through  85  illustrate  results 
Obtained  from  this  program.   Table  15  defines  the  various  symbols  uti- 
lized in  this  computer  code. 


Table  15.   input  data  and  variables  required  for 
the  IEM  lo20  Computer  Program  No.  k. 


SYMBOL  EXPLANATION 

1  i  IAX  Number  of  time  increments  At 

T(l)  Ith  time  increment,  in  seconds 

JI  IAX  Number  of  time  increment  data  for  given  set 

d(j)  Individual  time  increment  data  in  given  set, 

J  -  1  to  JMAX,  in  seconds 
POINT  Initial  data  pulse 

TnJi(l)  Number  of  negative  pulses  within  time  increment 

T(l),  starting  with  an  initial  negative  pulse 
IMP(l)  Number  of  positive  pulses  within  time  increment 

T(l),  starting  with  an  initial  negative  pulse 
IPP(l)  Number  of  positive  pulses  within  time  increment 

T(l),  starting  with  an  initial  positive  pulse 
IPM(l)  Number  of  negative  pulses  within  time  increment 

T(l),  starting  with  an  initial  positive  pulse 


Input  to  Computer  Program  No. 

The  data  input  to  Computer  Program  No.  k   is  the  output  of  Computer 
Program  No.  1,  along  with  the  time  increments  T(l).  The  data  input  is  read 
into  the  computer  as  FORMAT  (l1!-,  Ik,   8F8.3).   The  time  increments  T(l)  are 
read  into  the  computer  as  FORMAT  (l4,  Ik,   5E7-3). 
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Tac  total  number  of  tine  increments  IMAX  is  read  into  the  computer  as 
FORMAT  (l5)«   Table  l6  gives  the  time  increments  T(l)  and  the  number  of  time 
increments  IMAX  used  in  evaluating  the  data  sets  for  each  distance  between 
detector  levels  AZ. 

The  first  number  in  each  set  of  time  increments  T(l)  for  given  distances 
between  detector  levels  AZ  is  the  total  number  of  time  increments  IMAX  for  the 
data  set.   The  first  column  below  IMAX  gives  the  HEM  card  number,  and  the 
second  column  gives  the  subscript  I  for  the  first  time  increment  T(l)  on  the 
IBM  card.   The  remaining  columns  give  in  seconds,  T(l),  T(l-fl),  T(l4-2), 
T(l-f-3),  and  T(l+  4),  respectively. 

Output  of  Computer  Program  Ho.  h 

Table  17  gives  the  output  obtained  from  Computer  Program  No.  h.      The 
method  of  tabulating  IMM,  IMP,  IPP,  and  IPM  for  AZ  —  1/2  inch  is  indicated 
in  Set  1. 

Each  row  of  the  tabulated  concentrations  (iMMj  IMP,  IPP,  and  IPM)  gives 
the  number  of  data  pulses  counted  in  the  time  increments  T(l),  T(l-f-l),  T(l-J-2), 
T(l-f  3);  and  T(l+  k) ,   beginning  with  I  —  1.   Blank  spaces  in  the  block- 
tabulations  of  the  output  of  each  data  set  indicate  zeros. 

Table  l6.   Input  to  Computer  Program  No.  h. 


IMAX 

10 

Card 

1 
2 


T(I) 

0.020 
0.300 


T(l-fl) 

O.QkG 
0.500 


T(l-f  2) 

0.060 
1.000 


T(l+3) 

0.080 
2.000 


T(l-f^) 
0.100 

3.000 


Tn'ole   17.  Output  of  Computer  Program  ICO.   4. 

AZ  =  -V  inch 
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Logic  Diagram  for  Computer  Program  Ho.  4 


— >  - 


IM  =0 
I?  =0 
Sum  =  0.0 


Point  = 
Dcta  (J) 


•97 


(Punch 
IMM(I),IMP{I) 
IPP,  IPM(IJ 


1 

i 
/ 

Set 

IMMllJa  0 
IMP(I)s  0 

IP?(I)=0 
IPM{I)  =  0 

END 


3-^6 

K)RTRA2\T  Source  Program  for  Computer  Program  No.  4 

PROGRAM  FOR  INCREMENTAL  SUMMATIC/v 
DIMENSION  T(50)  ,D(600)  ,  IMM(  50)  ,  i  .>„J  (  50  )  ,  I  PM  (  50  )  >  IPP  (  SO  ) 

101  F  ORMAT (1 4 , 1 4 , F8 . 3  >  F8 . 3 , F8 • 3  v  F 0  -  3 , F8 . 3 , F8 . 3 , F8 . 3 , F8 „ 3 ) 

102  FORMAT (  14 t  14,14,14,14) 

103  FORMAT ( 15 ) 

104  FORMA T< I4»I4,E7.3»E7.3»E7.3»E7.3»E7.3»E7.3»E7.3>E7.3) 

209  RLAO  HO,  I  MAX 
1=1 

200  READ  iu4,K,L,T( I ) ,T( 1  +  1) ,T( 1+2) »T( 1+3) ,T( 1+4) 
1  =  1+5 

IF(  IMAX-I  )210,210,200 

210  READ  1U3,  JMAX 

201  J=l 

202  READ101»M,N,D(J),D(J+1) ,D( J+2) ,D( J+3) ,D(J+4) ,D(J+5) ,D( J+6) ,D(J+7) 
J  =  J  +  8 

IF( JMAX- J) 20 3,2 03 » 202 

203  1=1 

I  1=1 
J  =  l 
2  0  IM*0 
IP  =  0 
SUM=U.U 

1  POINT-D(J) 

2  J-J+l 

IF( JMAX-J)3>3,4 

4  TERM=D(J) 

I  F( TERM) 5, 5, 7 

5  SUM=SUM-TERM 

6  IM=1+IM 
60  TO  9 

7  SUM=SUM+TERM 

8  IP=1+IP 

9  DELTA=T( I ) 
IF(DELTA-SUM)10,2»2 

io  if(Ter:.>ii9ii,i2 

11  IM=IM-1 
IF(P0INT)13,13,14 

12  I?=IP-1 
IF(POINT)15,15,16 

13  IMM( I )=IM  +IMM( I ) 
IMP( I )=IP  +IMP( I ) 
GO  TO  17 

14  IPM( I )=IM  +IPM( I) 
IPP(I)=IP  +  IPP(I ) 
GO  TO  17 
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15 

16 

17 
18 

19 


400 
206 


300 


500 

205 
207 
208 


IMP( I )=IP  +IMP( I ) 

I  MM ( I )  =  I M  +  I MM ( I ) 

GO  TO  17 

I PP  C 1 )  =  I P  +  IPPU  ) 

IPM(I)=IM  +  IPM( I ) 

GO  TO  17 

1*1+1 

IF(  IMA;.-I  )  18,19, 19 

11=11+1 

1  =  1 

J=I  I 


GO 
IM  = 
IP= 
IF( 

1  =  1 
IF( 
PRI 
PR  I 
PRI 
PRI 


TO  20 

0 

0 

TERM) 6, 6, 8 


310 
211 


SENSE 
NT102* 
NT102* 
NT1029 
NT102> 
GO  TO  500 
PUNCH102» 
PUNCH102» 
PUMCH102* 
PUNCH102, 
1  =  1+5 
IF(IMAX-I 
IF  (SENSE 
1=  1 
I  MM ( I ) 
IMP(  I  ) 
IPP(I) 
I PM (  I  ) 
I  =  1  +  1 
IFC IMAX-I 
IF(SENSE 
END 


SWITCH  3 ) 206 , 300 

I  MM  (  I  )  9  IMM  (  1  +  1 )  9  I  MM  (  I  +2  )  »  I  MM  ( I +3  )  »  I  MM  (  I  +4  ) 
IMP( I) »IMP( 1+1) »IMP( 1+2) 9 IMP ( 1+3) »IMP( 1+4) 
IPP( I ) »IPP( 1+1) ,IPP( 1+2) ,  IPP( 1+3) ,  IPP( 1+4) 
IPM( I ) »IPM( 1+1) »IPM( 1+2) »IPM( 1+3) ,IPM( 1+4) 


)205,400,400 
SWITCH  1)  207* 


211 


=0 

=  0 

=c 

=0 


)310»203,203 
SWITCH  2)209,210 
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APPENDIX  E 


Description  and  Explanation  of  the  IEM 
1620  Computer  FORTRAN  Program  No.  5  for 
Calculating  the  Average  Residence  Times 
as  a  Function  of  the  Number  of  Data  Val- 
ues 


Tnis  computer  program  determines  the  average  time  for  a  particle 
to  pass  from  upper-to-lower  and  lower-to-upper  detector  level,  as  a  function 
of  the  number  of  data  values  obtained. 

Figures  35  a-^cl  3&  are  illustrations  of  the  results  obtained  from  this 
computer  code.   Table  18  defines  the  various  symbols  utilized  in  this  code. 

Table  18.   Input  data  and  variables  required  for 
the  IEM  1620  Computer  Program  No.  5. 

SYMBOL  EXPLANATION 

JMAX  Number  of  time  increment  data  for  given  set 

D(J)  Individual  time  increment  data  in  given  set, 

J  —  1  to  JMAX,  in  seconds 
M  Data  set 

N  Number  of  data  values  used  to  calculate  average 

time 
SUMA  Average  time  increment  for  given  number  of  data 

values 

After  the  data  in  each  set  is  read  into  the  computer,  the  average  time 
increment  is  computed  for  the  first  five  time  increments,  and  this  time  is 
punched  as  output  along  with  the  data  set  N.   Then  the  number  of  data  values 
N  is  increased  by  five,  and  the  process  is  repeated.  N  is  continously  in- 
creased by  increments  of  five  until  all  of  the  data  in  the  set  has  been  used. 
The  program  then  continues  to  the  next  set  of  data. 
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A  similar  code,  Computer  Program  No.  6,  computes,  as  a  function  of 
the  number  of  data  values,  the  variance  in  the  time  to  pass  from  one 
detector  level  to  the  other.   The  input  and  output  of  Computer  Program 
No.  6  is  included  in  this  appendix. 

Input  to  Computer  Program  No.  5 

The  input  to  Computer  Program  No.  5  is  the  output  of  Computer  Program 
No.  3>  obtained  with  Sense  Switch  1  ON  and  Sense  Switch  2  OFF.   The  input 
is  read  into  the  computer  as  FORMAT  {lh,   Ik,   KJ.k). 

Output  of  Computer  Program  No.  5 

The  output  of  Computer  Program  No.  5  is  given  in  Table  19 •  M  is  the 
residence  times  set,  where  odd-numbered  values  of  M  give  the  residence  tikes 
between  the  two  detector  levels  when  the  tagged  particle  passed  from  the 
lower  to  the  upper  detector  level,  and  even- numbered  values  of  M  give  the 
residence  times  between  the  two  detector  levels  when  the  tagged  particle. 
passed  from  the  upper  to  the  lower  detector  level.  N  gives  the  number  of 
residence  time  values,  in  order  of  appearance  in  the  residence  times  set, 
that  were  used  to  compute  the  average  residence  time  SUMA. 

Input  to  Computer  Program  No. 

The  input  to  Computer  Program  No.  6  is  the  output  of  Computer  Program 
No.  3,  obtained  with  Sense  Switch  1  ON  and  Sense  Switch  2  OFF.   The  input 
is  read  into  the  computer  as  FORMAT  (l4,  Xk,   E7.^). 


Table 

19.  Output  of 

Comput 

er  ?r 

ogram  l'.ro.  5 

:•: 

:: 

AVERAGE 

M 

IT 

AVERAGE 

i 

5 

.14640000 

6 

/ 

.47620000 

: 

10 

•12710000 

6 

.  j 

.44430000 

i 

15 

.13036666 

6 

.413133  3  3 

i 

20 

.12395000 

6 

■ 

.42895000 

i 

25 

.12576000 

6 

£   5 

.40332000 

i 

30 

.13123333 

6 

,0 

.33133330 

i 

35 

. 13054235 

6 

35 

.40297143 

i 

40 

.13067500 

6 

36 

.40300000 

i 

45 

.12664444 

i 

50 

.12712000 

7 

5 

.14900000 

i 

52 

.12617307 

7 

10. 

.13540000 

7 

15 

.17100000 

2 

5 

.24080000 

7 

20 

.18820000 

2 

10 

.19550000 

7 

25 

.19048000 

2 

15 

.19240000 

7 

30 

.19356666 

2 

20 

.20505000 

7 

33 

.19090909 

2 

25 

.20260000 

2 

30 

.20460000 

8 

5 

.38830000 

2 

35 

.21314285 

8 

10 

.37430000 

2 

40 

.21632500 

8 

15 

.32866666 

2 

45 

.21602222 

8 

20 

.30215000 

2 

50 

.21616000 

8 

25 

.30563000 

2 

53 

.21469310 

8 

30 

.31150000 

8 

34 

.32000000 

3 

5 

.10320000 

3 

10 

.10950000 

9 

5 

.27100000 

3 

15 

.09900000 

9 

10 

.35760000 

3 

20 

.10100000 

9 

15 

.33400000 

3 

25 

.10308000 

9 

20 

.36295000 

3 

30 

.10353333 

9 

25 

.35332000 

3 

35 

.10348571 

9 

30 

.37343330 

3 

36 

.10219444 

10 

5' 

.61640000 

4 

5 

.26920000 

10 

10 

.60630000 

4 

10 

.25860000 

10 

15 

.65160000 

4 

15 

.23886666 

10 

20 

.71670000 

4 

20 

.24155000 

10 

25 

.70312000 

4 

25 

.23120000 

10 

30 

.70056660 

4 

30 

.22920000 

10 

33 

.69724240 

4 

35 

.23600000 

4 

40 

.23900000 

11 

5 

.31280000 

4 

4  2 

.23950000 

11 

10 

.31550000 

11 

15 

.4 1646666 

5 

5 

.31620000 

11 

17 

.44470588 

5 

10 

.25930000 

5 

15 

.26113333 

5 

20 

.24760000 

5 

25 

.23996000 

5 

30 

.24480000 

5 

35 

.23648571 

5 

36 

.23676315 
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Output  of  Computer  Program  No.  6 

The  output  of  Computer  Program  No.  6  is  given  in  Table  20 ,   where  M 
and  N  are  defined  as  given  in  the  output  of  Computer  Program  No.  5-   The 
variance  SIC-SQ  for  the  number  of  residence  time  values  N  is  tabulated. 

Figures  37  and  38  give  illustrations  of  the  variance  of  the  residence 
times  as  a  function  of  the  number  of  residence  time  values  N  used  to  cal- 
culate the  variance. 


Table  20.   Output  of  Computer  Program  ITo.  6 
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2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
3 
3 
3 
3 
3 
3 
3 
3 

4 
4 

4 
4 
4 
5 
5 
5 
5 
5 
•5 
5 
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N 

5 
10 
15 
20 
25 
30 
35 
40 
45 
50 
52 

5 
10 
15 
20 
25 
30 
35 
40 
45 
50 
53 

5 
10 
15 
20 
25 
30 
35 
36 

5 
10 
15 
20 
25 
30 
35 
40 
42 

5 
10 
15 
20 
25 
30 
35 
38 


VARIANCE 

.00497780 

.00566743 

•003S3783 

.00355268 

.002S57S5 

.00387735 

.00403537 

.00423863 

.00451023 

.00460018 

.00444430 

.00883770 

.00664961 

.00457854 

.00574773 

.00569958 

.00630873 

.00682906 

.00734109 

.00691202 

.00731528 

.00699130 

.00099970 

.00093561 

.00192171 

.00146221 

.00122699 

.00139274 

.00132149 

.00134376 

.01185770 

.00799137 

.00731055 

.00883668 

.00876775 

.00896375 

.00892323 

.00886307 

.00863880 

.03651370 

.02399801 

.01810298 

.01614593 

.01448070 

.01779409 

.01688796 

.01787921 
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11 
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17 
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12 

10 

12 

15 

12 

18 

13 
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13 

10 

13 

15 

14 

5 

14 

10 

14 

15 

14 

17 

VARIAIICS 

.02674370 

.0218  3445 

.02643670 

.02732663 

.02697500 

.02514906 

.02945702 

.02861545 

.00135050 

.00706643 

.00657771 

.00656922 

.00584151 

.00593418 

.00573783 

.01336320 

.01079690 

.01204910 

.01170592 

.01193522 

.01035315 

.00935727 

.01566250 

.02335093 

.02241500 

.02592131 

.02347464 

.02610034 

.03795480 

.05159001 

.09939740 

.11630368 

.10579545 

.10401558 

.09838428 

.01233220 

.01539516 

.03829427 

.04332747 

.03045280 

.08997445 

.20039528 

.21437753 

.06373475 

.08850244 

.09570440 

.05018050 

.03023684 

.09274085 

.14066187 


Lo^ic  Diagram  for  Computer  Program  No.  5 
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Start 
Read 
J  MAX,  Dcta(J),M 


S  U  M  =0.0 
SU.\;A  =  C.O 
N  =  5 
J  =  I 


c 


Y  =  N 


1 


Calculate 

suma  = 
suma/y 


Punch 

M,    N,  SUM A 


N   =   N  +  5 


END 


Punch 

M,  J,  SUM  A 


Point    = 
Data  (J) 


SUM    = 

P  4-  SUM 


J  +■  I 


c 

c 

101 

102 

103 

200 

201 

206 
3JVDRTRAN  Source  Program  for  Computer  Program  No,    5 


PROGRAM  FOR  COMPUTING  THE  AVERAGE  TIME  AS  A  FUNCTION  OF  THE 

NUMBER  OF  DATA  POINTS  FOR  THE  BOTTOM  TO  TOP  AND  TOP  TO  BOTTOM  DATA 

DIMENSION  D(300) 

F0RMAT( I4»I4»E7.4) 

FORMAT ( 15 ) 

FORMAT (  15.15, F15. 3} 

1=1 

READ  102.JMAX 

J  =  l 

READ  l0l,M,N,D(J) 

J  =  J+1 

IF( JMAX-J) 1.201,201 

1  J=l 
SUM =0.0 
SUMA=( ,0 
N  =  5 

2  P=D(J) 
SUM=P+SUM 
J=J  +  1 

.  IF( JMAX-J)5,3.3 

3  IF<N-J)4.2»2 

SUMA=SUM/Y 
PUNCH  103.M»N,SUMA 
N  =  N + 5 
GO  TO  2 
5  J=J-1 
X=J 

SUMA=SUM/X 
P;^NCH  1G3,M,J,SUMA 
1=1  +  1 
GO  TO  200 
END 
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l^ORTRAU  Source  Program  for  Computer  Program  ITo.    6 


C      PROGRAM  FOR  COMPUTING  THE  VARIANCE  IN  THE 
C      AVERAGE  TIME  AS  A  FUNCTION  OF  THE  NUMBER 
C      OF  DATA  POINTS  FOR  THE  BOTTOM  TO  TOP  AND 
C      TOP  TO  30TT0M  DATA 
DIMENSION  DC  300) 

101  FCRMATt I4,I4,E7.4) 

102  FORMAT ( 15 ) 

103  FORMAT ( I59I59F15.S) 

200  READ  102,JMAX 
DO  201  J=1,JMAX 

201  READ  101,M,N»D( J) 

1  J  =  l 
SUM=0*0 
SUMSQ=0„0 
SIGSG=0.0 
N  =  5 

2  P=D(J) 
Q=(D(J)**2) 
SUM=P+SUM 
SUMSQ=Q+SUMSQ 
J=J  +  1 

IF( JMAX-J)5»3»3 

3  IF(N-J)4»2»2 

4  Y  =  N 

SIGSQ=(SUMSG-( <SUM**2)/Y> }/(Y-1.0) 
PUNCH  103»M*N»SIGSQ 
N  =  N+5 
•GO  TO  2 

5  J^J-1 
X  =  J 

SIGSQ=(SUMSQ-( (SUM**2)/X) )/(X-1.0) 
PUNCH  103»M»J»SI6SQ 
GO  TO  200 
END 
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APPENDIX  F 


Description  and  Explanation  of  the  IBM 
lo20  Computer  FORTRAN  Program  No.  7  for 
Computing. the  Average  Time,  the  Most 
Probable  Time,  the  Variance,  and  the 
Gamma  Distribution  Parameters  p  and  V 
of  the  Residence  Times  Data 


'This  computer  program  was  written  to  determine  the  average  time,  the 
most  probable  time,  the  variance,  and  the  gamma  distribution  parameters 
p  and  V  for  the  experimental  times  obtained  when  the  tagged  particle  passed 
from  one  detector  level  to  the  other,  moving  both  upward  and  downward. 
Table  21  defines  the  various  symbols  utilized  in  this  computer  code. 


Table  21.   Input  data  and  variables  required  for 
the  IBM  1620  Computer  Program  No.  7. 


SYMBOL 


EXPLANATION 


JMAX 

D(J) 

M,  P 

BA 

PR  and  V 

TP 

P 


Number  of  time  increment  data  for  given  set 
Individual  time  increment  data  in  given  set, 
J  =  1  to  JMAX,  in  seconds 

Data  set,  Number  of  data  in  set 

Average  time  for  data  set 

Gamma  distribution  parameters  p  and  u 

Most  probable  time 

Number  of  time  increment  data  for  given  set 


A  similar  code,  Computer  Program  No.  8,  computes  the  average  speed,  the 
variance  of  the  average  speed,  and  the  gamma  distribution  parameters  p  and  V, 
based  upon  the  Residence  Times  data.   The  input  and  output  of  Computer  Pro- 
gram No.  8  is  included  in  this  appendix. 
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Input  to  Computer  Program  ITo.  7 

The  input  to  Computer  Program  No    is  the  output  of  Computer  Program 
No.  3,  obtained  with  Sense  Switch  1  C   .1  Sense  Switch  2  OFF.   The  input 
is  read  into  the  computer  as  FORMAT  (.  '<,  Ik,   E7.4). 

Output  of  Computer  Program  No.  7 

The  output  of  Computer  Program  No.  7  is  given  in  Table  22.   Input 
data  information  tabulated  in  Table  22  are  as  follows.   The  set  number  M, 
as  explained  in  the  output  section  of  Computer  Program  No.  5  in  Appendix  E, 
is  given  for  each  Residence  Times  set,  1  to  152.   The  distance  between  the 
detector  levels  (WIDTH)  and  the  water  flow  rate  (INCHES  MERCURY)  are  also 
given. 

Tabulated  outputs  for  each  Residence  Times  set  are  as  follows. 

Average  residence  time  (AVERAGE  TIME) 

Most  probable  residence  time  (MOST  PROBABLE  TIME) 

Gamma  distribution  parameters  p  and  v   (PR  &  V) 

Variance  of  the  residence  times  (SIGMA  SQ) 

Speed  of  the  tagged  particle  based  upon  the  average 

residence  time  (VELOCITY  BASED  ON  T  AVG) 

Speed  of  the  tagged  particle  based  upon  the  most  probable 

residence  time  (VELOCITY  BASED  ON  T  PROB) 

Input  to  Computer  Program  No.  8 

The  input  to  Computer  Program  No.  8  is  the  output  of  Computer  Pro- 
gram No.  3,    obtained  with  Sense  Switch  1  ON  and  Sense  Switch  2  OFF.   The  input 
is  read  into  the  computer  as  FORMAT  (ik,    Ik,   F7.3). 


Table  22 


Output  of  Computer  Program  lTo.  7 


SET  NO    1        52 

W I D  T  H  = 

SIGMA  SQ= 

<Vv£RAGE  TIME* 

VELCCI  TY  BASED 

VFLOC1TY  BASED  ON  T  PRCB= 


DATA  POINTS 


.50000000 
.00444430   V= 
.12617307 
ON  T  AVG  = 


If   sES  MERCURY  = 
•C  ' 22383   PR= 
MOST  PROBABLE  TIME= 

^'626.  i  .'0   STD  DEV  = 
5.49757300   STD  DEV= 


2.50000000 

3.58203600 
.09004922 
2.09331510 
4.02970900 


SET  NO    2 
WIDTH= 

SIGMA  SO= 

average  time= 
velocity  based 
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DATA  POINTS 


.50000000 
.00699128   V= 
.21469810 
ON  T  AVG  = 


VFLOCITY  BASED  ON  T  PR03= 


INCHES  MERCURY55 
.03256334   PR= 
MOST  PROBABLE  TIME= 
2.32835150   STD  DEV= 
2.74522060   STD  DEV= 


2 . 5G000C  CO 
6.59324400 
.18213472 
.90696840 

1.26026360 


SET  N! 


3 


36.   DATA  POINTS 


WIDTH= 

SIGMA  SQ= 

average  time- 
velocity  BASED 

VELOCITY  BASED 


.5UOOO000   . 
.00134376   V= 
.10219444 
ON  T  AVG  = 
ON  T  PROB= 


INCHES  MERCURY5 
.01314905   PR= 
MOST  PROBABLE  TIME= 
4.39263400   STD  DEV5 
5.61511500   STD  DEV: 


2.50CC0CCC 
7.77199600 
.08904536 
1.75499710 

2.31157600 


SPT  NO    4 
W  I D  T  H  = 

SIGMA  S0= 

average  time= 
velocity  based 
velocity  based 


42.   DATA  POINTS 


.50000000 
.00863879   V= 
.23950000 
ON  T  AVG  = 
ON  T  PR03= 


INCHES  MERCURY5 
.03607011   PR= 
MOST  PROBABLE  TIME= 
2.03768260   STD  DEV: 
2.45734930   STD  DEV: 


6.63  934600 
.20342953 

.81013730 
1.12296700 


SET    NO         5 
WIDTH555 
SIGMA    SQ= 

average  time= 


38.   DATA  P( 
1.00000000 
.01787920   V= 
.23676315 


VELOCITY  BASED  ON  T  AVG  = 
VELOCITY  BASED  ON  T  PROB= 


:nts 

inches  mercury5 

.07551516   PR= 
MOST  PROBABLE  TIME555 
4.22363000   STD  DEV: 
6.20162300   STD  DEV^ 


2-500CCOOO 
3.13530600 
.16124797 
2.38531500 
5. 14263300 


SET  NO 


36. 


DATA  POINTS 


WIDTH*         1.00000000 

sigma  sg=     .02861542  v= 
average  time=     .40300000 
velocity  based  on  t  avg  = 
velocity  based  on  t  prob= 


INCHES  MERCURY5 
.07100600   PR= 
MOST  PROBABLE  TIME= 
2.48138900   STD  DEV5 
3.01210300   STD  DE\f- 


2.5000CC00 

67557600 

.33199395 

1.04157350 

1.53475730 


DATA  POINTS 


SET  NO    7        33 
WIDTH*         l.COOOOOCO 

sigma  sq=     .00578763  v= 
average  time=     .19090909 
velocity  based  on  t  avg  = 
velocity  based  on  t  prob= 


INCHES  MERCURY = 
.03031721   PR= 
MOST  PROBABLE  TIME= 
5.23809500   STD  DEV= 
6.22696500   STD  DEV= 


2.5000CCC0 
6.29705200 
.16059186 
2.0S739540 
2.94992450 
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Output  of  Compter  Program  No.  8 

The  output  of  Computer  Program  No  8  is  given  in  Table  23' 

Tabulated  outputs  for  each  residence  times  set  are  as  follows, 

Average  residence  time  (AVERAGE  TIME) 

Average  particle  speed  (AVERAGE  VEL) 

Gamma  distribution  parameters  p  and  V  (PR  &  V) 
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Table    25.         Output  of  Cc     v.  tor  Program  No.   8 
SET    NC=  1  5  2.  DAT/-    i    <IHTS 

width  =  -50    inches  mercury*  2.50 

v=     6.08539890  pr=           1.83532800  sgima  sq=  67.96605000 

average  time*    11.1687090r           std  dev=  3.24415240 

average  vel  =     5.5843543-            std  dev=  6.72964110 

set  nc=    2         53 -da . a  .  i  nts 

width  =  o50    inches  mercury*  2,50 

v=     .'63094580  pr=    8*42353800  sgima  sg=  3.35534000 

average  time*    5-31795200           std  dev*  1. 83175870 

average  vel=     2.6589760o           std  dev=  .3290iai6 

SET  NC=     3  36. DATA  POINTS 

WIDTH  =   .50     INCHES  MERCURY*   2.50 
V=      4.15123330   PR=      2.83539000   SGIMA  SC=     48.86154000 

average  times    11.77036700  std  dev=    6.99cio3oo 

AVERAGE  VEL=       5.83518350        ■        STD  DEV=      5.34C071C0 

SET  NC=     4  42. DATA  POINTS 

WIDTH  =  .50     INCHES  MERCURY=   2.50 

V=  .91904040   PR=      5.34654300   SGIMA  SG=      4.51583000 

average  time*    4.9i369ooo  std  dev=    2.12505990 

average  vel*     2.45634500  std  dev*     .49369549 

set  nc=    5         38. data  pcin's 

width  =  1.00    inches  mercury=  205g 

v=    1.32273230  pr=    4.21963300  sgima  sq=  7.33275900 

average  time*    5.58144600  std  dev*    2.71712320 

average  vel=     5.58144600  std  dev=    1.21371300 

set  nc=    6         36. data  points 

width  =  1.00    inches  mercury=  2o50 

v=     .56442770  pr=    5.24892170  sgima  sc=    1.67219450 

average  time*    2.96263700  std  dev=    1.29313350 

average  vel=     2.96263700  std  dev=     .23265245 

set  nc=    7        33.data  points 

width  =  1.00    inches  mercury=  2.50 

v=    1.29299350  pr=    4.82538400  sgima  so*    8.06729600 

average  time*    6.23921500  std  dev=    2.34029350 

average  ve_  =     6.23921500  std  dev=    1.42610990 

set  n0=    8         34. data  points 

width  =  1.00    inches  mercury*  2.50 

v=     .39266850  pr=    8.79659700  sgima  sq=    1.35633480 

average  time*    3.45414700  std  dev=    1.16461730 

average  »/el=     3.45414700  std  dev*     .23610753 

set  n0=    9         30.data  points 

width  =  1.00    inches  mercury*  2.50 

v=     .88564100  pr=    3.76749400  sgima  sq=    2.95507240 

average  time*    3.33664800  std  dev=    1.71903240 

average  vel*     3.33664800  std  dev*     .54374316 

set  no*   10         35. data  points 

width  =  1.00    inches  mercury*  2.50 

v*     .49079219  pr*    3.65043500  sgima  sq=     .87930590 

average  time*    1.79160530  std  dev*     .93771312 

average  vel*     i.7916C530  std  dev*     .15544073 


Logic  Diagram  for  Computer  Program  2To.    7 


<LL 


Calcula  to 
3AD  = 
(B-BA)     +  3AD 


Calculate 

P.N    = 

(SA)2     EDS 

' 

' 

Calculate 

V  = 

bds/b 

A 

:nd 


Ccicalate 
3DS  = 
3AD/P-I.0 


Cg!cuIq:c 
TP  = 
(PR-l.O)V 
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jjORTRAN  Source  Fro  gram  for  Computer  Program  No.  7 
C      PROGRAM  FOR  MU»SIGMA  SQUARED ,P ,NU > AND  TP 
DIMENSION  D(300) 

101  FORMAT  (14,  I4»  E7.4) 

102  FORMAT  (15) 

103  FORMAT ( 14 » F15.8 ,F15 *  8 ) 
I  =1 

200  READ  102 »  JMAX 
J=  1 

201  READ  101,  M,  N,  D(J) 
J  =  J  +  1 

IF  (JMAX-J)  202?  201,  201 

202  J  =  1 

P  =  1.0 
B  =  D(J) 

1  J  =  J  +  1 

IF  (JMAX-J)  3,  2»  2 

2  P  =  P  +  1.0 
B=B+D(J ) 

GO  TO  1 

4,5. 


3 

BA=3/P 

IF  (SENSE  SWITCH  3) 

4 

PRINT  103»M»P»BA 

GO  TO  6 

5 

PUNCH  103»M»P»BA 

6 

J  =  l 

P  =  1.0 

SAD=0.0 

7 

B=D(J) 

SAD=(B-3A)**2-i-BAD 

J=J  +  1 

IF  (JMAX-J)  9,  8,2 

3 

P«P+1.0 

GO  TO  7 

9  BDS=BAD/(P-1.0) 

IF  (SENSE  SWITCH  3)  10,  11 

10  PRINT  103»M»P»BDS 
GO  TO  12 

11  PUNCH  103»M»P»BDS 

12  PR=(BA)**2/(BDS) 
V=(BDS)/(BA) 

IF  (SENSE  SWITCH  3)  13,  14 

13  PRINT  103»M»PR»V 
GO  TO  15 

14  PUNCH  103>M,PR»V 

15  TP=(PR-1.0)*V 

IF  (SENSE  SWITCH  3)  16,  17 

16  PRINT  103,M,P,TP 
GO  TO  18 

17  PUNCH  103,MvP,TP 

18  1=1+1 

GO  TO  200 
END 
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K)RTRAN  Source  Program  for  Coraputer  Program  No.    8 


DIMENSION  D(150) 

101  FORMAT ( I4»I4»F7.3) 

102  FORMAT (  15  ) 

103  FORMAT(E10a4»E10«4> . 

1  FORMAT (7HSET  N0=I 5 > 10XF4.0 > 1 1HDATA  POINTS) 

2  FORMAT (7HWIDTH  =  F5 • 2 »4X15HINCHES  MERCURY=F6 .2 ) 

3  F0RMAT(2HV=F15.8»5H   PR=F15.8>11H   S6IMA  SQ=F15.8) 
t4  FORMAT (13HAVERAGE  TIME=F15.8» 14X8HSTD  DEV=F15.8) 

*5  FORMAT (13HAVERA6E  VEL=  F15 .8 > 14X8HSTD  DEV=F15.8> 

READ  102 » K 
105  I  3AD  103»CL»VELW 

READ  102»JMAX 

DO  700  J=1,JMAX 
700  READ  101,M,N»D(J) 

SUM1=0«0 

SUM=0.0 

DO  800  J=1,JMAX 

D( J)=1«0/D( J) 

SUM=SUM+D( J) 
800  SUM1=SUM1+D( J)**2 

CN=JMAX 

SGSQ= ( SUM  1- ( SUM**2  > /CN ) / ( CN-1 * 0 ) 

STDD1=SQRT(S6SQ) 

AVG=SUM/CN 

PR=<AVG**2)/SGSQ 

V=SGSQ/AVG 

VEL=CL*AVG 

STDD2=SGRT { CL* ( SGSQ**2 J / ( CN-1 • 0 ) ) 

PUNCH  1»K»CN 

PUNCH  2  >CL»VELW 

PUNCH  3»V*PR»SGSQ 

PUNCH  4>AVG»STDD1 

PUNCH  5>VEL»STDD2 

K  =  K+1 

GO  TO  105 

END 
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APPENDIX  G 


Description  and  Explanation  of  the  IBM 
lo20  Computer  FORTRAN  Program  No.  9  for 
Calculating  the  T  Statistic  for  Given 
Gamma  Distribution  Parameters  p  and  ~0 


This  program  determines  T,  the  upper  limit  in  the  integral 
,T 
f(t)dt  =  A 

where  A  is  the  area  under  the  distribution  f(t)  integrated  from  0  to  T. 
Table  2k   defines  the  various  symbols  utilized  in  this  computer  code. 


Table  2k.      Input  data  and  variables  required  for 
the  IBM  1620  Computer  Program  No.  9. 


SYMBOL 


EXPLANATION 


KMAX 
A(K) 

X(K) 

V 

N 
T 

GAMAP 
Y 


Order  of  Gaussian  quadrature  used  in  the 

program 

Christoffel  Numbers  for  the  order  of  Gaussian 

quadrature  used  in  the  program 

Zeros  of  the  Legendre  polynomial  of  degree  KMAX 

Parameter  v    of  the  gamma  distribution 

Parameter  p  of  the  gamma  distribution 

Upper  time  limit  of  the  integral 

The  gamma  function  P(p) 

The  value  of  the  integral  for  given  V,  N,  and  T 


Input  to  Computer  Program  No.  9 


The  input  to  Computer  Program  No.  9  is  given  in  Table  25-   The  order  of 
Gaussian  quadrature  KMAX  is  read  into  the  computer  as  FORMAT  (15)-   The 
Christoffel  numbers  A(K)  are  read  in  as  FORMAT  (EI5.8).   The  zeros  of  the 
Legendre  polynomial  of  degree  KMAX  is  read  in  as  FORMAT  (EI5.8).   The 


Table     25.         Input   to  Corn   iter  Program  Ho.    9 


KMAX 
5 


A(K) 

. +7S62S67 
•236926SS 
•56SS8689 
.23692688 
.47862867 


X(K) 
.53846931 
.90617935 
.00000000 
.90617985 
.53346931 


V 

IT 

T 

.01 

1 

0.01 

.01 

1 

0*02 

.01 

1 

0.03 

•  05 

1 

0.05 

.05 

1 

0.10 

•  05 

1 

0.20 

•  10 

1 

0.10 

•  10 

1 

0«20 

•  10 

1 

0*30 

•  20 

1 

0*40 

•  2u 

1 

0»50 

•  20 

1 

0*60 

•  30 

1 

0«60 

•  30 

1 

0.70 

•  30 

1 

0-80 

.40 

1 

0«80 

.40 

1 

0.90 

.40 

1 

1.10 

•  50 

1 

1.00 

.50 

1 

1.10 

.50 

1 

1.30 

GAKAP 
1.0 
1.0 
1.0 
1.0 
1.0 
1.0 
1*0 
1*0 
1  •  0 
1.0 
1.0 
1*0 
1*0 
1.0 
1*0 
1*0 
1.0 
1*0 
1.0 
1.0 
1.0 
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parameter  ~V  of  the  gamma  distribution  (V),  the  parameter  p  of  the  gamma 
distribution  (N),  the  T  statistic,  and  the  gamma  function  GAMAP  are  read 
into  the  computer  as  FORMAT  (E10.4,  15,  E10.4,  EI5.8). 

Output  of  Computer  Program  No.  9 

The  output  of  Computer  Program  No.  9  is  given  in  Table  26,   where  Y 
is  the  value  of  the  integral 

T 


I 


f(t)dt  -  y 

for  given  V,  N,  and  T. 
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Table  26.  Output   of  Computer  Program  ITo.    9 

V  H  T  Y 

•63212045E-00 
,3(  466< .>GE-0  0 
W9502I280E-G0 
•63212045E-00 
•86466460E-00 
-931633S6E-00 
.63212045E-00 
•86466460E-00 
<>9502I280E-0G 
•86466460E-CG 
.91791491E-0G 
«95021230E-00 
•86466460E-00 
<,90302786E-00 
•93051638E-00 
»86466460E-00 
.89460063E-00 
•93607194E-00 
•  86466-460E-00 
o83919667E-00 
«92572626E-GO 
»86466460E-0G 
.88544102E-00 
•95021280E-00 
.86466460E-00 
•88268068E-0O 
•94256696E-00 
•86466460E-00 
o88C56690E-00 
«956G6291E-00 
•86466460E-00 
•87889649E-00 
.96432545E-00 
e J6466460E-00 

.91791491E-00 

«95021280E-00 
*8G0S5175E-00 
*9G84234SE-00 
.95958613E-00 
•59399403E-00 
a90342348E-00 
•93271180E-00 
o80085I75E-00 
.90842348E-00 
•95958613E-00 
♦90842348E-0G 


„1000E-01 

1 

•1000E-0I 

.IGOGE-Gl 

I 

♦2GG0E-01 

.lOGOE-01 

1 

.3000E-01 

.50C0E-01 

1 

•5C  iOE-Ol 

.5000E-01 

1 

-1GOGE- jO 

.5000E-01 

1 
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Lo^ic  Diagram  for  Computer  Program  No.   9 
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Start 

Read 

KMAX,   Ak,Xk 


y 


Form 
EX  = 
Exp(XI) 


END 


Road 

V,  P,  T,  T  (P) 


Calculate 
X  I  = 


ro)(xk  +  0 


Calculate 

C  A  =  CA  + 

EX-  A,.  (X..+  I) 


P-! 


Punch 

KMAX,   Ak,  Xk 


3c! 
-s-icA  =  0.0 


Set 
K  =  I 


Sot 

:<  =  k  +  I 


Calculate 


Calculate 
Y  =  PRE-CA 


Punch 

V,  P,  T,Y 
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FORTRAN  Source  Program  for  Computer  Program  ITo.  9 

C      PROGRAM  FOR  GAUSSIAN  QUADRATURE 
DIMENSION  A(6) »X(6) 
101  FORMAT (£15. 3) 

103  FORMAT ( 15 ) 

104  FORMAT  <E10.4»I5»E10«4»E15*8) 

10  5  F ORMAT ( E10.4 » E10.4 ,E10.4, E 10. 4 > ElO *4» E10 .4 ) 

200  READ  103*KMAX 
PRINT  1G3,KMAX 

;<  =  i 

201  READ    1015AC<) 
PRINT    101, A  CO 
K  =  K+1 

I F ( KMAX-K ) 20  2  » 2  0 1,201 

202  K  =  l 

203  READ101,X(K) 
PRINT    101,X(K) 
K=K+1 

I F ( KMAX-K ) 204 , 2  03 , 203 

204  READ  104»V,N,T»GAMAP 
CA=0.0 

PRE=( (T/2«0)**N)/( (GAMAP)*(V>**N) 

1  K  =  l 

2  XI=-(T/(2.0*V))*(X(K)+1.0) 
EX=EXP(XI ) 

CA=(A(K)*<  (X(K)+1.0)**(N-1)  )*EX)+'CA 
K=K+1 

IF(KMAX-K)3»2,2 
3   Y=PRE*CA 

IF (SENSE  SWITCH  3)4,5 
4  PRINT  104,V»N,T,Y 
GO  TO  204 
5    PUNCH  104,V,N,T»Y 
GO  TO  204 
END 
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APPENDIX  H 


Discussion  of  the  C  Russian  Quadrature 
Approximation  of  the  Gamma  Distribution 


Certain  required  integrations  are  of  such  complexity  that  direct  inte- 
grations are  not  feasible.   It  is,  therefore,  necessary  to  fit  the  data  with 
empirical  expressions  and  solve  them  by  formal  analytical  methods.   Such  a 
technique  is  commonly  known  as  numerical  integration.   The  general  procedure 
tabulates  values  of  the  ordinate  y  at  n  predetermined  values  of  x.   The 
integral  over  the  region  is  then  the  sum  of  the  products  of  the  ordinate 
y  multiplied  by  the  difference  in  the  integration  limits  of  x, 

f(x)dx=  a1y(x1)+a2y(x2)  +  ...+any(xn)+E(x)    (1H) 
Jb 

where  E(x)  is  the  error  resulting  from  the  replacement  of  x  with  a  poly- 
nomial approximation. 

Gauss  utilized  unequal  spacings  of  specific  values  of  x.  to  use  in 
equation  (lH).   The  technique  developed  is  simple  to  use,  is  more  accurate 
than  the  well-known  equal- interval  numerical  integration  formulas,  and  is 
well- suited  for  digital  computation.   The  method  is  generally  known  as 
Gaussian  quadrature . ( 1 , 2 ) . 

Lapidus  (6)  presents  an  excellent  discussion  of  the  derivation  of  the 
technique  of  Gaussian  quadrature  numerical  integration.  Highlights  of  his 
discussion  are  as  follows: 

Certain  properties  of  orthogonal  functions  should  firsi:  be  presented. 
A  pair  of  real  functions  is  said  to  be  orthogonal  over  a  region  with  respect 
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to  a  weighting  function  w(x)  if 

J  v(x)gn(x)^n(x)dx  =  0  ,   -,  :n;  (2K) 

a  =  C,,  ,  n  =  m  . 

By  proper  choice  of  constants  involved,  in  g  (x)  and  w(x), 

Cn  -  1.0.  (5H) 

Functions  of  the  first  kind  and  degree  n  which  are  solutions  of 
Legendre 's  equation 

(l-x2)y"-2xy,+n(n+l)y  =  0  (4H) 

are  known  as  the  Legendre  polynomials,  indicated  symbolically  as  P  (x). 
The  Legendre  polynomials  are  orthogonal  in  the  interval  (  -1,  l)  when 
w(x)  ~  1.0. 

I  Pn(x)Pm(x)dx  =  0  ,       n^n;  (5H) 

"'  =  2/(2n+l)  ,  n  =  m;  n,m  =  0,1,2, .. . 

Pn(x)  is  orthogonal  to  any  polynomial  of  degree  less  than  n,  since 
Pn(x)  is  a  polynomial. 

Pn(x):crndx  =  0  ,      n  =  0,1,2,. . . ,n-l  .      (6H) 

rb 

Going  back  to  Equation  1H,   \  f(x)dx   will  be  exact  if  f(x)  is  a 


degree  of  2n-l  or  less.   This  is  because  there  are  2n  coefficients,  a^  and 

x^,  to  be  determined. 

Choosing  f(x)  to  have  the  form 

f(x)  =  xrP1T(x)  ,     r  =  0,1,...,1T-1,  (7H) 

P^j(x)  is  the  Legendre  polynomial  of"  degree  N. 
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If  Pjj(xk)  —  0  for  k  s  1,  2,  ...,  N;  ^  must  "be  the  roots  of  the 
Legendre  polynomial  of  degree  N. 

Choose  the  polynomial  f(x)  to  he  ,  at  x  s  x^,  Xo,  ..*,   Xjt;  except  at 
x,  where  it  has  the  value  1.0. 

One  polynomial  of  degree  N-l  which  fits  the  N  points  is 

fW-    X  Pl[(x)y(Xk)  (8H) 

If,  however,  f(x)  is  a  polynomial  of  degree  2N-1  or  less,  it  is  possible 
to  write 

f(x)  =   X     P]T(x)y(Xk)    ♦  %(x)^         (9H) 

A=/    (x-Xv)P^(xk) 
Integrating  both  sides  between  a  and  b, 


f(x)dx=  V  aky(xk)  +  E(x)  (10H) 


Thus; 


ak^   (  P1-(x)dx/(x-xk)P1![(xk)  (11H) 


M 


E(x)  =    f  P,T(x)xindx  (12H) 


If  E(x)   ==    0,   then  this  representation  will  be  exact  for  f(x)  being  a  poly- 
nomial of  degree  2N-1  or  less.      Recalling  that  if  (   a,  b   )   =    (    -1,   l), 


E(x)  =       \   PIJ(x)xmdx  =  0    ,     m  =■  1,2,3, . . .  ,1T-1 


225 


Therefore 


f(x)dx  =  ^  aky(xk)  .  (13H) 

This  also  implies  that  the  x^'s  in  Equation  10  H  are  the  roots  of  the 
Legendre  polynomial  of  degree  IT.   Choice  of  a  value  N  automatically  deter- 
mines X-,,  yj^j    ...,   x^,   and  the  a^'s  are  given  uniquely  "by  Equation  11H. 
While  only  N  ordinate  values  are  used,  the  representation  behaves  as  if  217 
ordinate  values  are  used  as  long  as  the  x^'s  are  the  roots  of  the  Legendre 
polynomial. 

Any  finite  interval  of  integration  (  a,  b)  may  be  changed  to  (  -1,  l) 
by  means  of  the  transformation 

x  =    2u-(a+b)  (14E) 

"b-a 


"b-a 


f(u)du=  —    \f(x)dx  (15K) 

2   J-, 

Table  27  gives  x.  and  a,  values  for  values  of  n  up  through  n  =  6  ( 2  ) . 

Since  x^  =  -x^^.^  and-  ak  —  ^-k+l'  only  fo-lf  of  "the  values  are  tabulated. 
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Table  27«  Parameters  utilized  in  Gauss'  mechanical 
quadrature  formula  ( 2  ) . 


n 


k 


av 


3 


1 

2 

1 
2 

1 

2 

3 

1 
2 
3 


O.577350269189626 

0.77^5966692^1^83 
0.000000000000000 

0.3399810^-358^856 
0.86113631159^53 

0.538^9310105683 
0.906l798^593866i+ 
0.000000000000000 

0 . 238619186083197 
0.3607615730^139 
0.932^6951^203152 


l.OOOOOOOCGOOOGO 

0.55555555555556 
0.88888888888889 

0.6521^515^6255 

0.3W5143^5137^5 

0. ^7862867049937 
0.23692688505619 
0.56888888888889 

0.4o79l3931+57269 
0.3607615730^81^ 
0.17132^9237917 
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ABSTRACT 

The  movement  of  a  single  particle  in  a  liquid  fluidized  "bed  was 
followed  using  a  radioactive  tracer  technique.   The  fluidized  bed  was 
made  up  of  precision  ground  stainless  steel  balls  fluidized  with  water. 
A  single  ball  bearing  was  made  radioactive  by  irradiation  in  the  Kansas 
State  University  TRIGA  Mark  II  Nuclear  Reactor,  and  introduced  into  the 
bed.   The  axial  motion  of  this  particle  was  followed  by  measuring  the 
passage  of  the  particle  past  two  horizontal  planes  in  the  bed. 

The  construction  of  the  fluidized  bed  and  the  detection  system  that 
was  developed  are  described  in  detail. 

Preliminary  data  consisting  of  over  20,000  data  values  were  obtained 
for  a  single  sized  particle  and  single  packed  bed  height.   The  data  were 
analyzed  to  give  several  parameters  that  can  be  used  to  characterize 
a  fluidized  bed.   These  parameters  include  the  average  particle  velocity, 
particle  velocity  distribution,  reflection  coefficient,  particle  traffic, 
residence  time  distribution,  dispersion  coefficient,  reflectivity  and 
transmissivity.   The  mathematical  analyses  necessary  to  describe  these 
parameters  are  presented.   IBM-1620  computer  programs  in  FORTRAN"  used  to 
reduce  the  data  to  obtain  these  parameters  are  given. 
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and 


=  -%  .  (97) 


z 


An  iterative  procedure  used  to  solve  for  E2  was  as  follows: 

The  quantities  (texp)m,L,  and  (x-x0)m  were  extracted  from  the  data  first. 
A  value  of  E  was  then  assumed.   The  quantities  ^m  and  tm  were  calculated, 
where 


~      (texp)mJJiz 

Cm    ~    S > 

1/ 


and 

t^  =  time  corresponding  to  (G/C«>)nax  • 
The  left-hand  side  of  Equation  °A  was  calculated  after  inserting  the  quan- 
tities into  Equations  95;  9&,   and  97  •   The  value  of  the  left-hand  side  of 
Equation  9^-  vas  plotted  versus  Ez,  and  the  Ez  corresponding  to  the  value  of 
the  left-hand  side  of  Equation  9^  being  equal  to  zero  was  the  least  squares 
analysis  solution  for  the  dispersion  coefficient  for  the  given  water  flow 
rate,  based  upon  minimizing  the  sum  of  the  squares  of  the  deviations  of  the 
times  corresponding  to  the  correlated  {0,/C^  )         and  the  experimental 
(-V^co'max  with  respect  to  Ez. 

Method  2.   The  second  method  minimized  the  sum  of  the  squares  of  the 
errors  of  the  experimental  values  of  I/I^  for  five  distances  between  detec- 
tor levels  with  respect  to  Ez  for  each  water  flow  rate  studied.   The  least 
squares  analysis  equation  used  for  the  second  method  was 


>z 


m=i  j=i 


(c/cjjm-  (i/ijjm 


=0,   (98) 


Ik2 


where  in  is  the  subscript  for  the  distance  "between  detector  levels,  AT.,    and 
j  is  the  subscript  for  the  correlated,  C/C^,  ,  and  experimental,  l/lw,  par- 
ticle concentrations  as  functions  of  time  for  each  distance  between  detector 
levels.   The  values  of  I/Iqj  were  obtained  from  Figs.  86  through  90.  From 
Equation  97 > 


J  El    _ 


^m 


z 


Differentiating     (C/Ca  )^m-(l/la} ^J 


with  respect  to    /C!\m  gives 


"3* 


(C/Cjta-  (I/I J 


'Jia 


>f  =  2   (c/cj^-  (lAjJ-^g 


(99) 

where  ..  .,  ,  y.iiF1  is  obtained  by  differentiating  Equation  85  with  respect  to 

fjm-      Thus, 


MC/Cjt1m 


2^2< 


=  2  \      (-n<V)e 


•jm 


cos  (ntfx^)  c  o  s(ntfx    ) 


r\  = 


-(x-Xo^/4^ 


33a 


4if*^r 


2  r,. 


(100) 


An  iterative  procedure  was  used  to  attempt  to  obtain  a  solution  for  Ez.   The 
procedure  was  as  follows: 

The  quantities  (i/I^).  >L,t.  ,x  ,  and  x  were  extracted  from  the  data 
first.  A  value  of  E  was  then  assumed.   The  quantities  ^  and  (C/C^  ) 


jm 


'a  'jm 


were  then  calculated  from  Equations  83  and  85.   The  numerical  value  of  the 
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left-hand  side  of  Equation  98  was  calculated,  utilizing  Equations  97  >  99> 
and  100.   The  value  of  the  left-hand  side  0^  Equation  98  was  then  plotted 
versus  Ez,  and  the  Ez  corresponding  to  the  value  of  the  left-hand  side  of 
Equation  98  being  equal  to  zero  was  to  be  the  dispersion  coefficient  based 
upon  minimizing  the  sum  of  the  squares  of  the  deviations  of  (c/CooJ-Cl/lc©) 
with  respect  to  Ez.   However,  a  plot  of  Ez  versus  the  left-hand  side  of 
Equation  98  gave  curves  asymptotically  approaching  zero  with  increasing  E2 
for  all  water  flow  rates  studied. 
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5.0  DISCUSSION  OF  RESULTS 

This  section  has  been  divided  i.  ,0   fou:     .3.  Each  part  contains  a  dis- 
cussion of  one  of  the  four  quantities  oxtr  cted  frora  the  experimental  data. 

Almost  20,000  individual  data  values  were  obtained  in  this  investigation, 
and  these  are  recorded  on  IBM  cards.   Since  it  is  impractical  to  present  all 
of  the  data  and  results,  only  representative  samples  have  been  shown  in 
Section  4.0  and  in  the  appendices.   A  complete  listing  of  the  data  and  results, 
in  the  form  of  computer  input  and  output,  is  available  in  Project  232^,  Supple- 
ment Number  One,  which  may  be  obtained  from  the  Kansas  State  University 
Engineering  Experiment  Station. 

5.1  Residence  Times 

Figures  35  through  38  illustrate  that  about  50  individual  data  values 
per  data  set  were  sufficient  to  study  the  characteristics  of  the  residence 
time  distribution,  with  adequate  assurance  that  the  samples  taken  represented 
the  population  distribution. 

Figures  39  through  k6   show  that  the  residence  time  data  can  be  repre- 
sented by  a  gamma  distribution  in  practically  all  of  the  data  sets  obtained. 

Figures  k"J   through  56  give,  as  functions  of  the  distance  between  detector 
levels,  the  trends  of  the  gamma  distribution  parameters  p  and  V  for  each  water 
flow  rate  studied.  Knowledge  of  their  trends  enables  prediction  of  the  charac- 
teristics of  the  fluidized  bed,  as  discussed  in  Section  U.3. 

The  residence  time  results  were  corrected  for  instrument  characteristics 
error.  Experimental  residence  times  of  less  than  0.2  second  were  subject  to 
a  major  correction,  and,  therefore,  the  data,  in  cases  where  very  short  resi- 
dence times  were  obtained,  were  least  reliable.  The  major  source  of  this 


